athematical ‘Tables 


and other 


Aids to Computation 





A Quarterly Journal 


Edited by 
.C. CRAIG ALAN FLETCHER EUGENE ISAACSON 
C. V. L. SMITH A. H. TAUB 
C. B. TOMPKINS, Chairman 





XII 


Nos. 61-64 
1958 


Published by 
THE NATIONAL RESEARCH COUNCIL 
Washington, D.C. 














olume XII January 1958 Number 61 


and other 


Aids to Computation 








Published Quarterly 
by the 
National Academy of Sciences—National Research Council 





Published quarterly in January, April, July, and October by the National Academy of Sci 
National Research Council, Prince and Lemon Sts., Lancaster, Pa., and Washington, D. C. 


Entered as second-class matter July 29, 1943, at the post office at Lancaster, Pennsylvania, under 
the Act of August 24, 1912. 











Editorial Committee 
Division of Mathematics 
National Academy of Sciences—National Research Council 


Washington, D. 0. 


C. B. Tompxins, Chairman, University of California, Los Angeles, California 
C. C. Crate, University of Michigan, Ann Arbor, Michigan 
Auan Fietcuer, University of Liverpool, Liverpool 3, England 


EuaenE Isaacson, Institute of Mathematical Sciences, New York University, 
New York 3, New York 


C. V. L. Surry, Ballistic Research Labs., Aberdeen Proving Ground, Aberdeen, 
Maryland 


A. H. Tavs, University of Illinois, Urbana, Illinois 


Information to Subscribers 


The journal is published quarterly in one volume per year with issues numbered 
serially since Volume I, Number 1. Subscriptions are $5.00 per year, single 
copies $1.50. Back issues are available as follows: 


Volume I (1943-1945), Nos. 10 and 12 only are available ; $1.00 per issue. 

Volume II (1946-1947), Nos. 13, 14, 17, 18, 19, and 20 only are available; 
$1.00 per issue 

Volume III (1948-1949), Nos. 21-28 available ; $4.00 per year (four issues), 
$1.25 per issue 


Volume IV (1950 and following), all issues available; $5.00 per year, $1.50 
per issue 


All payments are to be made to the National Academy of Sciences and for- 
warded to the Publications Office, 2101 Constitution Avenue, Washington, D. C. 


Agents for Great Britain and Ireland: Scientific Computing Service, Ltd., 
23 Bedford Square, London W.C.1 


Microcard Edition 


Volumes I-X (1943-1956), Nos. 1-56 are now available on Microcards and may 
be purchased from The Microcard Foundation, Box 2145, Madison 5, Wisconsin, 
at a cost of $20.00 for the complete set. Future volumes will be available on 
Microcards in the year following original publication. 


Information to Contributors 


All contributions intended for publication in Mathematical Tables and Other 
Aids to Computation and all books for review should be addressed to C. B. 
bee nam Department of Mathematics, University of California, Los Angeles 24, 
California. The author should mention the name of an appropriate editor for 
his paper if this is convenient. 








» 2064 
5 ESSA 


Coefficients and Roots of the Polynomials which 
Define the Derivatives of the Exponential 
of (—</T) 


By Edwin S. Campbell, E. M. Fishbach, 
and J. O. Hirschfelder 


The n-th derivative of exp (—«¢/T) can be shown to be of the form 
(1) d* exp [—«/T]/dT* = T exp [—¢/T]W,[¢/T] 


W.le/T] = > bl TY. 


The 6,,: are constant coefficients in a homogeneous polynomial of the m-th order in 
[e/T]. These polynomials attain added importance from their relation to the La- 
guerre polynomials [1 ] 


(2) Wole/T] = (—1)"{L,[e/T] — nL,«[</T]} 


L,[¢/T]: the m-th order Laguerre polynomial. 

The zeros of these polynomials are of interest since they locate the extreme values 
of the derivatives of exp (—«/T). Furthermore, an accurate, simple computation 
of the W, for values of (¢/T) which occur in physical problems sometimes requires 
use of the polynomial roots. This occurs whenever the direct evaluation of W,[«/T] 
by synthetic division introduces at intermediate steps of the calculation numbers 
which are larger than the value of the polynomial. In such cases the subtraction of 
these larger numbers can require the use of considerably more than the usual guard 
figures which allow for the effect of ordinary round-off error. Fortunately the occur- 
rence of this disastrous subtraction in synthetic division does not imply the occur- 
rence when the polynomial is calculated using its roots 


(3) W.L*/T] = II [e/T — rat] 


,1: the I-th root of W,. 


The polynomial coefficients, 5,,1, are integers. Table I lists all of their non-zero 
digits which were computed from the recursion relations 


(4) baa = (—1)""n! 
Dae = — (m+ 1 — 1) dna + On-ren(2 <1 Sm — 1) 
baw = 1. 
The values in the Table were checked numerically by the distinct relations 


(5) (n pa L) bat = —n(n aoe 1)d,—1,2. 
Received 16 September 1957. 
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SVoOnNIAUeP Whe 


COEFFICIENTS AND ROOTS OF POLYNOMIALS 


Taste I. 


Bas 


+ 
w 
NaS 


| 
to 
NPD 


+ 
-_ 
Not 


—4.032 
+1.4112 
—1.4112 


Exact Values of the Coefficients, ba: 


Pau 
0 


0 
0 


CrNNN one ooo 


Or N&wWWwh 


OM WPA a PORK WOUAaAAN COM WP Pan CORN WEP & 


n 
il 


12 


14 


15 


— 
COONAUMEP WN KF OODBNIAMPWNHE ow 


on’ 
Ne oO 


Bat 


+3.9916 8 
—1.9958 4 


—4.7900 16 
+2.6345 088 


+3.2931 36 
—1.3172 544 
+3.0735 936 
—4.3908 48 
+3.9204 
—2.178 


—1.32 


+6.2270 208 
—3.7362 1248 
+6.8497 2288 


+2.5686 4608 
—6.8497 2288 
+1.1416 2048 
—1.2231 648 


+1.3076 74368 
—9.1537 20576 
+1.9833 06124 8 
—1.9833 06124 8 
+1.0908 18368 64 
—3.6360 61228 8 
+7.7915 59776 
—1.1130 79968 
+1.0821 6108 
—7.2144 072 
+3.2792 76 
—9.9372 

+1.911 

—2.1 


SCHRAUANTDWOOOCOCH ONWUAI"SO0O™" & 


id 
CONF UADOOOOCOCSO”! 

















COEFFICIENTS AND ROOTS OF POLYNOMIALS 3 
Taste I. Exact Values of the Coefficients, b,,.-—Continued 
P., P n I Bat Pas n l Bas Pat 
7 16 1 —2.0922 78988 8 13 19 1 +1.2164 51004 08832 17 
Hi 2 +1.5692 09241 6 14 2 1.0948 05903 67948 8 18 
3 3 —3.6614 88230 4 14 3 +3.1019 50060 42521 6 18 
; 4 +3.9666 12249 6 14 4 4.1359 33413 90028 8 18 
7 5 —2.3799 67349 76 14 5 +3.1019 50060 42521 6 18 
7 6 +8.7265 46949 12 13 6 —1.4475 76694 86510 08 18 
; 7. —2.0777 49273 6 13 7 +4.4805 94531 72531 2 17 
; 8 +3.3392 39904 12 8 9.6012 73996 55424 16 
; 9 —3.7102 6656 il 9 +1.4668 61305 02912 16 
> 10 +2.8857 6288 10 10 1.6298 45894 4768 15 
3 11 —1.5740 5248 9 11 +1.3335 10277 2992 14 
12 +5.9623 2 7 12 —8.0818 80468 48 12 
8 13 —1.5288 6 13 +3.6264 84825 6 11 
9 14 +2.52 4 14 —1.1955 44448 10 
9 15 —2.4 2 15 +2.8465 344 8 
9 16 +1 0 16 —4.7442 24 6 
9 17 +5.2326 4 
8 18 —3.42 2 
- 17 1 +3.5568 74280 96 14 19 +41 0 
6 2 —2.8454 99424 768 15 
5 3 +7.1137 48561 92 15 20 1 —2.4329 02008 17664 18 
3 4 —8.2993 73322 24 15 2 +2.3112 56907 76780 8 19 
3 5 +5.3945 92659 456 15 3  —6.9337 70723 30342 4 19 
0 6  —2.1578 37063 7824 15 4 +9.8228 41858 01318 4 19 
7 +5.6514 78024 192 14 5 —T7.8582 73486 41054 72 19 
9 8 —1.0091 92504 32 14 6 +3.9291 36743 20527 36 19 
10 9 +1.2614 90630 4 13 7 1.3097 12247 73509 12 19 
10 10 —1.1213 25004 8 12 8  +3.0404 03432 24217 6 18 
10 11 +7.1357 04576 10 9 —5.0673 39053 73696 17 
10 12 —3.2435 0208 9 10 +6.1934 14399 01184 16 
9 13 +1.0395 84 8 11 —5.6303 76726 3744 15 
9 14 —2.2848 6 12 +3.8388 93222 528 14 
8 15 +3.264 4 13 —1.9686 63191 04 13 
6 16 —2.72 2 14 +7.5717 81504 il 
5 17. +1 0 15 —2.1633 66144 10 
4 16 +4.5070 128 8 
2 17 —6.6279 6 6 
0 18 1 —6.4023 73705 728 15 18 +6.498 4 
2 +5.4420 17649 8688 16 19 —3.8 2 
10 3 1.4512 04706 63168 17 20 +41 0 
il 4 +1.8140 05883 2896 17 
12 5 1.2698 04118 30272 17 21 +1 +5.1090 94217 17094 4 19 
12 6 +5.5024 84512 64512 16 2 —5.1090 94217 17094 4 20 
11 7 1.5721 38432 18432 16 3 +1.6178 79835 43746 56 21 
11 8  +3.0881 29063 2192 15 4 2.4268 19753 15619 84 21 
10 9 4.2890 68143 36 14 5  +2.0627 96790 18276 864 21 
9 10 +4.2890 68143 36 13 6 —1.1001 58288 09747 6608 21 
8 11 —3.1193 22286 08 12 7 +3.9291 36743 20527 36 20 
7 12 +1.6541 86060 8 11 8  —9.8228 41858 01318 4 19 
5 13 6.3622 5408 9 9  +41.7735 68668 80793 6 19 
4 14 +1.7478 72 8 10 —2.3647 58225 07724 8 18 
2 15 —3.3292 8 6 11 +2.3647 58225 07724 8 17 
0 16 +4.1616 4 12 —1.7914 83503 8464 16 
17 —3.06 2 13. +1.0335 48175 296 15 
12 18 «+41 0 14 —4.5430 68902 4 13 
12 15 +1.5143 56300 8 i2 
13 16 —3.7858 90752 10 
13 17 +6.9593 58 8 
13 18 —9.0972 6 
12 19 +7.98 4 
. 2” —4.2 : 
1 1 
10 + 
~ 
7 
5 
4 
2 
0 








22 


23 


24 


COEFFICIENTS AND ROOTS OF POLYNOMIALS 


TaBLe I. 


~ 


COnNAUPWHe 


Exact Values of the Coefficients, b, .—Continued 


-1240 


-5852 


.8271 


.2780 
-5436 
. 1879 


-6846 
-4750 
.3216 
-1072 
-4230 
-0633 
-5545 


- 1688 
06 


. 2044 
-1351 
-6162 
.5783 
.5783 
-8996 
2427 
7725 
3833 
.3972 
7782 
7513 
.3471 
- 1423 
8772 
-4539 
.2764 


Bar 


00727 77760 768 
00764 16648 8064 
02547 22162 688 
37366 43424 256 
53629 79081 8304 
07056 88146 37056 
74116 90722 42688 
37813 14435 072 
90192 22512 64 
19166 54740 48 
84545 32444 16 
53787 77036 8 
47299 8528 

22291 136 

$1585 28 

35879 04 


01673 88849 7664 
21841 27734 74304 
26444 47071 60064 
37740 74511 93344 
95853 70786 33676 8 
75122 24718 02060 8 
75644 71909 67500 8 
78755 63402 76428 8 
80740 90422 4256 
81152 51768 2176 
86816 61180 6208 
60742 37436 928 
06933 72551 68 E 
58754 79424 

80377 6896 

26792 5632 

10127 92 

35319 2 

838 

04 


84017 33239 43936 
56619 93225 35526 4 
24093 97515 96359 68 
92164 45652 93629 44 
92164 45652 93629 44 
48370 82246 85965 312 
06444 63820 08270 848 
01706 93739 53679 36 
37126 52754 52620 8 
22854 42125 75436 8 
83632 89978 23283 2 
39941 49220 98688 
84570 37862 2976 
24326 46420 48 

97298 31628 8 

86486 86860 8 

30843 7312 


+9.7826 84936 64 
—1.7162 60515 2 
+2.2582 3752 

















COEFFICIENTS AND ROOTS OF POLYNOMIALS 5 
Taste I. Exact Values of the Coefficients, b,;—Continued 
n l Baw Pas 
24 21 —2.1507 024 7 
22 +1.3965 6 5 
23 —5.52 2 
24 +1 0 
25 1 +1.5511 21004 33309 85984 25 
2 —1.8613 45205 19971 83180 8 26 
3 +7.1351 56619 93225 35526 4 26 
4 —1.3081 12046 98757 98179 84 27 
5 +1.3735 17649 33695 88088 832 27 
6 —9.1567 84328 91305 87258 88 26 
7 +4.1423 54815 46066 94236 16 26 
8 —1.3314 71190 68378 66004 48 26 
9 +3.1437 51422 44782 94732 8 25 
10 —5.5888 91417 68503 01747 2 24 
11 +7.6212 15569 57049 56928 23 
12 —8.0831 07422 27173 7856 22 
13 +6.7359 22851 89311 488 21 
14 —4.4412 67814 43502 08 20 
15 +2.3263 78378 98977 28 19 
16 —9.6932 43245 79072 17 
17 +3.2073 23132 7984 16 
. 18 —8.3851 $8517 12 14 
19 +1.7162 60515 2 13 
20 —2.7098 85024 il 
21 +3.2260 536 9 
22 —2.7931 2 7 
23 +1.656 5 
v3 —6 2 
25 +1 0 
: 26 1 —4.0329 14611 26605 63558 4 26 
2 +5.0411 43264 08257 04448 27 
3 —2.0164 57305 63302 81779 2 28 
+ +3.8648 76502 46330 40076 8 28 
5 —4.2513 64152 70963 44084 48 28 
6 +2.9759 54906 89674 40859 136 28 
7 —1.4171 21384 23654 48028 16 28 
8 +4.8080 90410 80256 27238 4 27 
9 —1.2020 22602 70064 06809 6 27 
10 +2.2704 87138 43454 35084 8 26 
11 —3.3025 26746 81388 14668 8 25 
12 +3.7528 71303 19759 2576 24 
13 —3.3679 61425 94655 744 23 
14 +2.4056 86732 81896 96 22 
15 —1.3746 78133 03941 12 21 
16 +6.3006 08109 76396 8 19 
17 —2.3164 00040 3544 18 
18 +6.8129 41295 16 16 
19 —1.5936 70478 4 15 
20 +2.9357 08776 13 
21 —4.1938 6968 11 
22 +4.5388 2 9 
23 —3.588 7 
24 +1.95 5 
25 —6.5 2 
26 +1 0 
27 1 +1.0888 86945 04183 52160 768 28 
2 —1.4155 53028 55438 57808 9984 29 
3 +5.8981 37618 97660 74204 16 29 
a —1.1796 27523 79532 14840 832 30 
5 +1.3565 71652 36461 97066 9568 30 
6 —9.9481 92117 34054 45157 6832 29 
7 +4.9740 96058 67027 22578 8416 29 
8 —1.7764 62878 09652 58063 872 29 
9 +4.6878 88150 53249 86557 44 28 





COEFFICIENTS AND ROOTS OF POLYNOMIALS 


Tasie I. Exact Values of the Coefficients, b, x—Continued 
Bat 


—9.3757 76301 06499 73114 88 
+1.4489 83610 16459 04935 936 
—1.7563 43769 89647 33255 68 
+1.6887 92086 43891 66592 
—1.2990 70835 72224 3584 
+8.0418 67078 28055 552 
—4.0209 33539 14027 776 
-6261 12828 32878 88 
-3140 94210 2248 
45844 796 
96515 36 
27525 6 
0328 
4 


83446 11713 86050 1504 
92652 25813 71167 70304 
96815 97852 60839 33798 4 
26699 95526 26748 6208 
92039 94631 52098 34496 
60563 95884 16608 73113 6 
74581 12129 80128 38297 6 
29679 20486 75481 43616 
86022 00135 20967 0656 
26046 44729 88708 2496 
15348 73194 36068 0448 
24312 76083 64633 088 
12115 65213 99623 68 
82512 90011 53536 

55008 60007 69024 

88129 65837 49888 - 
82998 37869 4848 

55222 92994 88 

94509 62864 

71207 0152 

73727 648 

54132 8 

7328 

7 

8 

6 


61993 73970 19545 43616 
46679 12355 82736 36106 24 
10056 05601 22313 62478 08 
00512 14546 93167 95203 584 
25640 18183 66459 94004 48 
00512 14546 93167 95203 584 
17090 32042 72110 21352 96 
78142 62588 21186 15531 52 
61249 32548 95126 28633 6 
02499 85010 87805 84140 8 
31590 65018 78937 36243 2 
61259 97752 89460 03968 
92398 69346 10309 0176 
43866 98326 24447 488 

59904 98804 46033 92 

68277 90969 26853 12 

11622 36250 18048 

12400 92647 1296 

19418 18455 68 





COEFFICIENTS AND ROOTS OF POLYNOMIALS 


TaBLe I. Exact Values of the Coefficients, b,,.-—Continued 
Bat 


—2.5101 36688 99593 6 
+5.3788 64333 56272 
—9.3140 50794 048 
+1.2885 05050 56 
1.4005 48968 

-1671 2414 

-1823 024 

-0693 6 

12 


-6525 28598 12191 05863 63084 8 
-8461 66467 27677 03502 26472 96 
-7948 77684 72915 94967 72354 048 
-0384 74790 64060 88677 37796 608 
-2500 17227 83279 15280 59135 5904 
.3750 14356 52732 62733 82612 992 
.5000 08203 72990 07276 47207 424 
.0267 89083 67478 06559 97960 192 
-1374 11089 00627 42266 60433 92 
.3206 25874 34797 31955 41012 48 
-3310 22886 24508 60355 52911 36 
-9158 66275 65580 56572 35251 2 
-2106 14933 44900 65275 79136 
.0648 60102 67214 89543 3216 
.5732 26744 89306 58699 6736 
.8326 67155 58166 16872 96 

-0609 31624 19644 35155 2 

.1500 68990 67169 16896 

.5441 01721 65505 856 

.1838 18919 42646 432 

.1995 68855 77729 6 

-0129 03023 85272 


-9011 44 
Se 


Legend: Wa(x) = 2 bes x! 


bat = Bayz X 10P.!, 1 S Bay < 10 


TaBLeE II. Roots of W, (¢/T) 


Tat 
2.00000 00000 0000 


1.26794 91924 3112 
4.73205 08075 6888 


0.93582 22275 2409 
3.30540 72893 3227 
7.75877 04831 4363 
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COEFFICIENTS AND ROOTS OF POLYNOMIALS 


PWN CONVAMPWNHHE BATANPWNHH AUNEPWHeH UARPWHHE PWNHeE wo 


CSCOCNIAU PRP wWHHE OCONan 


TaBLe II. 


_ 


Ne oe — 
oOuUnreF Oo FKP OMAK OO SK Awe OS PORNO COUN © 


Noe 
m ~IdS 


1 
4 


9 





Tal 


- 74329 19279 8143 
-57163 50076 4627 
.73117 87516 8905 
-95389 43126 8321 


-61703 08532 7825 
-11296 59585 7838 
-61083 31510 1760 
.39906 69712 0486 
-26010 30659 2081 


.52766 81217 1117 
- 79629 98096 4345 
.87664 15204 7699 
.91881 65667 0471 
. 23461 04290 8311 
-64596 35523 8068 


-46102 42198 0496 
.56358 61896 5431 
-35205 05025 3674 
.91629 72490 2042 
-42069 93830 2156 
. 19416 55480 0748 


09217 69079 5447 


-40938 35732 0319 


38496 31848 0312 


-95625 45561 6887 
- 18194 31010 4007 
-16170 96881 4582 
-07005 51268 3715 
.24973 55261 4898 
.58595 52436 5281 


-36817 84529 4174 
- 24335 79621 4047 
-64603 38413 8420 
-61688 25146 3485 
7. 
10. 
14. 
20. 
28. 


0. 


22178 65393 9663 
56732 08077 4184 
83591 45152 6107 
38218 19854 4899 
11834 33810 4993 


33452 86763 2476 


1.12825 33558 7666 
2.39586 99247 4731 
4. 

6.48735 30313 8081 
-42835 48133 3561 
13. 
17. 
23. 
31. 


16684 09879 2878 


10172 35803 6780 
69648 75668 4621 
57778 70883 6019 
68280 09748 3192 
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Roots of W,, («/T)—Continued 


Cc 








COEFFICIENTS AND ROOTS OF POLYNOMIALS 


TaBLe II. Roots of W, («/T)—Continued 
TaJl c 


-30652 67021 3005 
-03279 73987 7972 
. 18961 19419 6840 
- 79904 76060 5497 
-89491 11715 0275 
-52729 20080 7027 
11.77100 65654 9068 
15.74226 02870 2615 
20.63580 56686 4129 
26.82634 99493 7082 
35.27439 07009 6510 


- 28285 83482 3992 
-95232 60413 6462 
-01649 21385 7771 
-49235 40697 7799 
-40549 10200 1572 
-79281 39404 1242 
. 70738 86889 8909 
-22715 23637 8996 
-47199 66342 2982 
-64178 37524 0122 
-12005 86261 0650 
-88928 43760 9550 


WONIAUMEP WHE KFPOUOSNAUPEPWNEHE &~ 


eRe ee OOOO OO Ree OOOO °O OS 


Soutien 
NFO 


-26258 83981 7108 
.88355 03073 8774 
-86903 38151 9979 
.23241 86994 5120 
-99357 56070 7419 
- 18061 04941 0079 
.83280 82510 2000 
-00562 24002 5865 
-77961 36411 7585 
-27791 17554 4831 
- 70503 40989 2773 
-45278 49657 6721 
.52444 75660 1840 


1 
2 
3 
4 
5 
6 
7 
8 
9 


RR ee OOOCOCOCOCS 


- 24503 30150 9310 
-82408 22200 4784 

. 74187 24007 6449 
-00913 24590 1792 
-64163 40665 9064 
6.66134 14966 2717 
9.09830 27037 7053 
11.99370 34331 8574 
15.40498 80310 5145 
19.41499 28281 6615 
24.14975 79520 2928 
29.81810 51196 7132 
36.81962 42569 9721 
46.17743 00169 8345 


Ree eee OOOO CCOO 








COEFFICIENTS AND ROOTS OF POLYNOMIALS 


TaBLe II. 

l 

1 0 
2 0 
3 1 
4 2 
5 4 
6 6 
7 8 
8 11 
9 14. 
10 17. 
11 22. 
12 27. 
13 > B 
14 40. 
15 49. 
1 0 
2 0 
3 1 
4 2 
5 4 
6 5 
7 7 
8 10 
9 13 
10 16 
11 20 
12 24 
13 30 
14 36 
15 43 
16 53 
1 0 
2 0 
3 1 
4 2 
5 3 
6 5 
7 7 
8 9 
9 12 
10 15 
11 19 
12 23 
13 27 
14 33 
15 39 
16 47 
17 57 





Roots of W, (¢/T)—Continued 


Tat 


-22968 05054 2513 


-77214 
.63105 
.81514 
.33716 
. 21464 
-47116 
13833 
25891 
89205 
12262 
07931 
97497 
21658 
84622 


.21614 
- 72638 
-53359 
.64497 
.07097 
-82585 
-92850 
-40380 
. 28466 
-61517 
-45600 
. 89384 
-05986 
. 17069 
-64036 
.52915 


. 20410 
-68576 
-44719 
-49443 
-83615 
-48411 
-45372 
- 76498 
-44386 
.52434 
-05171 
. 08800 
.72155 
.08585 
-40115 
.08820 
. 22481 





49103 
30990 
45900 
40773 
27645 
39813 
19657 
00216 
34381 
01748 
14990 
35524 
37114 
17085 


03052 
82432 
31603 
09986 
81608 
55151 
41853 
82899 
10707 
32168 
60200 
70253 
29201 
45436 
51841 
11602 


91085 
75894 
86793 
08859 
60319 
54208 
29495 
43667 
10519 
30012 
87593 
07008 
63462 
98280 
37933 
96546 
18319 


7539 
6745 
1225 
3756 
5924 
4671 
5080 
2446 
6953 
2668 
4742 
0162 
8568 
5745 


3945 
5183 
7353 
1195 
8018 
0563 
0666 
5104 
0697 
6805 
3371 
3575 
6784 
3216 
5001 
5465 


7933 
9453 
8049 
6248 
9932 
4047 
8128 
5180 
8395 
1457 
8568 
7174 
7824 
6407 
0861 
0141 
1492 
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1380 
2090 
30900 
45000 
65600 
9540 
14300 
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TABLE II. Roots of W, («/T)—Continued 


CONAUMEANE ~ 


Tn, Cc en, 
0.19334 77686 7901 0 5 
0.64946 18107 9655 0 5 
1.37007 54867 1706 0 5 
2.36027 76257 7508 0 5 
3.62737 58656 5560 0 $ 
5.18117 49562 2015 0 5 
7.03443 06230 6566 0 5 
9.20350 66842 8214 0 5 

11.70932 47926 7879 1 69 
14.57876 11306 9283 1 1060 
17.84675 98274 1882 1 15900 
21.55965 87534 0514 2 12 
25.78070 10050 5240 2 17 
30.59981 44865 1918 2 23 
36.15265 02790 4671 2 33 
42.66288 77797 8564 2 45 
50.55778 33977 1778 2 64 
60.93200 77265 0073 2 93 
0.18366 51730 9616 0 50 
0.61681 63821 2673 0 50 
1.30079 94029 6677 0 500 
2.23994 93990 0886 0 500 
3.44047 19428 8150 0 500 
4.91064 58796 7473 0 50 
6.66116 06120 5821 0 50 
8.70559 71771 2802 0 50 
11.06110 86835 7575 1 612 
13.74939 41200 8341 1 9450 
16.79812 30059 5300 1 14100 
20.24308 49633 0246 2 10 
24.13156 72676 0267 2 15 
28.52794 79713 1414 2 20 
33.52361 76172 8890 2 28 
39.25629 42907 8263 2 39 
45 .95295 06486 3834 2 53 
54.04709 30247 0194 2 73 
64.64971 24378 1605 2 105 
0.17490 67523 8661 0 50 
0.58730 30806 3812 0 50 
1.23822 51018 3420 0 500 
2.13139 62600 7712 0 500 
3.27213 31335 1677 0 500 
4.66749 44658 8837 0 500 
6.32653 61976 7362 0 500 
8.26067 09520 1371 0 50 
10.48416 73812 0829 1 55 
13.01484 87721 5828 1 8470 
15.87508 70127 1499 1 1 26000 
19.09325 19076 0624 1 18 
22.70589 38881 7321 1 26 
26.76117 02293 7952 1 357 


11 
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TABLE II. 
l 
15 31 
16 36 
17 42 
18 49 
19 57 
20 68 
1 0 
2 0 
3 1 
4 2 
5 3 
6 4 
7 6 
8 7 
9 9 
10 12 
11 15 
12 18 
13 21 
14 25 
15 29 
16 34 
17 39 
18 45 
19 52 
20 61 
21 72 
1 0 
2 0 
3 1 
4 1 
5 2 
6 4 
7 5 
8 7 
9 9 
10 11 
11 14 
12 17 
13 20 
14 23 
15 27 
16 32 
17 koe 
18 42. 
19 48. 
20 55. 
21 64. 
22 


Roots of W,, («/T)—Continued 


Tal 


.32451 61370 
-48870 33461 
-39342 27457 
- 26881 38498 
-55442 09713 
-37703 78145 


- 16694 61390 
-56049 08612 
-18142 07583 
.03295 74520 
-11969 20402 
-44769 50426 
.02471 21681 
-86043 775386 
-96688 88709 
-35891 86988 
-05493 02175 
-07788 77154 
-45678 84662 
.22887 63492 
-44311 42999 
- 16593 02974 
-49140 82536 
-56112 99787 
-60828 92435 
-07827 57571 
- 11320 96312 


- 15967 90124 
-53602 44507 
- 12962 03228 
-94327 26194 
-98097 42342 
- 24798 60975 
- 75098 94501 
-49829 71042 
.50013 87408 
- 76904 61305 
-32037 80923 
- 17304 70973 
-35054 65523 
-88244 31429 
-80661 96348 
-17279 15883 


04832 56128 
52855 84501 
75689 07697 
96946 93367 
61735 31062 


0729 
4824 
7640 
6849 
1555 
5231 


1971 
5406 
5136 
5904 
1490 
1758 
1509 
0933 
4385 
9685 
9276 
4042 
6343 
4893 
0910 
6011 
2833 
8368 
3047 
7212 
7984 


6359 
2533 
4296 
9229 
9980 
3645 
8811 
4725 
1857 
3152 
6802 
6204 
0563 
7127 
6417 
7417 
1149 
9839 
4619 
7959 
9634 


75 .85754 84527 5272 
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ent 


490 
666 
898 
1210 
1650 
235 


50 
50 
500 
5C0 
500 
500 
50 
50 


765 
11300 
63000 

230 

318 

4310 
5830 
7810 
10400 
13900 
1860 
259 


50 
50 
500 


500 
500 
500 

50 


69 
1030 
14700 
10 
284 
3850 
5170 
68600 


19000 
15600 
20800 

2870 
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TaBLE II. Roots of W, («/T)—Continued 





n l Tn, Cc €n,3 
24 1 0.15301 84895 9942 0 5 
24 2 0.51360 83707 0681 0 5 
24 3 1.08218 80942 3753 0 50 
24 4 1.86121 64108 9897 0 50 
24 5 2.85418 83431 3952 0 50 
24 6 4.06569 47966 6047 0 50 
24 7 5.50154 23910 2483 0 50 
24 8 7.16891 89382 1722 0 5 
24 9 9.07661 44660 1974 0 5 
24 10 11.23531 54871 5978 1 63 
24 11 13.65799 92756 6470 1 933 
24 12 16.36046 86713 7067 1 13400 
24 13 19.36209 01881 5387 1 1 88000 
24 14 22.68683 57294 7187 1 2580 
24 15 26.36479 43150 0290 1 3470 
24 16 30.43444 26913 9403 1 46300 
24 17 34.94620 46462 8518 1 61100 
24 18 39 .96833 92266 7395 1 80000 
24 19 45 .59738 43134 3094 1 1 04000 
24 20 51.97849 53322 1339 1 1 36000 
24 21 59.35068 02360 3490 1 1 77000 
24 22 68.17050 51810 0240 1 23100 
24 23 79.60945 44050 8424 1 3140 
25 1 0.14689 16257 8146 0 5 
25 2 0.49299 48461 2182 0 3 
25 3 1.03859 20353 8078 0 50 
25 4 1.78584 95467 3969 0 50 
25 5 2.73784 45322 0507 0 50 
25 6 3.89860 82901 4203 0 50 
25 7 5.27322 15403 0644 0 50 
25 8 6.86794 56382 7667 0 50 
25 9 8.69039 61931 4098 0 5 
25 10 10.74977 14241 0897 1 58 
25 11 13.05715 36754 8216 1 852 
25 12 15.62591 13581 6414 1 12200 
25 13 18.47224 24144 4275 1 1 71000 
25 14 21.61592 28735 8341 1 2340 
25 15 25.08136 20017 9402 1 3140 
25 16 28.89913 24031 5677 2 21 
25 17 33.10826 73051 8389 2 27 
25 18 37.75986 10435 6086 2 36 
25 19 42.92302 51864 7370 2 46 
25 20 48.69545 32211 6803 2 59 
25 21 55.22399 77787 5012 2 76 
25 22 62.75044 41754 7154 2 98 
25 23 71.73671 59637 8361 2 129 
25 24 83.36839 49267 2194 2 174 
26 1 0.14123 67262 5809 0 5 
26 2 0.47397 45378 8442 0 5 
26 3 0.99838 34056 2151 0 50 
26 4 1.71638 16871 9240 0 50 
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COEFFICIENTS AND ROOTS OF POLYNOMIALS 


TABLE II. 
l Tal 
5 2.63069 31145 8470 
6 3.74487 77262 0273 
7 5.06340 83123 3855 
8 6.59177 56068 7320 
9 8.33662 63598 0514 
10 10.30594 30256 1368 
11 12.50927 80113 1609 
12 14.95806 12826 7794 
13 17.66600 89930 4484 
14 20.64967 47456 6110 
15 23.92920 78044 9273 
16 27.52942 09021 3584 
17 31.48133 78942 1101 
18 35.82451 67628 4751 
19 40.61069 00156 5943 . 
20 45.90978 68582 2976 
21 51.82061 58754 0514 
22 58.49167 48142 7643 
23 66.16744 93598 1048 
24 75.31508 13581 0590 
25 87.13389 48199 8148 
1 0.13600 12570 4444 
2 0.45636 93581 3000 
3 0.96118 10273 0921 
4 1.65214 28929 8375 
5 2.53167 76794 9454 
6 3.60294 55780 8908 
7 4.86990 60826 2376 
8 6.33740 31557 2642 
9 8.01127 62543 8048 
10 9.89850 34830 7608 
11 12.00738 61381 8993 
12 14.34778 78810 5506 
13 16.93144 76204 2762 
14 19.77239 41164 6035 
15 22.88750 43962 7852 
16 26.29727 10029 7123 
17 30.02688 18438 7846 
18 34.10778 47142 1721 
19 38.58003 51906 4629 
20 43.49597 58235 1797 
21 48 .92632 26588 5339 
22 54.97096 30253 0608 
23 61.77999 23077 3237 
24 69.60052 54573 4751 
25 78.90479 29547 6498 
26 90.90552 80995 1289 
1 0.13114 02048 4340 
2 0.44002 68410 2122 
3 0.92665 89979 7975 
4 1.59256 14353 9651 


Roots of W,, (¢/T)—Continued 
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Roots of W, («/T)—Continued 


8214 
8727 
2107 
5926 
1486 
8264 
2683 
0963 
2334 
9327 
5973 
2692 
8103 
4767 
6414 
0563 
2954 
3640 
7904 
6239 
1621 
5939 
9295 


6314 
2803 
2982 
6714 
0362 
4731 
0930 
5907 
3363 
5631 
9518 
9570 
0260 
4301 
4752 
8314 
5279 
4674 
9695 
1023 
6970 
0558 
5849 
6674 
3268 


4401 


TasB_eE II. 
l Tal 
5 2.43989 88598 
6 3.47148 98767 
7 4.69085 73686 
8 6.10229 80663 
9 7.71097 27209 
10 9.52302 14214 
11 11.54571 07983 
12 13.78762 26083 
13 16.25889 71483 
14 18.97155 08191 
15 21.93989 61709 
16 25.18110 69743 
17 28.71599 40224 
18 32.57009 65076 
19 36.77526 28337 
20 41.37202 26709 
21 46.41330 39270 
22 51.97058 20538 
23 58.14479 02456 
24 65.08757 79811 
25 73.04862 35108 
26 82.50512 46756 
27 94.68291 12582 
1 0.12661 47772 
2 0.42481 56894 
3 0.89453 69967 
4 1.53714 63597 
5 2.35458 60074 
6 3.34938 08943 
7 4.52468 00068 
8 5.88431 36574 
9 7.43286 73847 
10 9.17577 57820 
11 11.11944 12276 
12 13.27138 44455 
13 15.64043 65527 
14 18.23698 62111 
15 21.07330 14271 
16 24.16395 46279 
17 27.52639 39971 
18 31.18172 74331 
19 35.15582 50255 
20 39.48091 56524 
21 44.19798 25434 
22 49 .36051 71053 
23 55.04072 95577 
24 61.34056 90221 
25 68.41319 89167 
26 76.51079 93357 6897 
27 86.11542 16992 
28 98 .46569 76629 3931 
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10 
132 
1800 


16 
21 


338 
430 
5380 
6760 
8480 


1330 
171 
224 
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16 COEFFICIENTS AND ROOTS OF POLYNOMIALS 


TaBLE II. Roots of W, («/T)—Continued 


" l Tn, c €n,t 
30 1 0.12239 13636 8232 0 5 
30 2 0.41062 22202 4781 0 5 
30 3 0.86457 25660 9784 0 5 
30 4 1.48547 36509 1174 0 5 
30 5 2.27507 48206 5174 0 5 
30 6 3.23564 71256 6498 0 5 
30 7 4.37001 73777 7174 0 5 
30 8 5.68161 56855 1702 0 5 
30 9 7.17453 65183 2250 0 5 
30 10 8.85361 57459 8089 0 5 
30 11 10.72452 73673 4938 1 6 
30 12 12.79390 50612 0659 1 8 
30 13 15.06949 56117 3280 1 11 
30 14 17.56035 39915 4328 i 1540 
30 15 20.27709 38755 2512 2 10 
30 16 23.23221 33320 6006 2 14 
30 17 26.44052 46096 4200 2 18 
30 18 29.91973 14194 9155 2 224 
30 19 33.69122 07222 7026 2 284 
30 20 37.78117 59363 7412 2 3580 
30 21 42.22218 88741 5065 2 4440 
30 22 47 .05567 82251 9296 2 55200 
30 23 52.33567 93797 4686 2 68400 
30 24 58.13511 50396 8783 2 84800 
30 25 64.55692 28253 1261 2 1 04000 
30 26 71.75574 28866 9527 2 12900 
30 27 79.98619 96989 1509 2 16900 
30 28 89.73509 12115 8820 2 2020 
30 29 102.25357 28564 7799 2 260 


Legend: 
f.,i: the Ith root of the m-th order polynomial, W,. The estimated error bound is 
| tat — root | S én,: X 10-4 


The significance of the numbers in column C is: 
0: the root was computed explicitly 
1: the root was computed as fa,: = 10/ra. 
2: the root was computed as ra,2 = 20/7n,1 


The formulae of equation (4) can be verified by induction based on the actual 
formal construction of successive derivatives of exp (—«/T) [1]. After these rela- 
tions have been established they can be used with a simple induction argument to 
establish the correctness of explicit formulae which are convenient for other pur- 
poses, albeit not for calculations 


(—1)*"7(l + 1)? --- (wm — 1) 
be mat: (n—D! ; 

The roots were located by the method of false position [2], which depends upon 
evaluation of the polynomials for successive approximations to the roots. The evalu- 
ations were computed by synthetic division using double precision arithmetic on 
the ORDVAC computer at the Aberdeen Proving Ground. The polynomials were 
calculated in one of three forms as described in the legend to Table II. This was 
necessary since the original form of the polynomial was relatively insensitive to 
changes in approximations for certain roots; that is, the remainder was too small 
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TaBLeE III. Symmetric Function Check on Roots 


n bn a1 d, 


2 
6 
12 
20 
30 
42 
56 
72 
10 90 
11 110 
12 132 
13 156 
14 182 
15 210 
16 240 
17 272 
18 306 
19 342 
20 380 
21 420 
22 462 
23 506 
24 552 
25 600 
26 650 
27 702 
28 756 
29 
30 


mt 
»e OO 


dS 
ner A © 


L++t++i+i+i il 
3 


Legend: bani — 2 tas = dn X 107% 


i=1 


compared with the largest term in the synthetic division schema. When this oc- 
curred, one of the other forms of the polynomial was used in order to define the root 
more accurately. 

The error estimates of Table II were based upon locating the root between two 
approximate values which gave remainders of opposite sign. Although this method 
of locating roots is inherently self-checking, a further numerical test was made by 
computing the symmetric function summations of Table III. 
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18 BESSEL FUNCTIONS FOR LARGE ARGUMENTS 


Bessel Functions for Large Arguments 
By M. Goldstein and R. M. Thaler 


Calculations of Bessel Functions of real order and argument for large values of 
the argument can be greatly facilitated by the use of the so called phase-amplitude 
method [1]. In this method two auxiliary functions, the amplitude and phase func- 
tions, are defined in terms of the regular and irregular solution of Bessel’s equation. 
These auxiliary functions have the great advantage that they are monotonic func- 
tions of the argument; moreover, for large arguments these functions are slowly 
varying and, hence, easily amenable to computation and interpolation. 

The Bessel Functions of the third kind, so-called Hankel Functions, are defined 
as follows: 


(1.1) H,® (x) = J,(x) + +Y,(x) 
and 
(1.2) H,® (x) = J,(x) — iY,(x). 


These definitions suggest the alternative definition of the Bessel Function in terms 
of an amplitude and phase function, viz., 


(2.1) A, = | H,®(x) | = | H,®(=) | 
or 

(2.2) H,® = A,et** 

(2.3) H,® = A,ye-** 

so that 

(2.4) J,(x) = A, cos ®, 

and 

(2.5) Y,(x) = A, sin %, 


The phase function #,(x) obeys the first order differential equation 


d%, 2 
(3) = = Ae) 
Equation (3) represents the usual Wronskian condition. The amplitude function 
A, is seen to obey the following nonlinear second order differential equation 


(4) PE) 4 A) + (et — aie) — (7) AP = 0. 





Equation (4) follows trivially from Bessel’s differential equation and equation (3). 





Received 4 June 1957. This work was performed under the auspices of the U. S. Atomic En 


© Aer e 


— 





in a Ear 


NER. DEALS 


(5) 


whe 


(6) 


wh 


(11 


(1! 


(11 


(11 








" BESSEL FUNCTIONS FOR LARGE ARGUMENTS 19 
Defining B,(#) as 
(S) But) = J A,(0, 
where ¢ = 1/x, gives the following differential equations for #,(#) and B,(#) 


(6) or = -[B()}* 


(7) d a) 





+4 eo — 4)#]B,(t) — (B07 = 0. 


A boundary condition implied by equation (2.1) and the asymptotic behavior 
of H,(x) asx — © is 


(8) lim d*B,(t) 








ge ren 
so that near ¢ = 0, B,(¢) is approximately given by the solution of equation (7) with 
d*B,(t) _ 
- -* 
This yields 
(9) B,(t) = (i — (»? — fT. 


Equation (9) represents the so-called JWKB approximation [2] to the amplitude 
function, and 


(10) a) =14 ["(5 - feu ( + 4) 5. 


The JWKB approximation for ,(¢) consists of the substitution of equation (9) 
in equation (10) 

Power series expansion about ¢ = 0 for equations (6, 7) are readily obtained. 
One finds that 





(11) {B,(t)/t} = 1 + ast? + ag + agl® + ag + aie + --- 
where 
(11.2) a= ie 

oS 
(114) a = Zat— Fa 

8 oo 
(116) a= Toe +4 1s) 

_ 195 ,_ 611, , 1821, 315 

(11.8) = 504g * — 256% + 08 rT aby 
(11.10) ayy = 963 gs — 4199 44 4 29811 0, _ 2223 0, , 2835 


— a nu 8 8 
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with 
a — 3); 


a = (y 


and 


(12) ¢ [2 + (+ 5 | = 1 + dof? + dalt + det® + bet? + diol! + --- 


where 

(12.2) - se 

(12.4) — = -¢ 

(12.6) be = oa — JT attia 

(12.8) bs = ay at a + i ni 


Alternatively, the amplitude function may be expressed in terms of the follow- 
ing series 


(13) = = [1 — af® {1 + cot? + catt + cot® + col? + cof + --- }, 
where 
(13.2) cf = 0 
(13.4) — | * 
i 
(13.6) o= 16% +> 
- -* ot 1761 a? — 315 
(13.8) cs = + q28 * 16 Qa 
~ 7 ot 6975 ot 17469 ot as 
(13.10) C10 = ~—_— + Tas * aap + — a. 


Equation (13) represents an expansion about the JWKB approximation. The series 
in Equation (13) is semi-convergent. Substitution of Equation (13) into Equation 
(6) may be used to obtain an explicit “JWKB” series for ,(¢). However, the 
analytic form of the coefficients is so complicated as to be of little practical value. 
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“ The amplitude and phase functions for » = 0, 4, $, 3, 1 are plotted in figures 1 
and 2 respectively. In figure 1 the function 


[) 


is plotted against ¢. In figure 2 the function 


[ {e.) + (» +4) *\] 


is plotted against ¢. In terms of these functions the Bessel Functions J,(x) and 


Y,(x) become 
ie) « [PAO | Pt one &,(t), 


Y,(x) = [742] f sin ®,(2), 
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B,(t 
Ficure 1. The amplitude function, © vs. t for several values of » between zero and one. 


The amplitude and phase functions are defined by the relations 


B, y 
J,(x) = 2 | ~f cos ®, (t) 
Y,(x) = [| mi sin ®,(t), 


t=-, 


x 
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2 
-. 


i= 


For | » | S$ 1, the series equation (11 or 13) and equation (12) yield eight figure 
accuracy in J,(1/t) and Y,(1/t) for ¢ S 4, for the number of terms shown. This is 
in contrast with the usual asymptotic expansions which yield comparable accuracy 
only for ¢ $ 1/20, and which converge much more slowly in the region of good con- 


vergence. 


For » = 4, a = 0, so that the coefficients in equations (11) and (12) vanish. 
Therefore, for » = 4 the series reduce to single terms, yielding 




















(14.1) Iy(x) =. 2 sin x 
TX 
(14.2) Y,(x) = —/2 cos x 
RX 
or 
(14.3) A;(x) = yz 
WX 
(14.4) @(x) =x— ; 
tT ' ! i 
1.10 ‘ 4 
he 
eu 
3 
= 
é 
0.95 od 
l l l l 
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Ficure 2. The phase function ¢ [a + (@ + 4) ‘| vs. t, for several values of » between 


zero and one. 
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It is also clear that the series are most rapidly convergent for v close to 4. 

For computational purposes it is evident that greatest accuracy will result 
through the use of equations (11) and (12) for |» | < 1. The appropriate recur- 
rence techniques may be used to obtain numerical values for | » | > 1. 

Equations (11) and (12) may be regarded as starting series for the differential 
equations equations (6) and (7). For | » | < 1 equations (6) and (7) may be in- 
tegrated numerically starting at t ~ $ and give good accuracy to ¢ ~ 4. Thus 
the method described above can be used to generate values of J,(x) and Y,(x) 
for x => 2. 

The Bessel Functions of imaginary argument J,(x) and K,(x) are amenable to 
similar treatment. The phase-amplitude functions are defined as 


(15) i J C, sinh ¢ 
Tx 


= tf Or 
K, 4 x Cre ° 


These definitions are chosen so that 
(16) C,=1 


and 
3 = Xx, 


I(x) = = sinh x 


or 


BX 
K,(x) oa de g™. 


The phase function ¢, obeys the first order differential equation 


(17) HO — (1D)? 


and the amplitude function D, = tC,, obeys the nonlinear second order differential 
equation 


(is) SO 4 FO + GF - DED. - (DWT = 
The series expansions about ¢ = 0 are 


(19) {Pe = 1 — asf? + ag — agt* + agl* — at + ---, 


and 
ty (t) = 1 — bel? + belt — del® + def? — diol! + ---. 


Alternatively, in the JWKB type approximation, the amplitude function is given by 
the series 


D,{t) 


(21) ; 


= [1 + af®P {1 — col? + cat — col® + cof® — crt + --- j, 








24 BESSEL FUNCTIONS FOR LARGE ARGUMENTS 


where the a,, b, and c, are defined as in equations (11), (12), and (13) and 
1 ‘/1 i 
(22) ot) => + ['(5- fp) au 


The amplitude and phase functions for » = 0, 4, 4, 3, 1 are plotted in figures 
3 and 4 respectively. In figure 3 the function 


[D,(t) /t] 
is plotted against ¢. In figure 4 the function 


[t¢,(t) J 
is plotted against ¢. 
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Ficune 3. The amplitude function 22 





vs. t for several values of » between zero and one. 


The amplitude and phase functions are defined by the relations 


I,(x) = ew R= sinh ¢, (t) 
t ca 
Dt - 
K,(x) = ee hee [-#()) 
t x2 


ti=- 


x 


where 
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Ficure 4. The phase function é¢,(¢) vs. t, for several values of » between zero and one. 


In terms of these functions the Bessel Functions of imaginary argument I,(x) 


and K,(x) become 
Rlsj< EZ R 5 sinh ¢,(2) 


T 
K,(x) = [P| J25 e-Prtn 
t x2 
with 
t= = 
x 


The series, equations (19) and (20), can provide starting values for the nu- 
merical integration of the differential equations, equations (17) and (18). These 
series give seven figure accuracy for | » | < 1 for x = 8, for the number of terms 
shown. However, while straightforward numerical integration of equations (17) 
and (18) yields accurate values for K,(x) for x = 2, the values for J,(x) rapidly 
deteriorate. The reason for this is clear. At x = 8, e.g., the starting values are ac- 
curate to eight figures. The functions I,(x) and K,(x) so calculated may then be 
written as 


I(x) = I,(x) + pK,(x) 
R,(x) = K,(x) + o1,(x), 


where pK,(x) and oI,(x) represent the truncation error, and p and @ are constants. 
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As x decreases, K,(x) increases rapidly while I,(x) decreases rapidly. Therefore, the 
error introduced by the truncation, decreases in the case of K,(x), and increases in 
the case of I,(x). This difficulty is avoided by calculating p. That is, at x = 2, 
e.g., where K,(x) is still accurately given by the solution of the differential equations, 
equations (17) and (18) 

T(x) — I(x) 
(24) —— ar 


At x = 2, I,(x) is easily obtained by standard methods, e.g. the defining series, or 
integral representation. 
M. GOLDSTEIN 


R. M. THALER 
Los Alamos Scientific Laboratory, 
Los Alamos, New Mexico 
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Numerical Experiments in Potential Theory 
Using the Nehari Estimates 


By U. W. Hochstrasser 


1. Introduction. It is well known that the Dirichlet problem can be solved by 
using a set of harmonic functions which are orthogonal over the boundary [7]. 
Prior to the availability of high speed computers, however, the method was not 
frequently used since generating sets of orthogonal functions by hand computa- 
tions is very laborious. With the advent of high speed computers interest in the 
method has been renewed and several papers have appeared recently [1, 2, 3]. 

In using the method, the numerical analyst is particularly interested in esti- 
mates of the error committed by replacing infinite orthogonal expansions by finite 
ones. Nehari [1] recently has given such estimates for a number of Dirichlet prob- 
lems. The purpose of the present paper is to report the results of some numerical 
experiments which were designed tc investigate how close these estimates come to 
the actual errors in certain particular examples. An existing computer program 
[4] was used for the orthogonalization of a given set of vectors on the SEAC com- 
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puter at the National Bureau of Standards. Several boundary value problems, each 
involving the same irregular pentagon, were solved. For two of these problems the 
exact solution is known so that one can compare the exact errors with the esti- 
mates. In all cases the number of orthogonal harmonic functions used to approxi- 
mate the solution was varied. The effect of using a sophisticated choice of the 
harmonic functions was also studied. 

2. Dirichlet problems for an irregular pentagon. In the experiments an irregular 
pentagon inscribable in a circle was used in all the cases as the region on which the 
Dirichlet problem was defined (see Figure 1). 


! 
! 
+ 
| 
! 
! 
! 
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For this region the following Dirichlet problems were considered : 
(a) The boundary values are given by 


U(x, y) = Re (cosz) = Re cos (x + ty) = cos x cosh y. 
(b) The boundary values are given by 


1 ai x+4 
z+4 (x + 4)? + y* 


(c) The boundary values are equal to one on one side and zero on all the others. 

(The solutions of these problems are the harmonic measures for the sides.) 

The solution of the first Dirichlet problem changes sign in the interior and on 
the boundary of the region. In the second case the solution is always positive on the 
whole region. The set of functions which were used to approximate the solution in 
the cases a and 6 consisted of harmonic polynomials obtained by orthonormalizing 
1, Re z, Im gz, Re 2’, Im 2*, - - - on the region. It was expected that with a few func- 
tions one would obtain a better approximation in the second case than in the first. 

In order to find approximate values for the harmonic measures of the sides of the 
pentagon one can start from the harmonic measure of the circumscribed circle. Let 
the segment 2a be that part of the circumference of this circle which is cut out by 
the side of the pentagon with boundary values equal to 1. Then the harmonic meas- 
ure 49 which has the value 1 on the segment 2a and 0 on the rest of the circumference 





U(x, y) = Re 
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[ (x — x2) (x — x1) + (y — 92) (y — 1) 
= arc cos 


Vie — 1)? + (y — nn)* IL (e — m2)? + (y — 92)*] 
XiX%2 + ive 
# arc cos | 7+ 219 + 


where ¢ is the angle belonging to any point in the region and (x, y1), (x2, y2) are 
the coordinates of the end points of that side of the pentagon with boundary values 


Fic. 2 


equal to 1. In addition to these harmonic measures the harmonic polynomials 1, 
Re z, Im 2, Re z*, Im z*, -- - were used for obtaining the set of approximating func- 
tions. These sets were orthonormalized with respect to the inner product 


(S.8) = [Fe rele, y)ds, dst = det + dy* 


where the integration is along the boundary c of the region. 

If m functions 1, Re z, Im z, Re(z?), Im (z?) --- are used as basic set one obtains 
after orthonormalization a system of n orthonormal harmonic polynomials u,(x, y). 
The approximate solution of the Dirichlet problem is given then by the sum 


u(x, y) = >| / U(X, Y)u,(X, ¥ dS |ua(x, y), (dS? = dX? + dY?). 


For the numerical computations the contour integrals appearing in these com- 
putations must be replaced by an integration formula involving finite sums. Since 
the contour consists of straight lines, splitting the integral into five integrals cor- 
responding to the five sides of the contour was suggested. 

For their evaluation an eight point Gaussian-type integration rule was chosen 
involving the end points of the interval of integration as given by Radau [5]. 
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This rule is exact for polynomials up to degree 13. It corresponds therefore to a seven 
point ordinary Gaussian integration rule. Since it involves the endpoints, the com- 
putational work for the eight point Radau rule and the seven point Gaussian rule is, 
however, the same for a closed contour. Radau’s formula has the added convenience 


of using the corner points of the pentagon. Thus the contour integral is replaced by 
a sum over 35 products of the form 


[ Gt, as ~ EL wala +(1—Ajeen Age + (1 — Ayes) 
where (x;, yi) are the coordinates of the vertices of the pentagon 
(xo = Xs, Yo = Ys) 
and the weights w,; as well as the coefficients A; are given by the integration rule 
wiz = 4 V (xi — xen)? + (i — ven)? )|=60fG = 2,3, ---, 7,4 = 1,2, ---,5 
wa = {V(x — xen)? + (9% — ¥en)? 
+ V (xia — x2)? + (ven — yi-2)*} A (x1 = 4, ¥-1 = a) 
j 1 2 3 4 5 6 7 











H; .03571 428 .21070 422 .34112 270 .41245 880 =-H, =H; =H; 
A; 0 -06412 9925 .20414 991 .39535 039 1—A,1—As;1—As 


The orthonormalization is carried out by the Gram-Schmidt process: If f; is a 
set of m vectors with N components then a system of orthonormal vectors u; with 
components u%, is found by the formulas 


“4 = aufi 


Ue = Gnfi + Goofe 


Un = Gnifi + Gaof2 + --- + Gnn—ifa—1 + Gnntn 
where the ai; are determined by the orthogonality and normalization condition 
(us, uz) = 6 = 154 = 7 4,j = 1,2, ---,%. 
=O;t4j 


Since in the continuous case the orthonormalization is with respect to the con- 


tour integral / f(x, y)u(x, y)ds, the inner product is defined by 


N 


(fi, uj) = > Wrf ath in N = 35, wy = Wi,yWe = Wi,2 +++ Wr = Wir, 
k=l 


Ws = W2,1 “Pp Wig = W2,7 
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where the weights w, are given by the integration formula. Here the components 
of the vectors f; are the values of the corresponding harmonic functions of the basic 
set 1, Re z, Im z, Re 2? --- , in the cases a, b; or of the set mo, Rez, --- , in case c, 
at the 35 points given by the integration rule. 

The approximate solution u(x, y) of the Dirichlet problem is given then by the 
sum 


u(x,y) = > difi(x, y) 


dy = D) (U, us)ans #=1,2,-- -,m. 


k=l 


U is a vector the components of which are given by the boundary values of the 
Dirichlet problem in the points defined by the integration rule; f;(x, y) is the har- 
monic function corresponding to the vector fi. 

u(x, y) has the property that among all the possible approximations of the form 


> Cufe (x, y) 


k=l 
the deviation from the exact solution U,(x, y) is least in the least square sense: 
ifA = U.— Dodigfi(x,y) = U.—u 
i=] 
then (A, A) = (U, — u, U. — u) = minimum. 


Nehari [1 ] has given the following estimate for the deviation in the region where 
the Dirichlet problem is defined 


Af as 
2x? Je (X — x)? + (Y — y)?* 


— 5 uelx, 9) |. 





(Uae, 9) — Ladle, 99)? 5 (0,4): | 


3. Numerical results. For each of the three cases a, b, and c the number n of 
harmonic functions in the basic set was varied in cases a and 6 from 3 to 11, in case 
c from 4 to 12. Due to the limitations in internal storage capacity it was not possible 
to go higher than 12. The functions were numbered in the following way 


fi=1, fe=Rez=x, fs=Imsz=y, f, = Rez? =x?—y*, fs = Imz? = 2xy 
fe = Rez’, f; = Im2*, fg = Rex, fy = Ima, fi = Rez’, fi = Ime. 


The variation of was done in order to observe the convergence of the approxima- 
tions to the exact solution of the Dirichlet problem. Due to the limited accuracy of 
the numerical computations (the orthonormalizing code works with floating decimal 
numbers containing only 8 digits) it must be expected that the improvement of 
the approximation with increasing m will level off until the addition of another 
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harmonic function f; to the set will not improve the approximation at all. In the 
three cases this point has not yet been reached, as one can see from table 1, in 
which the norm of the vector consisting of the deviations in the points given by 
the integration rule is listed. 


TABLE 1. Norm of Error Vector A 
Harmonic Functions 


Sis fo, fa fi tofs fi to fr fi tofs fitofu 


18.007 9.565 8.822 2.500 2.061 

351 .214 -157 -107 -076 

-176 .100 -094 081 -057 

-126 077 -026 -023 -022 

.219 -158 149 115 -113 

-099 -059 -038 -037 

cs* -086 .078 


*In these cases %» was used in addition to fi, --- , f,, and c; is the case where 
the boundary values are one on the side 7 (see Figure 1) and zero on the rest of the 
boundary. One observes that in the case a the norm decreases more rapidly be- 
tween n = 7 and 9. The norm for m = 8 is 2.517 and shows that the decrease occurs 
if Re 2 is added to the set. The explanation for this behavior is probably that the 
function Re z* = x‘ + y' — 6x*y? with its changes of sign along the lines x = 
(+ 1 + v2)y (all combinations of signs are possible) fits particularly well the 


prescribed boundary conditions. 
In table 2 the maximum error is given for the points on the boundary defined 
by the integration rule. 


TABLE 2 
Maximum Error in the Boundary Points Defined by the Integration Rule 


Number of Max. bdry 
functions 3 5 7 9 11 value 


5.739 4.825 .480 - 988 12.915 
. 148 109 -087 -066 -595 
.065 — .078 -068 -052 , 
-055 — .029 -031 -028 
115 — .095 -091 095 
-077 .032 -033 -034 
cs* -070 .079 


From this table one sees that with eleven vectors the maximum relative error 
in cases a and 6} is approximately 10% on the boundary. (Since the deviation has 
been evaluated only in a finite number of points it is of course possible that this 
error may be in some points somewhat larger.) In cases c; the relative error is in 
most cases smaller than 10% due to the use of an additional harmonic function u» 
in the basic set. One remarks also that the maximum relative error for the first few 
approximations to the harmonic measures is much smaller than the error in the 
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cases a and b. This is due to the use of the corresponding harmonic measure uo for 
the circumscribed circle which gives a fairly good first approximation. From these 
‘observations one can conclude that from the numerical point of view the basic set 
of harmonic functions which has to be orthonormalized should not be chosen arbi- 
trarily. The present computations show that a sophisticated choice of these func- 
tions creates considerably better approximations for the same number of ortho- 
normal functions. 

The Nehari estimates for the deviations in the point x = y = 0 are given for 
the cases a and 3 in table 3. 


TABLE 3. Deviations at x = y = 0 and Nehari Estimates. 


Harmonic functions fitofs fitofs fi tof; fitofo fitofu 


Case 
a Nehari estimate 5.19 2.74 2.41 .643 .528 
actual error (abs.) -855 -495 .646 .0163 .0133 
6 Nehari estimate .1010 .0612 .0428 .0276 .0195 
actual error (abs.) .0099 .0015 .0050 .0017 .0002 


From this table one sees that the estimate is still considerably larger than the 
actual error and that it converges slower than the actual error with increasing 
number of harmonic functions. If one compares, however, this estimate to the 
maximum error on the boundary given by table 2 it is evident that the Nehari 
estimate gives better error bounds than the one which can be derived from the 
maximum modulus principle for harmonic functions. In table 4 the Nehari esti- 
mates for the cases a and b are given for a number of points in the interior of the 
pentagon and for 11 harmonic polynomials. : 


TABLE 4 
Case (y = 0) x=-2 x=-1 x«=0 «x=1 


Nehari estimate 1.798 .539 .528 .528 .614 1.631 
actual error .052 .094 .0133 .0105 .0181 .406 
Nehari estimate .0663 .0199 .0195 .0195 .0226 .0601 
actual error .0091 .0023 .0002 .0021 .0038 .0086 


Table 4 shows that the Nehari estimates become worse as the boundary of the 
region is approached if they are compared with the estimates given by the maxi- 
mum modulus principle. With respect to the actual error, however, the Nehari 
estimates are further off near the middle of the region than near the boundary. 

In table 5 the Nehari estimates for the harmonic measures of the sides in the 
case where 12 harmonic functions were used are given for x = y = 0 together 
with the maximum error in the points of the boundary given by the integration 
rule. 

In these cases a comparison to the exact error is no longer possible since the 
exact solution is not known. The values show however that the Nehari estimate 
gives a better bound for the error at the origin than the one suggested by the 
maximum error on the boundary. 
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TABLE 5 
Nehari Estimates at x = y = 0 and Maximum Error for Harmonic Measures 


Case Ci C2 Cs; Cs Cs 


Nehari estimate at .0146 .00556 .0239 .00945 .0201 
x=y=0 
Maximum error — .0520 — .0279 — .0945 .0335 .0788 


4. Concluding Remarks. The best approximation in each of the seven cases gave 
values which are exact within at most ten per cent for many points of the region. 
Similar computations by P. Davis and P. Rabinowitz [6] show that with the given 
limited accuracy of the code the addition of further harmonic polynomials beyond 
Im z* would probably no longer improve this result considerably due to the in- 
fluence of the round-off errors. In order to improve the accuracy one probably has 
to use at least double precision in the computations. 

The numerical examples given on the previous pages show that the Nehari 
estimate gives a considerably better bound for the error in the interior of the region 
than the maximum modulus principle. With respect to the actual error the Nehari 
estimate is, however, still off by a factor which varied in the examples from 5 to 
almost 100. It may be that the situation improves if more harmonic polynomials 
are used with a double precision code, but this can be decided only by experiments 
using computing machines which have a larger high speed memory than the SEAC. 
Since the use of orthonormal functions for solving problems in potential theory 


seems particularly attractive if the solution has to be known only in certain points 
of the region it would be interesting to see how much the present results can be 
improved by using double precision codes on the faster large computers like the 
IBM 704 or Remington Rand 1103AF. 

The author wishes to thank Dr. P. Davis for many valuable suggestions and 
G. Blakley, J. P. Menard, and B. Urban for performing the necessary computations. 
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Approximate Relations between 
Series and Integrals 


By P. A. P. Moran 
If we wish to evaluate numerically an integral of the form 
(1) t= [0 f(x)ax 


we may replace it by the sum 


(2) s= ¥ if(nh +2). 
h=_—© 


Here h is the interval width and 6 a position constant. If f(x) tends to zero suffi- 
ciently rapidly at x = + © and is reasonably continuous we may expect (2) to give 
a good approximation to (1). In practice it is often found that (2) is a much closer 
approximation-to (1) than would have been guessed simply by looking at the graph 
of f(x). For example if 


f(x) = (2n)-te#" 


and h = 2 so that the ‘panels’ of the integration formula are very wide compared 
with the scale of variation of f(x), we find that S, which is a periodic function of 6, 
never diverges from J = 1 by more than 0.015 whilst if k = 1, the divergence is 
never greater than 0.52 X 10-* (Aitken [1], Yates [14]). On the other hand a func- 
tion like f(m) = x~'(1 + n*)— gives a much lower accuracy (0.37 X 10-* for h = 1). 

If f(m) tends to zero like e, say, all its derivatives will also tend to zero at 
x = +, and use of the Euler-Maclaurin sum formula would suggest that S = I 
exactly. The fallacy in the argument, explained by Goodwin [5], lies in the fact 
that the Euler-Maclaurin series is an asymptotic one only. Another reason for ex- 
pecting a good approximation arises from the Gregory formula for numerical in- 
tegration (Whittaker and Robinson [13], p. 143) which expresses the integral as a 
sum of ordinates together with correction terms involving the differences at each 
end. 

Another approach to this situation is provided by Poisson’s summation formula 
(Fettis [4] and Salzer [11]). This is 


[F seas = 0440) + ¥ scan} — 23 g(2kei-) 
0 k=1 k=l 


gle) = | $(0) cos mid, 
When f(x) is even it will often happen that g(x) is very small for large x and the 
approximation is then very good. Fettis, [4], gives a number of illustrations of the 
use of this fact. 
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This phenomenon of finding such surprising accuracy is not confined to integrals 
of functions which behave like e~*’. In a problem concerning dilution series of viruses 
(Moran [9 ]) it became necessary to sum the series 


(3) > (p™ ast pm) 
ee) 


where 0 < p < 1. This series is in fact periodic in 6 where = exp —2’. Its mean 
value for varying 6 can easily be found to be 


[tex [—3(2*)] — exp [—4(2*) ]}dx = 2 — log 3 (log 2) = 0.450375, 


and an actual exploration of the field shows that the sum never diverges from this 
value by more than 0.000004. 

Since any integral over a finite range can be turned into one over (— ©, «) this 
suggests a useful method of evaluating finite integrals. As an example consider the 
integral 

_ [exp —u* 
Je Xk +e 


which is tabulated by Goodwin and Staton [6]. For the value X = 0.2 they give 
J = 1.5948. Putting « = exp x* we get 


du 





© 2 _ 3 3 
ee) Ea Ae 


X + exp 2 


Taking intervals of width 0.5 in m we get an estimate 1.58069. A better method is 
to put u = e* so that 





J= e exp (—e* + 4) au 
- X + e& 

and again taking intervals of width 0.5 we get an estimate 1.594812 whose accuracy 
is at least 0.0001. One advantage of this method often is that when we want to 
evaluate an integral for a large number of values of some parameter in it (as X in 
the Goodwin-Staton integral) we can do most of the calculation of the ordinates 
once and for all. A second advantage is that the accuracy increases extremely 
rapidly with decreasing h. It is also clear that the problem of deciding what par- 
ticular transformation into an integral on the range (— ©, ©) is best, is one which 
would repay more study. 

The difficulty in using this technique in practice is to obtain an estimate of the 
error. One method is to calculate S for § = 0, 0.2, --- 0.8, say (or 0.1, 0.3, --- , 0.9) 
and examine the variation by drawing a smooth curve. Since the mean value of S 
is equal to J, J must lie between the estimated largest and smallest value of the 
sum. The trouble with this is that it requires five to ten times as much work. How- 
ever, it may be worth while to do this if the integral is required for many values of 
some parameter since it should be only necessary to estimate the error for two or 
three values of the parameter. 

It may also, in some cases, be possible to obtain a still better estimate by the 
method invented by L. F. Richardson and called by him the ‘deferred approach to 
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the limit’ (for some interesting examples of this see Salzer [12] and Bolton and 
Scoins [2]). We calculate the sum S for a series of values of h (ho, 2ho, --- , ho 
say) and use simple polynomial extrapolation to obtain an estimate for h = 0. Some 
care is necessary in this. The remainder is a function not only of h but also of the 
position of the grid. It could happen that an exact estimate was obtained for a 
particular value of 4 and that for smaller # the estimate moves away again from the 
true value. Moreover, the error may well be oscillatory. If, however, f(m) is an even 
function and 6 is put equal to zero it seems almost certain in most cases that the 
convergence will ultimately be monotonic and this method may be useful. It must 
also be remembered that in many cases the error will decrease faster than any power 
of h. 

In some particular cases an estimate of the error can be calculated explicitly. 
This has been done in a valuable paper by Goodwin [5] who considers integrals of 
the form 


/ “ exp (—x?)f(x)dx 


where f(x) satisfies certain regularity conditions in the complex plane. In this case 
a good estimate for the error is, when f(x) is taken as even, and 6 = 0, 


E(h) = 2exp (—#%*) [f(y — inh) exp (—y*)dy 
and if f(m) does not increase too rapidly as m tends to infinity this is approximately 
E(h) + 2xt exp (—x*h-)f(ixh"). 


This result with f(x) = 1 can be used to evaluate more complicated integrals 
such as the integral of the multivariate normal distribution over the range 


(x1 > 0, +++ 5%, > 0) 


when the correlations have certain prescribed forms (Moran [8]), and also certain 
incomplete integrals of the bivariate and trivariate normal distribution with ar- 
bitrary lower bounds (Das [3]). 

The trouble with these methods when applied to other functions is to evaluate 
the error. If we have an even function then it is clear from (2) that 35/86 is zero 
when 56 = 0 and 6 = $h. Since S is a periodic function it will usually happen that 
one of these is a maximum and the other a minimum. In fact we can expand S in 
a Fourier cosine series and then it will usually be found that the first cosine term is 
much larger than all the rest. Thus we can be practically certain in most cases that 
if f(m) is an even function its integral will lie between the limits 


Y if(mh) and YS f((n + Ph). 
n=O n= 


This can be proved to be true for the integral (Das [3]). 


texp —}? [° exp —4x? 


2x epee 
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and since this is equal to 


(2n)-4 "exp —(de*)dx, 


Das has found two very sharp inequalities for the probability integral. However, if 
we add these sums together we get a vastly more accurate result for the same amount 
of labour. 

It is also of some interest to notice that if 6 is taken as a random variable 
with a uniform distribution on the interval (0, kh) the mean value of S is the in- 
tegral and its variance about this mean, like the error term above, can be found 
as a series in the Fourier transform of f(x) (Moran [7]). 

The above ideas can also be applied to integrals on the range (0, 27) of func- 
tions which are periodic with period 27. In this case the sum over a finite num- 
ber of points is often remarkably close to the integral. If, in addition f(x) is an 
even function it is sometimes possible to obtain very sharp inequalities from above 
and below (Moran 10 ) for such functions as 


Ip = (2)— [exp (x cos 6)dé 


and this technique may be useful in calculating such functions as 


- Io(x) dx. 


I am much indebted to a referee for some useful criticisms which have sub- 
stantially improved this note. 


P. A. P. Moran 
The Australian National University, 


Canberra, A. C. T. 
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A Search for Analogues of the Mathieu Groups 
By E. T. Parker and Paul J. Nikolai 


I. Introduction. This paper describes a search for simple groups of a certain type 
using the UNIVAC Scientific Computer, Model 1103A. 

Except for alternating and symmetric groups, the Mathieu groups of degrees 
11, 12, 23, and 24 are the only known finite, quadruply transitive groups. These 
four groups, together with the Mathieu group of degree 22, are the only known ab- 
stract finite simple groups belonging to no known infinite system. ([1] see “‘Ma- 
thieu groups” in index). The curious properties of these groups provide considerable 
motivation for a large-scale search by digital computer for analogous groups. 

It can be noted that the degrees = 11 and » = 23 are primes of the form 
p = 2q + 1, with g also a prime. The analogues to be considered are groups G with 
the following properties: 

1) G is a transitive permutation group of degree p = 2q + 1, p and g primes; 

2) Gisa simple group, and not merely the group of order p; 

3) G is a proper subgroup of the alternating group of degree ». Such groups 
have been studied by G. A. Miller [2] and first footnote thereof for earlier papers, 
and K. D. Fryer [3]. In these papers it is reported that no G exists for degrees 47, 
59, and 83. Even for these fairly low degrees the computations are of such length 
that details were not published. For a fixed degree p = 2g + 1 the groups to be 
examined are relatively few. For each group tested, it was possible to construct a 
pair of generators in a form admirably suited to digital computation; this form is 
probably not in previous literature. 

The program was run for all degrees p = 2g + 1, p and g primes, 23 < p < 
1823. The values of are 23, 47, 59, 83, 107, 167, 179, 227, 263, 347, 359, 383, 467, 
479, 503, 563, 587, 719, 839, 863, 887, 983, 1019, 1187, 1283, 1307, 1319, 1367, 
1439, 1487, 1523, 1619, and 1823. No G exists for any of these degrees, with the sole 
exception of the Mathieu group of degree 23. Degree 23, for which the facts are 
well known, yielded an invaluable check on the validity of the program. Inclusion 
of the next three degrees provided a check both on the program and on the tedious 
computations of earlier investigators. 

In Section II are the arguments which lead to description of a pair of generators 
of a group to be tested, and to the criteria by which such a group may be rejected 
as the alternating group. In Section III the code for the UNIVAC Scientific Com- 
puter is described, with comments on running the program. In Section IV the 
results are summarized. 

Il. Mathematical arguments. As noted in the Introduction, the groups G under 
consideration have properties: 

1) G is a transitive permutation group of degree p = 2q + 1 > 23, p and q 
primes; 

2) G is simple, and not of order p; 

3) G is a proper subgroup of the alternating group of degree p. 

Received 27 ne wow! The program was written and run during the summer of 4956, 
while the authors were employed by Remington Rand UNIVAC, Division of Sperry Rand Co 


tion, St. Paul, Minnesota. The authors vite express their thanks to Remington Rand UNI AG 
for providing’ time on the UNIVAC Scientific Computer, model 1103A. 
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The alternating group of degree p is disregarded as uninteresting in this search. 
Any group containing odd permutations has a normal subgroup of index 2, hence 
is composite if of order greater than 2. The group G is required to be simple (prop- 
erty 2), and of order > p > 23 (properties 1 and 2). The motivation for demand- 
ing property 3 is established. 

The group G, being transitive of degree » (property 1), is of order divisible by 
p ([4] p. 170, Corollary). Thus G has a Sylow subgroup P belonging to the prime 
p. This subgroup P is of order ~ only, hence cyclic. Were P in the center of its 
normalizer in G, then the commutator subgroup K of G would be of order prime 
to p ([4] p. 327, Theorem II). A transitive group G of prime degree is primitive. 
The subgroup K is a normal subgroup of G, hence either the identity or transitive 
on the # letters (([4] p. 196, Theorem X). If K is the identity, then G is abelian and 
in turn of order ~, which is ruled out by property 2. On the other hand, K transitive 
on the # letters implies that K is of order divisible by ~, conflicting with the as- 
sumption that P is in the center of its normalizer. The only alternative is that G 
has an element b such that b-'Pb = P, but b-'ab + a, where a is some non-identical 
element of P. 

It is expedient for computational purposes to identify the p letters permuted 
by G with the integers reduced modulo ~. Without loss of generality a generator a 
of P may be chosen as x — x + 1 (mod #). Being in the holomorph of P, 6 is of 
the form x — tx (mod p), ¢ a fixed integer + 0 or 1 (mod p). Since p — 1 = 2g, 
q an odd prime, 6 is of order g, 2, or 2g. The element 5 fixes only the zero residue 
class and permutes the remaining 2g letters in cycles of equal length. Hence the 
element 6 of even order is an odd permutation, excluded by property 3. Accordingly 
b is of order g. As the holomorph of P has a unique subgroup of index 2, b may be 
chosen as x — r°x (mod p), r a primitive root of p. Thus chosen, } consists of two 
cycles, one permuting the quadratic residues of p, the other the quadratic non- 
residues. 

The symmetric group of degree 2g + 1, g an odd prime, has its Sylow g-sub- 
groups of order g*, generated by two disjoint cycles of length g. If G had as an 
element a single cycle of length g, then G would be p — g + 1 = g + 2 times transi- 
tive ({4] p. 207, Theorem I); noting that a transitive group of prime degree is 
primitive. In turn G would contain the alternating group [5]. Thus G has its 
Sylow q-subgroups of order g only; let Q denote the Sylow g-subgroup generated 
by b. 

Assume that Q is in the center of its normalizer. Then K is of order prime to gq. 
The normalizer in G of P is of order pq; hence the normalizer in K of P is of order 
only p. The second derived group of G is then of order prime to , thus the identity 
alone. The group G must then be solvable ([4] p. 56), in conflict with property 2. 
Thus G has an element c such that c' Qc = Q, c! bc € b; then c* be = bY, w  O, 
1 (mod g). Two possibilities arise: either c fixes set-wise the two cycles of }, or ¢ 
maps each cycle into the other. The search may with no real loss be restricted to 
minimal extensions of {a, b}, so that c may be assumed to be of prime order. (Should 
a G be found, it would then be of interest to drop the restriction that c be of prime 
order, repeating the computations with c’ normalizing Q, c a power of c’. This was 
not done, as no G was found with c of prime order.) Any element c of prime order 
interchanging the two cycles of 6 must be of order 2, consisting of g transpositions; 
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such an element is not in the alternating group. Thus c may be chosen of prime 
order, fixing set-wise the cycles of quadratic residues and quadratic non-residues 
constituting 5, and transforming the cycles into like powers of themselves. Accord- 
ingly c is of the form, (x* — rx**) .(—x? — nx**) (mod p), where exponent & is 
fixed, x ranges over all residues (mod ~), r and m are a fixed quadratic residue and 
non-residue of ; —1 is a quadratic non-residue of , since q is odd and in turn 
pb = 2q + 1 = 3 (mod 4). The holomorph of the group of order g has a unique 
subgroup of order gh, h any divisor of g — 1; thus {b, c} = {b, c:}, where c; is any 
power of any conjugate of c, provided that c and ¢; are of the same order. Thus it is 
necessary to consider only those c of the form (x? — rx**) -(—x? — —x*) (mod 9). 
This form of the generator c was referred to in the Introduction. Other mathe- 
matically equivalent notations [1, 3] would be very awkward to simulate on a 
computer. A further observation reduces essentially to half the number of groups 
to be examined. Note that y — —y (mod ) transforms a into a“ and {b, c} into 
{b, c*}, where c and c* are respectively (x? — rx%)-(—x? — —x*) and 
(x? — rx) .(—x? — —x*) (mod p). The powers of 4 include all quadratic 
residues of p. Accordingly, for a given p and exponent k, all groups to be examined 
have c of the form (x? — 4‘x**) .(—x? —- —x*) (mod p),0 <i < (p — 3)/4. 
For each value of 7, the permutation c belonging to 7 + 1 is obtainable by multiply- 
ing by 4 (mod p) the maps of all quadratic residues of p in the previous permu- 
tation. 

For each prime divisor \ of g — 1 it is necessary to select only one value of k 
yielding c of order \. Thus the number of groups to be examined for a given degree 
p = 2q + 1 is (p + 1)/4 times the number of prime divisors of g — 1. 

For a given p and k, the c with 7 = 0 is of the simple form y — y* (mod 9), 
provided that & is odd, this condition being easily satisfied. As noted above, p = 3 
(mod 4). Further, g = 1 or 2 (mod 3), and g = 1 (mod 3) implies = 0 (mod 3), 
which is impossible for p a prime > 23. Hence g = 2 (mod 3) so that p = 2 (mod 
3). Thus p = 11 (mod 12). By quadratic reciprocity, —3 is a quadratic non-residue 
of p. The multiplicative group of residues prime to p being cyclic of order 2g, —3 is 
a primitive root of p Thus c with i = 0 is (0 — 0) -[(p — 3)’ -— ((p — 3)*)*] 
(mod p), 1 < j < p — 1. This observation made unnecessary the direct genera- 
tion in the computer of large blocks of k-th powers. 

Any of the groups {a, b, c} is generated by a and c alone. For {a, c} is transitive 
on the # letters, the normalizer in G of P is of order pg, and P is not in the center 
of its normalizer in {a, c}. 

Each group {a, c} is simple ([3] Theorem 4.1), satisfying property 2. Thus 
any {a, c} of the type considered is either the alternating group or a G having the 
three properties. 

Since = 2 (mod 3), p — 2 is composite. Thus any group {a, c}, with an ele- 
ment consisting of a single cycle of length a prime x < 9, is the alternating group 
of degree p. The argument is given above for the special case = g. 

The program was designed to construct the permutation c for each group, then 
to decide whether any element a‘c (1 < i < p — 1) has any power consisting of 
a single cycle of odd prime length  < p. (x = 2 is ignored, since the alternating 
group has no element with only one cycle of even length.) Either a group was re- 
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jected as alternating; or, when testing all a‘c failed to show the group alternating, 
an output of c was provided. 

Ill. Description of the Code for the UNIVAC Scientific Computer. For each 
prime value p = 2g + 1, 23 < p < 1823, the value g — 1 was factored. Then for 
each prime divisor \ of g — 1, an odd integer k was found such that = 1 (mod 
q), k = 1 (mod q). These 88 9, & pairs together with about 300 program instructions 
comprised the input. 

The routines, with numerical designations, are described below for a given 
b, k pair. 


1. Generate and store the quadratic residues of p. Generate and store in ad- 
dresses y, 1 < y < p — 1, (high-speed magnetic core), the numbers y* 
(mod p) ; store p in address p. Addresses 1 to p are thus loaded with c (i = 0) 
for the first group of the block of (p + 1) /4 groups for the pair p, k. A copy 
of c is also stored on the drum for possible output and for proceeding to the 
next group (Routine 5). Routine 2 is now entered. 

2. Multiply the permutation in addresses 1 to p by the permutation 


a:x— x + 1 (mod p). This routine is rapidly accomplished by using the Re- 
peat instruction of UNIVAC Scientific to transmit the contents of address y to 
address y — 1,2 < y < p. The initial contents of address 1 are then placed at 
address p. Addresses 1 to ~, having been loaded with a‘c, now contain a**'c, 
0 <i < p — 2. Proceed to Routine 3. 


3. Count and store the lengths of cycles of the permutation stored in addresses 1 to p. 
The principal idea of this routine is the following: when an address has been 
accounted for in a cycle, the contents of this address are marked by a minus 
sign. Each cycle is traced through its course by proceeding from an address 
to the address indicated by the contents of the preceding. A sign test indi- 
cates the completion of a cycle. The next cycle is then begun at the first 
address whose constants are positive. The process necessarily terminates with 
the contents of addresses 1 to # all negative. Routine 4 is entered next. 

. Test the set of cycle lengths to decide whether the group must be alternating. 
This routine makes use of the criterion developed near the end of Section 2. 
Reformulated for computation, the proposition is: If a‘c has a cycle of prime 
length x < , and no other cycle of length divisible by x, then {a, c} is the 
alternating group. « = 2 can be ignored, since the alternating group con- 
tains no element with only one cycle of even length. 


A technique previously used by E. Lehmer ([6] p. 105) was employed to de- 
cide whether a given odd number is prime. Consecutive bits of the contents of 
consecutive addresses were loaded (as part of the program data) with one or zero 
according as consecutive odd numbers 1, 3, --- , 1823 are prime or composite. 
Having an odd number Z to be tested, 1/2 (Z — 1) is divided by 36, the number 
of bits per word. The quotient then indicates the address and the remainder the 
bit in that address containing the information as to whether Z is prime. 

Each cycle length L; is checked as it is found by Routine 3. If L; is a prime 
such that (p — 1)/2 < L; < >, enter alternative B described below. In this case 
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all other cycle lengths are <>, hence indivisible by p. If L; = p, enter alternative 
A; the group cannot be rejected on this permutation. 

If Routine 3 is completed without the group having been rejected, the cycle 
lengths are summed. If this sum were not ~, the machine would stop. This error 
stop never occurred. Now for any prime cycle length Z among the L;, the product 
1; L; # 0 (mod Z*) indicates rejection of the group, and alternative B is entered. 
If the group is not rejected on any prime cycle length Z, enter Alternative A. 

Alternative A: Return to Routine 2 except when all a‘c, 1 < i < p — 1, have 
been checked by p — 1 passes through routines 2, 3, and 4. If all a‘c have been 
tried, print a typewriter character indicating output, and punch a tape of the 
permutation c. Proceed to Routine 5. 

Alternative B: Print a typewriter character indicating rejection of the group 
as alternating. Proceed to Routine 5. 


5. Generate c for the next group in the block of (pb + 1) /4 groups for this p, k pair. 
Multiply the contents of each address y of c, 1 < y < p which is a quadratic 
residue of » by 4 (mod ). Write this new c over the old one; enter Routine 
2 as before. 


The program was run for 15,924 groups, requiring somewhat less than 6} hours 
on the UNIVAC Scientific computer; that is, an average of about 1.5 seconds was 
consumed per group. For the groups of high degree, the typewriter indication of 
rejection of a group as alternating was nearly rhythmic. Some calculation of in- 
struction times made it clear that most groups of high degree were rejected on ac, 
the first permutation tested. 


IV. Results. It was demonstrated by the program that for p = 2q + 1, p and 
q primes, 23 < p < 1823, the only G of the type sought is the Mathieu group of 
degree 23. An output tape was prepared for exactly one other group, namely p = 
47, with c the initial permutation x — x** (mod 47). This occurred through no 
program or machine error. In this isolated case all a‘c fail to satisfy the criterion 
for rejection. Longhand multiplication of a few permutations showed the group 
to be alternating. Actually it is a bit surprising that all but one group were rejected 
by this criterion. 

The authors somewhat expected to construct at least one G with this program. 
It seemed plausible that the G might be quadruply transitive. Instead, the Mathieu 
groups remain isolated more than before. The results of the program yield good 
evidence for the conjecture: An unsolvable transitive group of degree p = 2q + 1, 
p and q primes, is alternating or symmetric, except for p = 7, 11, and 23. (There is 
exactly one exception for each of 7 and 23, and two for 11.) 

ADDED IN PROOF. Paul J. Nikolai, currently with Wright Air Development 
Center, Wright-Patterson AFB, Ohio, has obtained computer time to continue 
the problem. The machine used is an 1103A with 12,288 words of core storage. 
The large storage capacity rendered unnecessary several time-consuming instruc- 
tions of the earlier code. In addition, the authors improved the logical organization 
of the program. The result was a speed better than twice that reported above. 
As for degrees 23< $< 1823, all groups generated were rejected as alternating. 
The further degrees run are p=1907, 2027, 2039, 2063, 2099, 2207, 2447, 2459, 
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2579, 2819, 2879, 2903, 2963, 2999, 3023, 3119, 3167, 3203, 3467, 3623, 3799, 3803, 
3863, 3947, 4007, and 4079. 


The Ohio State University, 
Columbus, Ohio 


1. Ropert D. CARMICHAEL, Introduction to the Theory of Groups of Finite Order, Ginn and Co., 
New York, 1937. 
2. G. A. MILLER, ‘Transitive wot 4 degree p = 2g + 1, p and gq being prime numbers,” 
one Jn. of Math., v. 39, 1908, p. 210-21 
085 5 KD Ferme, “A class of Rovotins od groups of prime degree,” Canadian Jn. Math., v. 7, 
1 
’ Ww. BuRnsIDE, Theory of Groups of Finite Order, 2nd ed., Cambridge Univ. Press, Cambridge, 
tot Dover Publications, nc., New York, 1955. 
. G. A. MILier, “Limits ‘of the degree of transitivity of substitution groups,” Am. Math. 
on * Bull., v. 22, 1915, p. 68-71. 
6. EMMA LEHMER, “Number theory on the SWAC,” Proceedings of Symposia in Applied 
ates, v. 6, McGraw-Hill Book Co., Inc., New York, for the Am. Math. Soc., 1956, p. 
103-1 


On the Accuracy of Implicit Difference 
Approximations to the Equation of Heat Flow 


By Wolfgang Wasow 


1. Introduction. Implicit finite difference methods are gaining favor in the 
numerical solution of initial value problems for partial differential equations, because 
their computational stability is comparatively insensitive to the mesh ratios. The 
difference equation problems of importance in the applications are mostly too com- 
plicated to permit a theoretical analysis of their approximating qualities.It is a 
reasonable working hypothesis that the properties of the simplest linear difference 
equations with constant coefficients are indicative of the situation to be expected 
in many more involved problems. 

The present paper is devoted to a study of the convergence, the stability and 
the truncation error of implicit difference approximations to the initial value prob- 
lem defined by the differential equation 


(1.1) L[u] = u: — uz = 0, 0<x<z, 0<t<T, 
where the subscripts indicate partial differentiations, and the subsidiary conditions 
(1.2) u(x,0=f(x), u(0,t) = u(x,t) = 


in which f(x) is a given function. The approximating difference equation to be con- 
sidered is 


(1.3) LanlU] = kO{U(x,t + k) — U(x, t)} — sh-*{ U(x + h,t + k) 
— 2U(x,t +k) + U(x —h,t+k)} — (1 — s)k-*{ U(x 4+ &, t) 
— 2U(x, t) + U(x — h,t)} = 0, 0<x<z, @.< 4 SF. 


eceived 1 May 1957. Sponsored by the United States Army under Contract No. DA-11- 
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The function U(x, ¢) must also satisfy the same subsidiary conditions as u(x, ¢), 
i.e., 


(1.4) U(x,0) = f(x), U(0,#) = U(x,t) = 0. 


It will be assumed throughout that +/h is an integer. Furthermore, h and k will be 
restricted to intervals0 <hkh <i <1,0<k< ki <1. 

The cases s = 0, s = 1 have been extensively studied in the literature (See, 
e.g., [1], [2], [3]). A paper in process of publication by M. L. Juncosa and D. 
Young [4] contains an appraisal of the truncation error for s = }. This appraisal 
is less sharp than the one contained in theorem 2 of the present paper. P. Lax and 
R. Richtmyer [5] have applied their general theory of finite difference approxima- 
tions for initial value problems to an equation slightly more general than (1.3). 


. 4 
Their arguments are based on the norm { i ¢*(x, t)dx> . The elegant methods by 
0 


means of which they study the convergence and stability can doubtless be extended 
to yield appraisals of the truncation error as well, but since uniform or pointwise 
appraisals are more satisfactory to the numerical analyst than bounds on the mean 
square norm, it is hoped that the results of the present paper are of some interest. 

2. The Convergence of U(x, t). If f(x) is integrable and }°@_: }, sin rx is its 
Fourier sine-series the function 


(2.1) u(x,t) = > exp (—r*t)b, sin rx 
r=] 


is, for ¢ > 0, a solution of the differential equation (1.1), as can be readily verified. 
Under rather mild additional conditions on f(x) this function fulfills also the 
boundary conditions (1.2), at least in some weak sense [3]. Here, the stronger 
hypothesis 


(2.2) > | a | < © 
r=l 


will be made, which is, for instance, satisfied if 
(2.3) f(0) = f(x) = 0, 
and f(x) is continuous and has a piecewise continuous derivative. 


The method of separation of variables leads to the formal trigonometric series 


(2.4) U(x, t) ~ >> E,*!*b, sin rx 


r= 


for the solution of (1.3) and (1.4). Here E, is defined by 


h? — 4(1 — s)k sin? (rh/2) 
i? + 4sk sin® (rh/2) 


The derivation of (2.4) is straightforward and will be omitted. 
If 


(2.5) E, = 





(2.6) 2(1 — 2s)k — h? < 0, 
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and only then, it follows that 
(2.7) | Z| <1, 


for all r. In this case the series in (2.4) is uniformly convergent, thanks to condition 
(2.2), and represents therefore the solution of the difference problem. 
In order to prove the convergence of U(x, t) to u(x, t) it will now be shown that 
(2.8) lim £,*/* = exp (—r*). 
hk 
If k is so small that | kr#(1 — s) | < 1 and | kr*s | < 1, it follows from part a) of 
the Lemma in section 5 that 


| log E, + 4h-*k sin? (rh/2) | < c(4h-*k sin® (rh/2))?, ¢ a constant. 


The right member is a bounded function of h, and it can therefore be concluded 
from this inequality that 


lim (t/k) (log E, + 4h-*k sin® (rh/2)) = 0, 
e0 


uniformly in h. Hence, 
(2.9) lim E,*/* = exp { —4h~*t sin* (rh/2)}, 
70 


and (2.8) is now obtained by letting h tend to zero. 

Since, under the hypotheses made, the series (2.4) converges uniformly, for all 
small 4 and k, the passage to the limit indicated by (2.8) can be performed termwise 
and leads to the series (2.1). Thus the following theorem has been proved. 

THEOREM 1. If the Fourier sine coefficients b, of f(x) satisfy the condition 


DL lo| <@, 
r=1 


and if h and k tend to zero in such a way that 2(1 — 2s)k — h*® < 0 then 


lim U(x, t) = u(x, t), 
hk 
where U and u solve the problems (1.3), (1.4), and (1.1), (1.2), respectively. 

For s < 3, the inequality (2.6) restricts the size of the ratio k/h®. For s > 4 
theorem 1 proves the convergence of U to u for any form of the passage to the limit 
h— 0, k— 0. This goes beyond the frequently made statement that k/h? may have 
any constant value in this case. Observe that s is not restricted to the interval 
O<s<il. 

Coro.iary. Denote by u(x, t) the solution of the difference-differential equation 
problem obtained from (1.3) by substituting U, for k{ U(x, t + k) — U(x, t)} and 
by imposing the subsidiary conditions (1.4). Then 

lim U(x, t) = u(x, t), 
k0 
provided (2.2) is true. 

To prove this corollary it suffices to refer to (2.9) and to show that the uniformly 

convergent series obtained from (2.4) by substituting exp { —4h~*t sin* (rh/2) } for 
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E, satisfies the difference-differential equation and the boundary conditions (1.4). 
This can be done by simple verification. The condition (2.6) is, of course, true for 
all sufficiently small k, if h is kept fixed. 

It is easily seen that an analogous corollary could be stated concerning the dif- 
ference-differential equation obtained if the difference quotients with respect to x 
in (1.3) are replaced by the corresponding derivatives. Here, however, the condition 
s > 3 has to be imposed in order to maintain the inquality (2.6). The solution of 
this problem is the series obtained from (2.4) by substituting 


1 — (1 — s)kr 
1 + skr’ 





for E,. 

Such difference-differential equations are sometimes used on analog machines. 
The arguments of this section justify these procedures mathematically. 

3. The Truncation Error. In this section the order of magnitude of the truncation 
error U — u will be determined (see theorem 2 below). The condition (2.6) is known 
to be necessary at least for s = 0, if convergence for all reasonably regular initial 
functions f(x) is desired. Very likely it is necessary for all 0 < s < 4. Various heuris- 
tic arguments lead to the conjecture that s = } is the most favorable choice as far 
as the order of magnitude of the truncation error is concerned. For instance, if 
Ly,x(u] in (1.3) is developed in powers of 4 and k, where wu is a solution of (1.1), 
the resulting series begins with terms of order 0(% + h*), except if s = 4, in which 
case the leading terms are 0(k? + h*). The results of this section will confirm the 
expected exceptional role of the value s = 4. Since the purpose of implicit methods is 
to gain more freedom in the choice of & and k, the range 0 < s < 3, in which the 
condition (2.6) is necessary, is of little interest. Throughout this section the in- 
equalities 


(3.1) k<s<i 


will therefore be imposed. 

In view of theorem 1 the convergence of U to u takes place no matter how h and 
k tend to zero, if (3.1) is true. However, the order of magnitude of U — u is very 
likely affected if h is esponentially small with respect to k. Here the inequality 


(3.2) h/k> p> 0, p a constant 


will be assumed to be true. Without too much complication the arguments apply 
also to the more general case that h*/k > p > 0,8 > 0. The restriction to (3.2) is 
motivated by the desire to avoid involved formulas. The condition (3.2) is consid- 
erably weaker than the hypothesis k/h? = const frequently met in this connection. 
It is also milder than the assumption k = 0(h/log h) made in [4]. 

Finally, the conditions on f(x) must be tightened in this section. It will be as- 
sumed that 


(3.3) |b, | < Kr-*. 
This is, for instance, the case if (2.3) is true and f(x) possesses an integrable third 
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derivative in 0 < x < 7, for three integrations by parts show that then 


(34) 6, = 2 [* f(x) sin rade = 2 C(-1)97"(e) -— 7”"0)] 
TO T 


_ 2 rf" 7") cos rxdx. 
Tr 0 


The following theorem will now be proved. 

THEOREM 2. Let} Ss<i,h/k>p>O0k<et<T,h<h < p/2. Assume 
that f(x) possesses an absolutely integrable bounded third derivative and that f(0) = 
f(x) = 0. Then 

| U(x, t) — u(x, t) | Sct (h*? + k*), 
where a = 1 for s # 4, and a = 2 fors = 4. Also, 


c=M sup (|f’(x) | +|f’"(x) |). 
OsS2S3f 


The constant M depends on s, p and T only. 
The truncation error will be appraised by splitting the series in 


(3.5) U(x, t) — u(x,t) = >> (E,*/* — exp (—r%))b, sin rx 


r=] 


into four parts Vi, V2, Vs, V4 corresponding, respectively, to the ranges r < 1, 
11 <7 < f2,%2 <1 < 1s, 73 <r. The integers r;,7 = 1, --- , 4, are defined below. 
a) The range r < 1;. The numerator 
h? — 4(1 — s)k sin? (rh/2) 
of E, vanishes if 


: h 
7 = 2/h arc maw EVAC EE) 9 


Let ro be the smallest positive determination of the right member of (3.6)if the 
latter is real, i.e., if h/(2~/ (1 — s)k) < 1. Otherwise, set 7 = ©. If 7 is finite it 
lies in the interval 


(3.6) 


cog SE bolls Pare Ce 
(3.7) Va — sk Sm 4/2 Va— sk 


because 1 < (arc sin g)/q < 2/2, for 0 < g < 1, if the principal value of arc sin g 
is taken. Define 7; by 


(3.8) 1 =([Vep/h], O<06<1. 


Here [z ] designates the largest integer not exceeding z. 

The letters c;, 7 = 1, 2, --- , below denote positive constants that may depend 
on s, p. and T, but not on &, k, x, t, or f(x). 

In order to appraise V; let r < 11, and observe that, by (3.2), (3.7) and (3.8), 


(3.9) ty S 1%/N2. 
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Also, if z, is defined by 

(3.10) 2, = 4kh~ sin* (rh/2), 
it follows from (3.8) that 

(3.11) 2 < kr;? < 86, 
Since 


1— (1 — s)s, 
1 + sz, ; 


an application of the Lemma of section 5 shows that 

(3.12) E,‘* = exp [—t(2, — R)/k], 

where R satisfies the inequalities 

(3.13) —t1%,? £ RS }s%,*, 

If S is defined by 

(3.14) S=r — 2k? + Re, 

formula (3.12) can be written 

(3.15) E,*/* = exp (—r*t + St), 

so that 

(3.16) E,*/* — exp (—r*t) = exp (—r*t) (exp St — 1) = St exp ((S* — r*)Z), 


£E, = 


with S* having a value between 0 and S. The quantity S can be appraised by in- 
serting for z, its expression from (3.10) and by applying the inequalities (3.13) and 


0 < 7? — sin’r < art. 


This leads to 


(3.17) ~ckrt'<S< a rh? + 4stkr! 


and therefore, if (3.2), (3.8) and (3.11) are recalled, to 
(3.18) St<S< anh + bstkrt? < ne + 4st, 


Applying the inequalities (3.17) and (3.18) to (3.16) the appraisal 
(3.19) | Z,'/* — exp (—r*t) | < tcs(h? + k)r* exp (—3r%t/4), rin 
is obtained. 

The expected improvement in the order of magnitude, when s = $ can be based 
on part c) of the Lemma in section 5, according to which 


—c2,' < R <0, 
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so that (3.17) can be replaced by 
1 1 
_ 276 _— 4h? — #2 
ak oS ry li < 74 
Inserted into (3.16) this leads to 
(3.20) | £,*/* — exp (—r*t) | < tes(r*h? + r°k*) exp (—r#t), s=4, rion 


The quantity V; is now readily appraised by means of (3.3), (3.19), and (3.20). 
One has 


(3.21) | Vil = > (E,*/* — exp (—r*t))b, sin rx 
r=1 


< tes(h? + k)K Dor exp —3r%/4 
r=l 


< tes(h? + k)K [ r exp (—3r%t/4)dr = caK(h® + k), 


and for s = 3, in addition, in view of (3.20), 


(3.22) | Vil < tak { [ Vipemp (ar [ ye (—rar} 


< tcgK (h? + &*), s =}. 


b) The ranger; < r < ro. The quantity EZ, is a decreasing function of r and posi- 
tive for r < 1. Hence, E,‘/* is a positive decreasing function of r in r < ro. Let re 
be defined by 


(3.23) re = [ro] 
(so that re = ©, if 7s = «), then 
O< E* < EH, nersn 
or, in view of (3.19), 
| E,** — exp (=r) | < EX + exp (—1:%) < | Eff — exp (—1:%) | + exp (—1:%) 
< {tcs(h? + k)ri* + 1} exp (—11°t/2). 
By virtue of formulas (3.3) and (3.5) this implies 


| Val < {tea(h® + Brit + 1} exp (—ri4/2)K Sort 


The last summation is less than [ r~*dr = 2r,*. Now we assume that h; < 0p/16, 
rT. 
so that r;? > 0p/h — 2~/ @p/h + 1 > 0p/2h. Then the inequality for | V2 | becomes 
(3.24) | Ve| < {cr(h® + k) + (8h/(0p)) exp (—t0p/4h)}K 
< cet IK (h? + k), for 71; < 7 < fe. 
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c) The range rz <r < 13. The set r > 12 is empty, if 7 = ©. In the opposite 
case, E, is a decreasing negative function in 7 <r < x/h. Let 
%3 = a/h 
‘Then, for re < 7 < 7s, 


@ _ (4(1 — s)k sin? (rh/2) — <3 
tsk | = tik =— 
| E, | | E, | ( 4sk sin? (rh/2} + 2 


< (Ae —s)- “" < (Hs _ ) 
Pg rks +1 ~\rks + 1 : 


Here we have made use of the monotonic dependence of E, on the variable 











4k sin 2rh/2. 


As r*ks > roks > s/(1 — 5), part b) of the Lemma in section 5 can be applied to 
the last bound for | Z,‘/* |, provided s > 4. This leads to 





1 — 1/r*ks 


/ tied 
(3.25) |Et*| < (; + 1/Pks 


tlk 
) < exp (—2t/r?k?s), te <7 Srp. 


By continuity, formula (3.25) remains true for s = 4. Hence, using (3.3) and (3.25) 








T3 T3 #1 
(3.26) > 6,E,*/* sin rx | < K >> r-* exp (—2t/r*k?). i 
r=. TT. iq 
The function of 7 in the last summation assumes its maximum for r = 7, = +/ 4t/3k?. 5 


It is an increasing function of r for r < r», a decreasing one for r > r,, and the value 
of the maximum is less than ¢~*/*k*. Therefore, (3.26) can be replaced by 








(3.27) : i b,E,*/* sin rx | < K [ r~* exp (—2t/r*k*)dr + Kt-*/2R* 
r=r2 0 
< K(k?/4t + t-*/*k*) < cyt Kk’, for t > k. 
On the other hand, 
(3.28) > b, exp —r*t sin rx | < K [ exp (—r%t)r-*dr 
r=r2 rs 








= — (1/21) exp (—rtr-*) 





— (2/t) / exp (—r%t)r-“dr < t°Kr,-* < ¢°Kk?, 
rT? rT: 


because rs > % > ((1 — s)k)~*. Combination of (3.27) and (3.28) proves that 
(3.29) | V3 | < ¢ut"Kk?, t>k. 





d) The range r > 13. For r > rs = x/h, 






(3.30) Inl< HS lel<k [ dr < caKht. 


r= jh 


hh 
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The desired bound for the truncation error is now obtained by adding the 
inequalities (3.24), (3.29), and (3.30) to (3.21) or (3.22), respectively. This com- 
pletes the proof of theorem 2. 

The order of magnitude of this appraisal cannot, in general, be improved. One 
consequence of this theorem is that the particular choice s = 4 improves the order 
of magnitude of the truncation error only if h = 0(+/k). 


As was mentioned before, it is possible to carry out the arguments of this section 
with the milder hypothesis 


W/k>p>O0, B>0, 


taking the place of (3.2). However, higher powers of t-' will then appear in the bound 
for the truncation error. If k shrinks so much more slowly than h that 


| log |"/k < p < @, 


then the method of proof breaks down decisively. 

4. Stability. The term computational stability refers to the extent to which 
small errors in the data and in the stepwise computation of U(x, t) affect the final 
result. In the case of linear homogeneous difference equations the cumulative mod- 
ification of the solution by such errors is the sum of the departures caused by errors 
on each line ¢ = 0, k, 2k, --- . Thus, a knowledge of a bound on U(x, ¢) in terms of 
the initial function f(x) conveys information concerning the stability of the pro- 
cedure. In the previous section it was decisive that f(x) had certain smoothness 
properties. If the round-off errors are to be treated as uniformly bounded, but 
otherwise arbitrary, quantities, no condition beyond (2.3) and boundedness may be 
imposed on f(x) in the present section. 

In view of the last remark it is natural to base the arguments here on the finite 
trigonometric interpolation series for f(x), i.e., on 


N-1 


(4.1) f(mh) = > B, sin mrh, 


where N = x/h and 


(4.2) B, = (2/N) > f(oh) sin prk, 


It can be directly verified that the expression 


Nn-1 


(4.3) V(mh, nk) = >> B,E,* sin mrh 
r=] 


solves the difference problem (1.3), (1.4) for U(x, #) at the grid points x = mh 
t = nk,m = 1,2, ---,N— 1; = 1,2, --- , Hence, 


V(mh, nk) = U(mh, nk). 
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Insertion of (4.2) into (4.3) and an interchange of summations leads to 


(4.4) U(mh, nk) = > gomaf (oh) 


P=] 
with 
N-—1 


(4.5) Zomn = (2/N) >> E," sin mrh sin prh. 


r=1 


The quantity (4.5) will now be appraised under the assumption that (2.6) and 
hence (2.7) is true. Moreover, the restriction 


(4.6) O0<s<1 


will be imposed. Under this hypothesis EZ, is a decreasing function of r, for r < 
N — 1. Therefore, EZ," is a decreasing function of r for all odd m, while for even n 
the quantity E," decreases as long as E, > 0 and increases thereafter, for r < N — 1. 
In order to take advantage of this observation (4.5) is written in the form 


(4.7) Sonn = (1/N) Ecos (m — p)rh — (1/N) > E-* cos (m + p)rh, 


and each of the two right hand terms is subjected to a summation by parts. Recall- 
ing the identities 


N—1 N—1 Tr N-1 
» > a,b, — > { (on _ Gr) >> ,} + an py Ss b, 
r=0 r=0 


u=0 u=0 


sin ((r + 1)0@/2) cos (r@/2) 


spew sin 46 





ul _ sin (r + 3)0 
2, 008 0 = 


’ one finds that 


(4.8) (1/N) SE cos (m — p)rh 


= | —(1/N) | (Eras — E) > cos (m — ean} 


b=0 


N-1 


+ (1/N) Ex" >> cos (m — p)uh 


w=0 


N-1 
, ae | E — Er, | + | Es? | 
< r=0 


~# | sin $(m — p)h | 





Now, if 2 is odd, 
N—1 N-1 
DE — Eyl = 2 (EP — By) = BP — Ev $2, odd. 
r=0 r=0 


If E, is non-negative for 0 < 7; < WN the last formula is correct for even n as well. 





IMPLICIT DIFFERENCE APPROXIMATIONS 


But if Z, vanishes for = 1 one has, instead, with re = [ro], 
Nn-1 [r_) —1 
Lif -Bal= 2 (2 - By) +lBw- Bl 
r=0 r=0 
N-1 
+ 2 (Ea — EY) = ES - E+ |B - E|—- Bnt Ev <3, 
rerstl 


even 7. 


Inserting the result of the last two inequalities into (4.8) and taking account of the 
fact that both m and p range between 0 and N = x/h it follows that the first right 
hand term in (4.7) is numerically not greater than 4 | m — p |-*. The second sum- 
mation in (4.7) is similarly found to be bounded, numerically, by 


4h/(x | sin }(m + p)h |), 


where m and p range from zero to x/h. For m + p < 2/h that bound is clearly not 
greater than 4(m + p)-', which is less than 4 | m — p |-*. For 


m+ p = (2x/h) — y, 0<¥ < x/h, 


the bound 4y~" = ((2x/h) — m — p)~ is obtained, instead. But this last quantity 
is also not greater than 4 | m — p |-', since (2%/h) — m — p = ((x/h) — m) + 
((x/h) — p) = | ((x/h) — m) — ((x/h) — p) | = | p — m |. Collecting these 
results it is seen that (4.7) implies the inequality 


8 
| Sr - 
| Zomn | [m—pl 


On the other hand, 
| Lomn | < 2, 
because of (2.7) and (4.5). Combination of the last two inequalities yields 


10 


Let f be defined by 
f = lub. | f(x) |. 


OSzar 
Then (4.4) and (4.9) imply 


1 


N-1 
< Coe ee 
| U(mb, nk) | < 107 2) ESIC 


+ dx te dx 
< 107 {f went l ait 


< 10f(1 + 2 log (*/h)) < Cf log (1/h), 
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where C is a constant. This proves 
Tueorem 3. If f(0) = f(r) = 0, and | f(x) | <finO0 < x < xr;andif 
2(1 — 2s)k — hk? < 0, 0<s <1, 
then | U(mh, nk) | < Cf log 1/h, for0 < m < x/h, 0 < n. The constant C depends 
on hy only. 
As a consequence of this theorem the departure caused at a grid point by the 
cumulative effect of all errors at grid points below that ¢-level is less than 


e Cif log (1/h)k7, 


where ¢ is an upper bound for the error committed at one point. If h; is small C is 
about equal to 10. Round-off errors of this size are considered quite manageable in 
numerical work. 

It is an open question whether the result of theorem 3 is the best possible, i.e., 
if there exists a constant C* such that corresponding to arbitrarily small values of 
hand k satisfying (2.6) a function f(x) can be found for which 


| U(mh, nk) | > C*f log (1/h) 


at some point x = mh, t = nk. It is well known that the factor log (1/h) can be 
omitted if h, k and s satisfy the inequality 


s(li—s)k-—h? <0, 


for then the difference equation (1.3) satisfies a maximum principle. But this in- 

equality is more restrictive than (2.6). For s = 4, for instance, it implies the unde- 

sirable condition k/h? < 1. The boundedness statement proved in [5 ] iscompatible 

with values of U of the order 0(4-*) and is therefore weaker than that of theorem 3. 
5. A Lemma. 


1— (i-—s)z 
1 + sz 
a) Jf min(sz, —(1 — s)z) > —0, > —1, then 
| &.(z) | < cs, Cy @ constant. 
b) f0 <s < 1,2 > Oand (1 — s)z < 1, then 
®,(z) < }s*2*. 
c) If 2/2 < 4 < 1, then 
—c* < (2) < 0, C2 a constant. 


Lemma: Let ®,(z) = log + 2. 


Proof. a) It is readily verified that 


(5.1) &(s) = | eric if —(1 —s)s > —1, ss > —1. 


Under the assumption of part a) one has | 1 + r| > 1 — 6, in the interval of in- 
tegration. Hence, 


| (2) | < (a — 4) [ lelér < ost 
—s8)Z 
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b) Since —(1 — s)z < 0, sz > 0 it follows from (5.1) that 


#.(2) S Is + T 


c) This result is an immediate consequence of the identity 


dr< t rdr < 4(s2)?. 


a/e 2 
(2) = -2f to. 
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An Opes Formula for the Numerical Integration 


First Order Differential Equations (II) 
By Herbert S. Wilf 


In a previous paper [1 ], referred to below as J, a set of formulas was derived for 
the numerical integration of systems of first order differential equations. In what 
follows we will consider the questions of convergence and stability of the method 
and higher order formulas. 

We quote here, for reference, the final results of J, namely the propagation 
formulas: 


(1) =wt Sf (x, yo) + 8f(x1, 91) — f(x2, y2*) J 


(2) yo* = Syo — 4y1 + AL f(x, yo) + 2f (x1, 91) J 
for the solution of 
(3) y’ = f(x,y) 


where ¥o is given. 


1. Convergence. The solution of (1) and (2) for y: proceeds by taking as an 
initial guess the y.* from the preceding point. One then applies (2) and (1) succes- 
sively until consecutive values of y, agree to sufficient accuracy. 


Received 16 August 1957. 
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The convergence is governed by the following theorem. 
THEOREM. Suppose that on the interval [xo, x2], 


Hi: exists and 1s everywhere continuous. 
H2: | < M. 
oy 
Further suppose that the mesh interval h has been chosen so that 


: V2i —3 _ .7926 
HS: 0<4< 3 — = =. 





Then the iterative process defined by (1) and (2) converges to a solution y; 
of the equations (1), (2). 

Proof. Letting y; denote the solution of (1), (2) and 1, y;.r denote the r-th 
iterated values of y, and y2 respectively, we have 


yest = tH [Sf(n yo) + Bley a2) — Sey Hargs)] 
and therefore, 
(4) on — maven = (SCS Cen 90) — fl 914)] — Liles 92") — fl I) I). 
Now, . 
(5) len 9) — fen 92) = (1 — Ho) 2 (eum) 


for some m, in (1, yi.) and 


(6) f(x, yo") — f(X2, Voorn) = C92" — V3.r41) x, (x2, 2) 


~ x (x2, »)|C-40 — Yr) + 4h(91 — Yr) x (x1, »)| 


= 4(y1 — yr) x (x2, m) [a (x1, 1) — i] 


Inserting (6), (5) in (4) we get 


(7) (1 — M1941) = (91 — Yr) HAs x (x1, 1) 


" ao (xn »)| x beacead, 97 i]} 


for some 2 in (y2*, ¥3,-41)- 
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To prove convergence, i.e., that (y1 — 1) —> 0, we need only show that the 
coefficient of y: — y1,, in (7) is, in absolute value, < 1. We have 


lyn — Y1r41 1 S| 91 — ner | = (8M + 4M(hM + 1) 


& 
12 


2 2 
=I — ns | (eae +22"), 


Hypothesis H3 assures us that a(t + Mf) < 1 which proves the theorem. 


2. Stability. Let y denote the true solution of (3), z denote the calculated solu- 
tion from (1), (2) and 7 = y — z be the error. 
Now y satisfies 


(8) y= y+ be [Sf (xo, y1) + 8f(x1, 91) — f(x2, y2*) ] + Ti 


(9) ye* = Syo — 4y1 + 2h (xo, yo) + 2f(x1, 71) J 


where 7; is the truncation error given by (8.1) of J, and z satisfies 
(10) a1 = so + [5G (a to) + f(s, 1) — fen, 2*)] + hn 


(11) Za* = S2o — 42, + 2h[ f(xo, 20) + 2f(x1, 21) ] + Apr 


where p, p2 are round off error in the evaluation of f. 
By subtraction, we get 


292 292 
(12) moss {1 — hg + EI = m4 -Ye\ er 


where we have assumed 


a) ¥ 


F is constant (=g) over the interval in question. 


b) Round off error is negligible. 
c) Truncation error is constant (=T7). 


The solution of (12) is 
(13) in = A"No + (x* ~- 1) 


where 
1 — (hg)?/6 
ei The + Hey 
The initial error 4 will therefore decrease in magnitude when | \ | < 1 which 
occurs for — © < hg < 0 and for hg > 2. 
We may then summarize the discussion of this and the preceding sections by 
noting that convergence and stability occur simultaneously only when 


(15) —.7926 < hg < 0. 

If af/ay is positive, errors introduced at any stage will grow in magnitude, and 
the integration had better be performed in the reverse direction. This behaviour is 
characteristic of numerical integration techniques. 


y T 
hg(1 — hg/2) 
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3. A Higher Order Formula. The formulas corresponding to (1) and (2) for an error 
of order # are as follows: 


(16) 91.= yet Pe [OVC 90) + 197 eu, 1) — Sflens 924) + Slee 99")] + OC8) 


(17) 93" = 0 + 5 Cleo 90) + Afl2n 90) + Sen 12") 


(18) ys* = 9y1 — 8y0 — ShL f(x0, yo) + 2f(x1, 91) — f(x2, y2*) J- 


These formulas are used to find ; as follows: 
a) Guess yi = y2*(x — h), yo* = ys*(x — h). 
b) Calculate improved values of y3*, y2*, yi in that order. 
c) Repeat from b) until y; has converged before proceeding to the next 
point. 
Note that no starter formulas are required since initially we may guess yo = 
¥1 = y2* = ys* and proceed from b) above. 
Nuclear Development Corporation of America, 
White Plains, New York 


1. Hersert S. Wixr, “An open formula for the numerical integration of first order differential 
equations,” MTAC, v. 11, 1957, p. 201-203. 
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Note on the Computation of the Zeros of Functions Satisfying 
a Second Order Differential Equation 


By D. J. Hofsommer 


It has been pointed out by P. Wynn [1] that, if a function satisfies a second 
order differential equation, this fact may be used with advantage in the computa- 
tion of its zeros. In his note he only pays attention to Richmonds formula which, 
incidentally, was already known to Schréder [2]. We will elaborate his idea to con- 
struct another iteration formula. 

Let f(x) be the function, the roots of which are to be computed. Let a be sucha 
root and let x be a first approximation. If the approximation is sufficiently close, 


(1) a=x— f/f — ¥(f"/f) (f/f)? — 3 BFP)? — FP S/F) * + OL(F/P") I. 


This series may be used either for direct computation in taking enough terms, or 
for obtaining an iterative process if only few terms are retained. If f(x) satisfies the 
homogeneous differential equation 


(2) Sf" =2Pf'+Qf+25, 
substitution in the series (1) yields 
(3) a=x—f/f' + (P+ S/f’) (f/f)? 


— 3(4P? — P’+Q + 10PS/f' + S'/f' + 6S/f*) (f/f')* + OL(F/f')*] 
Received 3 May 1957. 
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a=x— f/f — (P+ SP) (S/F + CFP) 


as = 1 3 





From (4) and (5) it follows that the iterative processes 
(6)  Xegn = He — f (Xx) /f (xe) — LP (me) + S(xe)/f (xe) IS (xe) /f (me) F 


or 
(7) Sey = Xe — : 
f (xx) /f (xn) — P(x) — S(xe)/f’ (xe) 


are cubically convergent. 
Wynn’s process gives 





1 
SF )/f) — Pl) — Sl) /F ee) — FOCI) 
if f(x) satisfies (2). 
From equality (7) it is seen that the last term in the denominator may be left 
out without affecting the order of the iteration process. 


For comparison we have calculated by aid of (7) the first three zeros of Jo with 
one iteration step. The results are 


Xo x1 x, (Wynn) 


2.405 2.404 825 557 694 2.404 825 557 697 
5.520 5.520 078 110 286 5.520 078 110 286 
8.654 8.653 727 912 904 8.653 727 912 914 


Of course one iteration step amounts to using the development (4) or (5). If in 
such a case the required accuracy is almost, but not quite, obtained, application 
of (3) is preferable to taking a second iteration step. 

In such a way we have calculated the zeros of P3; (cos ¢). Using Smith’s table 
[3] the zeros yg» were found in hundredths of a degree. Then new values for gn 
were calculated by aid of (3) or (4) using the NBS tables [4]. In this way inter- 
polation was avoided. We found 


Xe+1 = Xe 





m ¢m(radians) ¢m(radians) 


-570 796 326 8 


3 259 745 5 
9 491 725 5 


-48 
-40 


1 
1 
1 
1.31 

1.235 724 047 9 

1.151 956 859 3 

1.068 190 338 7 

0.984 424 715 0 0 

0.900 660 291 7 0.064 126 781 0 


COnAME WHO 


The approximation in hundredths of a degree could have been obtained also by 
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starting with the approximation for the zeros of P,(cos ¢) 
on * (3 = a n odd 


Cn = (3 - ym te, n even 


and applying one Newton step xx41 = x, — f(xx)/f’ (xx). 
Finally we remark that if f(x) is a member of a set f,(x), it often happens that 
fn’ (x) satisfies a recurrence relation 


fn! (x) = A(x)fas(x) + B(x) f(x). 


It then may be useful to transform 
Sn! (xt) /fa(%%) = A (Xe) fare) /fa(xe) + B(xx), 


especially if f,(x) is tabulated and f,’(x) not or if the tables of f,(x) are more 
extensive than those of f,’(x). 


Mathematical Centre, 
Computation ment, 
Amsterdam, Netherlands 


1. P. Wynn, “On a cubically convergent process for determining the zeros of certain func- 

- MTAC, v. 10, 1956, p. 97-100. 
. E. SCHRODER, “Ueber unendlich viel Algorithmen zur Auflésung der Gleichungen,”” Math. 

Pm... v. 2, 1870, p. 317-363. 

3. E. R. Smite & ARCHIE Hicpon, ‘Zeros of the Legendre polynomials,” Jn. Science, lowa 
State College, v. 12, 1938, p. 263-274. " (MTAC, Rev. 132, v. 1, 1944, p. 153-514.) 

4. NBSCL, Table of Sines and Cosines to Fifteen Decimal Places at Hundredths of a Degree 
(NBS Applied Math. Series, no. 5), Govt. Printing Office, Washington, 1949. (MTAC, Rev. 662, 
v. 3, 1949, p. 515-516.) 


A New Mersenne Prime . 
by Hans Riesel 


On September 8, 1957, the Swedish electronic digital computer BESK estab- 
lished that the Mersenne number Mj; = 2*!7 — 1 is a prime. The result was re- 
calculated on September 12. The method used for testing this number was a Lucas’ 
test, and the running time was 5’ 30". The new prime has 969 decimal digits and 
on the BESK was found to be 

23217 on 1 


2 59117086 013 20262 77762467 67922441 53094181 888 75531 254 27303 
974 92316 18740192 665 86362 086 20120 951 68004 83406550 695 24173 319 41774 
416 89509 238 80701 74103777 09597512 042 31306 66240829 16353517 952 31118 
615 48622 65604547 691 12759 58487756 105 68757 93119101 771 14088 262 52153 
849 03583 04011850 721 16424 74746182 303 14713 98340229 28807454 567 79079 
410 37288 235 82070 589 23510 684 33882 986 88861 665 86502 80927692 080 33960 | 
586 93087 905 00409 503 70987 590 21190 18371991 62099400 256 89351 131 36548 | 
829 73911 26567973 03241986 51725011 641 27035 09705427 77347797 234 98216: 
764 43446 668 38311 932 25400 996 48994 051 79024 16240565 19054483 690 80961 
606 16257 43042361 721 86333 941 58524 264 31208 737 26659 196 20617 535 35748 
892 89459 96291951 83082621 860 85340 093 79328 394 20261 866 58614 250 32514 
507 73096 274 23537 682 29386 49407127 700 84607 71242118 23080804 139 29808 
705 75047 138 25264 571 44837 93711250 320 81826 126 56664 90842516 994 53951 
887 78961 365 02484 057 39378 59459944 433 52311 882 80123 66040626 246 86092 
121 50349 937 58478 229 22371 443 39628 85848593 82157388 21232393 687 04616 
067 73629 093 15071 


Received 8 November 1957. 
Ormiingsgatan 67, 
Stockholm, Sweden 
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REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 
1[C, D].—J. Peters, Ten-place logarithm table: 

v.1. Ten-place logarithms of numbers from 1 to 100,000, together with Appen- 

dix, Mathematical Tables, edited and computed by J. Peters and J. Stein (with 

22-place logarithms of the trigonometric functions by G. Witt), 

v. 2. Ten-place logarithms of the trigonometric functions from 0° to 90° for every 

Thousandth of a Degree, 


Auxiliary Tables for Ten-place logarithms table, 


Published by the Reich Office for Land Survey by J. Peters, Under Scientific 
Supervisor, Berlin, 1922, 1919, 1919, New Revised Edition, with English trans- 
lations, F. Ungar Publishing Co., New York, 1957, xvi + 607 + xxciii + 195, 
vii + 901, 76 + 53 p., 27 cm. Price $50.00 per set. 


This is a very welcome (photographic) reprint of volumes long out of print, 
which have been much sought after, perhaps more for the Appendix, familiarly 
known as Peters-Stein, than for the main tables. The reproduction is excellent, the 
binding seems durable, and there is an English translation by C. J. Hyman of the 
German introductory and explanatory text. A supplementary list of the errors 
which have been noticed since the original printing has been prepared by C. J. 
Hyman; there is even a translation of this into German. 

The reviewer would like to point out that three errors which are associated with 
him (p. 51 of the Auxiliary volume) were actually found by the late C. R. Cozens 
and by Dr. J. W. Wrench, Jr. They were originally pointed out in this journal [1], 
and this fact was noted by the reviewer in his communication with the publisher. 

According to R. C. Archibald [2], J. T. Peters (1869-1941) was perhaps the 
greatest mathematical table maker of all time. An astronomer by training, he was 
largely responsible for the official German nautical and astronomical tables. In 
addition, he prepared a long series of mathematical tables, at varying numbers of 
decimal places, all of the highest standard of accuracy and of format. 

For those who want a set of logarithmic tables, of upper middle class precision, 
with trigonometrical functions for arguments in degrees and decimals, the present 
volumes are unequalled. They include the following: 

In volume I there is a table of logarithms to 10D of the integers from 1 to 1,000 
and from 10,000 to 100,000, occupying about 600 pages. First differences are tab- 
ulated. 

In volume II there is a 10D table with positive characteristics, of log sin x, log 
tan x, log cot x, log cos x, for x = 0(.001)45°. First differences are tabulated, except 
for small angles, where other methods of interpolation are to be used (see b) be- 
low). These occupy 900 pages. 

In the third volume there are auxiliary tables for use with 1), 2), which occupy 
about 70 pages. This volume also includes the translations of the textual matter, 
which occupies 53 pages. 

a) Tables for improvement of first difference, given the second difference 
(0(10) 1000) and the phase (0(.01)1) (20 pages). 
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b) Tables of S = log sin x — log x and T = log tan x — log x, with first dif- 
ferences, for x = 0(.001)2.1°, to be used for interpolation (21 pages). 

c) Miscellaneous relevant conversion tables (22 pages). 

The various introductions give some account of the sources and construction 
and checking of the tables; there are worked examples to show the uses of the 
main and auxiliary tables. Included are Grimpen’s values of logiop to 84D and 
log.p to 82D, for primes up to 113. 

In addition to these tables there is, in Volume I, an Appendix (of about 200: 
pages) containing the basic numerical data which Peters and his collaobrators used 
in their work. This collection, made in collaboration with J. Stein (which is supple- 
mented by a 22D table of the logarithms of the trigonometrical functions due to 
G. Witt), is invaluable to all serious computers. We shall describe its contents in 
some detail. 

1) Many decimal values of such constants as x, e, y, M = loge, together with 
various powers and multiples (12 pages). 

2) Exact values of the first ten powers of the integers up to 308 with certain 
higher powers of the primes under 100 (23 pages). 

3) Reciprocals, to 32D, of the first sixteen powers of the integers up to 100 (22 
pages). 

4) Exact values of ! for m up to 60, and their expressions in prime factors 
(2 pages). 

5) Reciprocals, to 54D, of n! for m up to 43. Also e, e~, sin 1, cos 1, sinh 1, cosh 1 
to 52D (1 page). 

6) Common logarithms of m!, to 18D, for up to 1200 (8 pages). 

1) for » up to 60,0 < k < n, and 
their expressions in prime factors (14 pages). 

8) Exact values (as rational fractions in lowest terms) of the Bernoulli num- 
bers B,, for m up to 90. Also logi B, to 10D and B,/B,. to 4S for m up to 250 
(S pages). 

9) Exact values of the tangent numbers 7, and the Euler numbers £,, for 
n up to 30. Also logi 7, and log B, to 10D, T,/T,-1 and B,/B,1 to 4S, for n 
up to 50 (2 pages). 

10) Various sums of reciprocals, all to 32D: 2-* + 3-* +4" +5" + --- and 
2*+4°*+6"+--- uptom = 100;4*% + 6* + 8* + --- and3* +5" + 
7* + --- upton = 50;3* —4*+ 5" — 6%+ --- and3* —5*+7>" - 
9 + --- upton = 53 (5 pages). 

11) Early coefficients, to 24D, of terms in power series for direct and logarithmic 
trigonometrical functions of ax (or related auxiliary functions such as tan 4x — 
4/x{(2 — x)“ — (2 + x)}) (3 pages). 

12) Complete factorization of integers up to 10192 (29 pages). 

13) Various many decimal natural logarithm tables: log, », to 48D, for integers 
up to 146 and thereafter for primes up to 9973; radix table of log, (1 + a-10~*), 
to 48D, for a = 1(1)9, = 4, 5, 6; log, 2, log, 3, log. 5, log. 7 to 272D (26 pages). 

14) Various many-decimal common logarithm tables: radix table of logi 
(1 + a-10™), to 28D, for a = 1(1)9, m = 1(1)14; logio m for m = 2(1)9, to 28D; 


7) Exact values of binomial coefficients ( 
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logy ”, to 61D, for integers m up to 100 and thereafter for primes up to 1097; radix 
table of logw (1 + a-10), to 61D, for a = 1(1)9, m = 3, 4, 5, 6 (10 pages). 

Witt’s Tables. 1. log sin a, log cos a, log cos a, log tan a, to 22D, for a = 0(10’)45° 
(positive characteristics) . 

II. log sin 8, log cos 8, log tan 8, for 8 = 0(1’)10’ and T = log tan 8 — log 6” 
together with the first three differences of log cos 8 and T. It is shown how to get 
21D values of the trigonometrical functions and their logarithms, at 1’ interval, 
using these tables and the addition formulae (33 pages). 


Joun Topp 


California Institute of Technology 
Pasadena, California 


1. Errors noted by J. W. Wrencu, MTAC, v. 2, 1946, p. 138 and v. 6, 1952, p. 255, and by C. 
R. Cozens, v. 1, 1943, p. 57. 


2. RAYMOND CLARE ARCHIBALD, ‘‘Mathematical Table Makers. Portraits, Paintings, Busts, 


Monuments, Bio-Bibliographical Notes,” Scripta Math. Studies, No. 3, Scripta Mathematica, 
New York, 1948, p. 68. 


2[F ].—A. GLopEN, Table de factorisation des nombres N* +- 1 dans l’intervale 20000 
< N < 30000, 100 p., 27 cm., hand-written, placed in the UMT Fine. 
[F]}—.A. GLopEn, Table de factorisation des nombres N* + 1 dans l’intervalle 30000 
< N < 40000, 100 p., 29 X 21 cm., deposited in the UMT Fine. 


These manuscripts considerably extend the tables of factorizations of N* + 1. 
All odd prime factors less than 800 000 are listed. If the remaining factor is less 
than 64 X 10" it is prime and the factorization is complete. The primes greater 
than 800 000 thus found are listed in the volume 30 000 < N < 40 000. Factoriza- 
tions of N* + 1 have been published for N < 6000, with many entries blank or 
incomplete. [For m < 1000, see MTAC, v. 2, 1947, p. 252. For 1000 < N < 3000 
see MTAC, v. 2, 1947, p. 211; v. 3, 1948, p. 118. For 3000 < N < 6000 see MTAC 
v. 4, 1950, p. 224; v. 11, p. 274. For 6000 < N < 10 000 see MTAC, v. 5, 1951, 
p. 28. For 10 000 < N < 20 000 see MTAC, v. 9, 1955, p. 26. For 30000 << N< 
40 000 see MTAC, v. 6, 1952, p. 102; v. 7, 1953, p. 33.] 


J. L. SELFRIDGE 


University of California 
Los Angeles, California 


3[F].—A. Scuoiz & B. SCHOENEBERG, Einfiihrung in die Zahlentheorie, Sammlung 
Géschen Band 1131, Walter de Gruyter and Co., Berlin, 1955, 128 + 20 p. 
15 cm. Price DM 2.40. 


This is a second edition of the booklet of the same title, by A. Scholz alone, 
which was first published in 1939, and which is a compact source of information 
in the subject of elementary number theory. Scholz died in 1942 and the new edi- 
tion was prepared by B. Schoeneberg. In the new version some misprints and errors 
have been corrected and the material has been brought up-to-date, e.g., recent 
numerical work on Mersenne numbers has been included. However, the chapter on 
“Algorithmisches Rechnen”’ has been omitted altogether; this is a great loss to 
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those interested in computation, for it contained accounts of many useful ideas in 
this area, in which Scholz was a real expert. The booklet is divided into four chap- 
ters: Divisibility properties; Congruences and residue classes; Quadratic Residues; 
Quadratic Forms. 


OLGA TAussKY 


California Institute of Technology 
Pasadena, California 


4[G ].—Emma Leumer, Roots of the reduced quintic, sextic, and seventh degree period 
equations, 22 p., 37 XK 28 cm., tabulated from cards punched by SWAC. De- 
posited in the UMT Fine. . 


These tables list the generalized Gauss sum S, = )>2=} exp (2 4 im*/p) for each 
prime p = kf + 1. This sum is the principal root of the reduced period equation of 
degree k. The other k — 1 roots are computed using similar sums of cyclotomic 
cosines [1], and are listed. The tables extend for k = 5 to p = 10061, for k = 6 
to p = 4093 and for k = 7 to p = 5153. The evidence for and against various 


conjectures about the distribution of S; is discussed in [1]. 


J. L. SELFRIDGE 


University of California 
Los Angeles, California 


1. Emma LeEaMer, “On the location of Gauss sums,” MTAC, v. 10, 1956, p. 194—202. 


5[K].—J. Arrcuison & J. A. C. Brown, The Lognormal Distribution, University of 
Cambridge, Department of Applied Economics Monographs, 5, Cambridge 
University Press, Cambridge, 1957, xviii + 176 p., 25 cm. Price $6.50. 


This useful monograph gives a thorough discussion of the history, theory, and 
applications of the lognormal frequency function, bringing together material 
hitherto to be found only in widely scattered publications. Though the authors’ 
principal motivation lay in applications to economic research, applications to other 
fields have not been slighted. Four new tables are contained in the appendix. The 
first gives the coefficient of variation, the skewness, the kurtosis, the ratio of the 
mean to the median, the ratio of the mean to the mode, the frequency up to the 
mean, and the Lorenz measure of concentration all to 4D or 4S for the standard 
deviation, ¢, of the normally distributed variable = 0(.05)1(.1)4. The next two 
tables give values of two functions defined by the series, 


ws (m — 1)*# 
a(t) “1+ 2G - 1)(m + 1) --- (n — 3 + 2j)j! 





a(t) = Yu(2t) — v,(* a 1) 


n—1 


which are required in finding minimum variance unbiased estimates of the mean 
and variance of the lognormal distribution from a sample of n. These are tabulated 
to 4D for nm = 10 (10) 200 (100) 500, 750, 1000, 5000, «© for both functions; for 
# = .05 (.05) 2 for y,(#), and for = .05 (.05) 1 for x,(#). The fourth table lists 
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solutions of the equation u* + 3u = k, required for moment estimation in the three 
parameter case, to 4D for k = 0(.2)10(1)24. The corresponding estimates of the 
variance are also given to 4D. An additional pair of tables gives for 65 artificial 
samples of sizes 32, 64, 128, 256, and 512 the results of various methods of estima- 
tion of the mean and variance. 


C.c4.¢ 


6[K].—Smney SreceL, Nonparametric Statistics for the Behavioral Sciences 
McGraw-Hill Book Co., New York, 1956, ix + 312 p., 23.7 cm. Price $6.50. 


This useful text contains twenty-one tables in an appendix. None of these are 
new but, like the text which brings together in one place discussions and illustra- 
tions of a good selection of the more important nonparametric tests which previ- 
ously had been available from rather scattered sources, this appendix makes the 
associated tables readily available. The sources from which these tables are either 
copied or adapted are given in each case. 

In addition to tables of areas under the normal frequency curve, percentage 
points of the Student-Fisher ¢ and x?, factorials, binomial coefficients and squares 
and square roots of the kind ordinarily found in statistics texts, the following are 
more special: 


Table D. The cumulative binomial probabilities, 


>(*) (4)* to 3D for x = 0(1)N and N = 5(1)25. 


Table E. 1, 5, 10, 15 and 20 percentage points for the Kolmogorov-Smirnov 
one-sample statistic D to 3D for N = 1(1)20, to 2D for N = 25, 30, 35 with a 
large sample approximation formula for N > 35. 

Table F. Upper and lower 5% values for the number of runs in sequences of 
alternatives for the numbers of the two kinds of elements, m, m2 = 2(1)20. 

Table G. Critical values of T in the Wilcoxon matched-pair signed-ranks test 
at the .025, .01 and .005 levels for the one-tailed test (.05, .02 and .01 for the two- 
tailed test) for the number of matched pairs N = 6(1)25. 

Table H. A page of the tests on the median due to J. E. Walsh with one and 
two-tailed significance levels to 3D for sample sizes N = 4(1)15. 

Table I. For the 2 X 2 contingency table with row entries A, B and C, D using 
the Fisher-Yates test, .05, .025, .01, and .005 critical values for C given A (or for 
D given B) for C+ D = 2(1)A + Band A + B = 3(1)15. 

Table J. Probabilities to 3D of values of U as small as that observed in the 
Mann-Whitney test for the sample sizes m < m = 3(1)8. 

Table K. Critical values of u in the Mann-Whitney test for the significance 
levels .001, .01, .025, .05 for the one-tailed test (.002, .02, .05, .1 for the two-tailed 
test) for m = 1(1)20 and m,. = 9(1)20. 

Table L. Critical values of Kz in the Kolmogorov-Smirnov two-sample test at 
the .01 and .05 levels for the one-tailed and two-tailed test for the sample sizes 
m = m, = 3(1)30, 35, 40. 

Table M. Approximate formulae for critical values of D in Kolmogorov-Smirnov 
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two-sample test at the .001, .005, .01, .025, .05, and .1 significance levels for the 
two-tailed test for large samples. 

Table N. Probabilities to 3 or more D for x,? as large as that observed in the 
Friedman two-way analysis of variance by ranks for N sets of K matched samples. 
For K = 3 and N = 2(1)9 values are given for selected values of x,? covering the 
effective range in each case. For K = 4 values are given for similarly selected 
xr’’s for N = 2, 3, 4. 

Table O. Probabilities to 3D for selected values of H as large as that observed 
in the Kruskal-Wallis one-way analysis of variance by ranks for all sets of sample 
SIZES i, M2, Mz UP tO Mm, = M2 = nz = 5S. 

Table P. Critical values of the Spearman rank correlation coefficient to 3D at 
the .01 and .05 levels for the one-tailed test for sample sizes, N = 4(1)10(2)30. 

Table Q. Probabilities to 3D or 2S for values of the Kendall rank correlation 
coefficient S as large as that observed for S = 0(2)36 with sample size N = 4, 5, 
8, 9 and for S = 1(2)45 with NV = 6, 7, 10. 

Table R. Critical values of s to 1D in the Kendall coefficient of concordance at 
the .01 and .05 significance levels for k sets of N rankings. For N = 3, values are 
given for k = 8, 9, 10, 12, 14, 15, 16, 18, 20; for N = 4 for k = 4, 5, 6, 8, 10, 15, 
20; for N = 5, 6, 7 for k = 3(1)6, 8, 10, 15, 20. 

Se {a< 


7[K].—NBS Applied Mathematics Series, No. 48, Fractional Factorial Experiment 
Designs for Factors at Two Levels, U. S. Gov. Printing Office, Washington, D. C., 
1957, iv + 85 p., 26 cm. Price $.50. 


This table gives designs for experiments with from 5 to 16 factors, each of 
which may occur in two levels. These involve fractional replication ranging from 
1/2 to 1/256. The designs are given in terms of blocks. For each design we are told 
which main effects are measurable, which two-factor interactiona are measurable 
for blocks only, and also for blocks and rows. In addition we are told which com- 
binations are block confounding and which are row confounding. The number of 
designs given is so large that it is to be expected that an experimenter will find 
exactly the design needed for any desired information on two level factors involv- 
ing not more than 16 factors. 

MARSHALL HALL, JR. 


The Ohio State University 
Columbus, Ohio 


8[K ].—D. Tercuroew, ‘Tables of expected values of order statistics and products 
of order statistics for samples of size twenty and less from the normal distribu- 
tion,” Ann. Math. Stat., v. 27, 1956, p. 410-426. 


These tables represent the end product of a project undertaken originally by 
the Institute for Numerical Analysis, National Bureau of Standards, to extend, 
both in terms of sample size and decimal place accuracy, existing tables of the 
means, variances, and covariances of order statistics for samples from a normal 
population. Tables I and II give expected values of the order statistics and their 
products to 10D for samples of size m = 1(1)20. Tables III and IV give certain 
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auxiliary functions required in the calculation, and in any further extension to 
larger sample sizes, to 25D. Tables I and II to 20D can be obtained on 925 punched 
cards by writing the Director, Numerical Analysis Research, University of Cali- 
fornia, Los Angeles 24, California, and requesting table 5024. A charge of three 
new cards for each one used will be made. 


C. A. BENNETT 


Hanford Laboratories Operation 
General Electric Company 
Richland, Washington 


o[K ].—P. K. Baatracuarya, “Joint test for the mean and variance of a normal 
population,” Calcutta Stat. Assn., Bull., v. 6, 1955, p. 73-90. 


Let x, ---, x, be a sample from N(m, o*) and let n@ = >> x;, ns¢ 
= >> (xi — mo)*, (n —1)8? = >> (x; — 2)%, st = (& — me) V/n. The problem of 
testing simultaneously Hi:m = mo against m > mo, and He: ¢ = oo against ¢ > oo 
is regarded as a four-decision problem. We reject Hi when ~ > c, and H; when 
s* > k, where c and k are fixed by requiring that Pr (reject H;) = Pr (reject H2) 
and Pr (accept both) = 0.95, when m = mo and o = oo. A second procedure is 
obtained by replacing Z by ¢ and s* by s,*. For both procedures we are given to 
3D the probabilities of the four decisions, for m = 5(5)20, o?/o0? = 1(.5)3.5, and 
—/n(m — m) = 0, .1, .4, 1, and 1.5. 


J. L. Hopcgs, Jr. 


University of California 
Berkeley, California 


10[K ].—P. B. Patnaik, ‘Hypotheses concerning the means of observations in 
normal samples,” Sankhyd, v. 15, 1955, p. 343-372 (appendix). 


The non-central ¢’ distribution is given by 


where Hh, (x) - | : exp 34(u + x)%du. The solution of P(#’ < x) = .05, 
0 PV 


P(t! > x) = .05 is tabulated to 2D for » = 4, 5, 6 = 0(.5)4, and to 3D for »y = 
6(1)15, 20, 6 = 0(.5)4. The tabulation is based on the tables by Johnson and 
Welch [1]. An approximation for the upper and lower percentage points is given 
in terms of the central ¢-distribution. A number of statistical hypotheses for which 
the ?’ statistic can be used are discussed. 


INGRAM OLKIN 


Michigan State University 
East Lansing, Michigan 


1. N. L. Jonnson & B. L. Wetca, “Applications of the non-central ¢-distribution,” Bio- 
metrika, v. 31, 1940, p. 362-389. 
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11[K].—F. N. Davi & N. L. Jounson, ‘Some tests of significance with ordered 
variables,”” Roy. Stat. Soc., Jn., ser. B, v. 18, 1956, p. 1-31. 


A survey of test procedures based on medians, quartiles, interquartile ranges, 
etc., is made. Numerous small tables, computed by the method of moments, of the 
approximate percentage points are presented. In particular, to test the hypothesis 
of equality of the means of two normal universes with equal but unknown variance, 
the test criterion T= (Xu, — Xu,) (Ni + N2)/Ni(Xu, — Xi.) + No(Xu, — Xi) 
is proposed in which Ni, N2; X u,, X uz; Xu,, Xuz3 X1,, Xi, are the sample sizes, the 
medians, and the upper and lower quartiles in samples from the two universes 
respectively. 1%, 2.5%, and 5% points for T are given to 3D for Ni, Ne = 11(4)51. 


I. R. SAVAGE 


Center for Advanced Study 
in the Behavioral Sciences 
Stanford, California 


12[K}].—C. I. Butss, W. G. Cocuran, & J. W. Tuxey, “A rejection criterion 
based upon the range,” Biometrika, v. 43, 1956, p. 418-422. 


The data consists of k independent samples of m measurements each. Under the 
null hypothesis, the measurements are assumed to be normally distributed and 
have the same variance. The alternative hypothesis of interest is that one of the 
measurements is too large or too small to be consistent with the null hypothesis. 
That is, the problem is that of rejecting an outlying observation. To perform the 
proposed test first compute the range for each sample. Let T equal the largest of 
the k ranges divided by the sum of all the ranges. If T is significantly large, the 
null hypothesis is rejected and the aberrant observation determined by inspection 
and rejected. Table 1 contains 3D upper 5 per cent points of T for k = 2(1)10 
and m = 2(1)10. Table 2 contains 2D upper 5 per cent points of (k + 2) T fork = 
10, 12, 15, 20, 50, and = 2(1)10. 

J. E. Wats 


Military Sppeetions Research Division 
Lockheed Aircraft Corporation 
Burbank, California 


13[K].—H. A. Davnn, “Revised upper percentage points of the extreme studentized 
deviate from the sample mean,”’ Biometrika, v. 43, 1956, p. 449-451. 


The upper 100 a percentage points of the extreme studentized deviate in a 
random sample of n observations, based on an independent estimate of the variance 
with »v degrees of freedom, are given to 2D for a = .10, .05, .025, .01, .005; and to 
1D for a = .001, for m = 3(1)12 and » = 10(1)20, 24, 30, 40, 60, 120, «©. Accuracy 
to one unit in the last place recorded is claimed. This is a revision and slight exten- 
sion of the original tables by Nair [1]. 


C. A. BENNETT 
Hanford Laboratories Operation 
General Electric Company 
Richland, Washington 


1. K. R. Narr, “Tables of percentage points of the ‘Studentized’ extreme deviate from the sam- 
ple mean,” Biometrika, v. 39, 1952, p. 189-191. [RMT 1182, MTAC, v. 8, 1954, p. 84.] 
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14[K ].—Cuartes W. Dunnett, “A multiple comparison procedure for com- 
paring several treatments with a control,” Amer. Stat. Assn., Jn., v. 50, 
1955, p. 1096-1121. 


This paper is concerned with the case where the numerical results of an experi- 
ment performed to compare p treatments with a control can be summarized in the 
form of a set of numbers Zo, 71, --- , ,, and s, where the Z’s are the means of the 
p + 1 sets of observations which are assumed to be independently and normally 
distributed, Z referring to the control, and s is an independent estimate of the 
common standard deviation. A procedure is presented for making confidence state- 
ments about the true values of the differences #; — Z, the procedure having the 
property that the probability of all » statements being simultaneously correct is 
equal to a specified value P. Tables are available that enable one to use the pro- 
cedure for P = .95 or .99, p = 1(1)9, and degrees of freedom = 5(1)20, 24, 30, 
40, 60, 120, «. When the number of observations in each set are equal the tables 
yield exact probabilities for the 3D values given for either one-sided or two-sided 
confidence limits while if the number of observations are not equal the value of P 
must be approximated. Illustrative examples are given along with a theoretical 
discussion of the procedure. A section dealing with the problem of optimum alloca- 
tion of observations between control and treatments is also included in the paper. 


C. F. Kossack 
Purdue University 
Lafayette, Indiana 


15[K ].—Geratp J. LizepeRMAN & GeorGE J. Resnikorr, “Sampling plans for 
inspection by variables,” Amer. Stat. Assn., Jn., v. 50, 1955, p. 457-516. 


Inspection procedures are known as sampling inspection by variables when 
items of a sample drawn at random from a lot are classified into defective or non- 
defective on the basis of a variable quality characteristic; the decision to accept 
or reject the lot is a function of these measurements (as opposed to the decision 
being based on the number of defectives as in sampling by attributes). Measure- 
ments of this quality characteristic are assumed to be independent, identically 
distributed normal random variables. 

Sampling plans are presented giving the sample size » and the maximum allow- 
able estimated proportion of defectives ». They are classified according to Accept- 
able Quality Level (AQL) and sample size. Tables are presented for the purpose 
of transforming the data into uniform minimum variance unbiased estimates of 
the percentage defective in the lot. 

Table I presents inspection plans based on a known standard deviation for 
AQL’s = .04, .065, .1, .15, .25, .4, .65, 1, 1.5, 2.5, 4, 6.5, 10, 15 all in per cents. 
Minimum values of (U — L)/o are given below which no inspection is necessary 
in order to reject the lot. U and L denote upper and lower specifications of the 
quality characteristic and ¢ denotes the known standard deviation. Table II gives 
the estimates of the lot percentage defective to 3D for 


Gia San <(. *—)' and Cx! = z= Hf “Y= 0(.01)4 


co n—1 
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in which o is the known standard deviation. Table III gives sampling plans based 
on an unknown standard deviation for the same set of AQL’s and Table IV gives 
estimates of the lot percentage defective to 2D for C, = (U — X)/Sand C, = 
(L — X)/S = 0(.1)3.9 where S* is the unbiased estimate of the variance. Tables V 
and VI are strictly analogous tables for use with the average range, 2, the argu- 
ments in Table VI being Cy” = k(U — X)/R and C,” = k(X — L)/R with the 
values of h given to 3D in Table V for each sampling plan. Figures 1 up to 16 
present OC curves for all the sampling plans based on unknown standard deviation 
(curves for sampling plans based on average range and known standard deviation 
are essentially equivalent) . 
Cyrus DERMAN 


Columbia University 
New York, N. Y. 


16[K ].—STANLEY REITER, “Estimates of bounded relative error for the ratio of 
variances of normal distributions,” Amer. Stat. Assn., Jn., v. 51, 1956, p. 
481-488. 


Let m1, m2 be sizes of independent random samples, respectively. Let @ = m:/me. 
Enter Table 1 with @ = 1, 1.5, 2, 5, 10, 100, and for a required bounded relative 
error K (in the sense of Blackwell and Girshick [1] K = 1.01, 1.02, 1.05, 1.10, 
1.25, 1.30, 1.40, 1.50) find 8* to 4D to be used in the “best” estimate 8*(é,7/6,") of 
o:?/o2?, where m01"/o;? and m262*/o:? are required to be independently distributed 
as x? with m and m, degrees of freedom. The confidence coefficient a associated 
with the estimate is then read from one of seven charts which show a as a function 
of m, m, and K (K as above, excluding K = 1.01). 


H. A. FREEMAN 


Massachusetts Institute of Technology 
Cambridge, Mass. 


1. Davin BLacKweE.Lt & M. A. Grrsuick, Theory of Games and Statistical Decisions, John 
Wiley and Sons, Inc., New York, 1954, p. 316-323. 


17[K ].—Epcar J. Gipert, “The matching problem,” Psychometrika, v. 21, 
1956, p. 253-266. 


Each of two decks of cards consists of s suits of c cards per suit. One deck is 
thoroughly shuffled, and each card is then compared with the card in the same 
ordinal position in the other deck. If the two cards compared are of the same suit 
a “matching” is recorded. The probability of 4 or more matchings is tabulated to 
5D for s = 2(1)11 for c = 1, 2; for s = 2(1)8 for c = 3; for s = 2(1)5 forc = 
4,5; for s = 2, 3 for c = 6, 7; and for s = 2 for c = 8(1)12; in each case for all 
possible values of h. The fitting of various standard statistical distributions as an 
approximation to other cases is discussed. It is concluded that the normal distribu- 
tion is frequently wide of the mark, and that the Gram-Charlier type B is the most 
generally accurate. 


T. N. E. GREVILLE 
Social Security Administration 
Washington, D. 
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18[K ].—A. Ciirrorp CoHEN, Jr. ‘Maximum likelihood estimation of the dispersion 
parameter of a chi-distributed radial error from truncated and censored samples 
with applications to target analysis,” Amer. Stat. Assn., Jn., v. 50, 1955, p. 
1122-1135. 


The radial error (from the true mean) in a p-variate normal distribution with 
covariance matrix oI, is distributed, essentially, as x? with p degrees of freedom. 
The maximum likelihood estimate of ¢ is to be computed on the basis of a truncated 
sample of n observations in the range 0 < r < 1, e.g., as in a circular target. Set- 
ting $ = 1ro/o, the maximum likelihood estimate is obtained by solving 


P 1} p_ 9 (&) | 
2 re _S- = = 
E rejart = 2 BP - oe 
for & and converting to the equivalent ¢ where #,(f) is the standard x* density 
at the point & for p d.f. and J,(&) is the upper tail area of the same distribution. 
Table I gives G,(£) for p = 2and p = 3 to 5D for — = .0(.1)4(.5)5. 


Leo Katz 


Michigan State University 
East Lansing, Michigan 


19[K ].—A. Hurtson, “A method of assigning confidence limits to linear combina- 
tions of variances,” Biometrika, v. 42, 1955, p. 471-479. 


Consider estimation of a linear function of r unknown population variances 
01, ‘** , o,*. The estimation data consists of the corresponding statistics, s;*, --- , s,?. 
The s? are mutually independent and s?/o? has a x? distribution with f; degrees of 
freedom (i = 1, --- , r). The problem is to obtain approximate confidence limits 
for >-* \w#, where the constants \,; are arbitrary but specified. These confidence 
limits are obtained by determining an approximate expression for the statistic 
y = y(s?, ---, 5, P) with the property that Pr[ >> \us2/ > dw? < y] = P. 
Let denote the 100P percentage point of the normal distribution with zero mean 
and unit variance and Vin = (>, As™s2"fi)/ (> As2)™. Then, to terms of 
order 1/f?, y = ho + hi + he + hs. Here Mo = 1,h = EV2V2, 


he = 24/2Va, — 4Ve3/3Vn] — 2Vn/3Vnn, 
hs = V2Valtl— 4V02(3Ver)7 + 9Vesv2 Vn)? — Veo( Vax) 
— 23(Vs2/3)*Vu)*] + S[— 8Vn(3V2)* 
+ 2Va — 16(Vs2/3)*( Van) + 5Ves(v2 Var)~. 


For r = 2, the \,; of the same sign, and f; > 16, two decimal accuracy appears to 
be obtainable from the relation y = ho + Ii, + he. Tables 1-4 contain 2D values 
of y for r = 2; (Aisi? + Aose*) = 0(.1)1.; fi, fo = 16, 36, 144, o; P = .01, .05, 


95, .99; and Ade > 0. 


J. E. WatsH 


Lockheed Aircraft Corporation 
Burbank, California 
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20[(K ].—C. RADHAKRISHNA Rao & I. M. CHAKRAVARTI, “Some small sample tests 
of significance for a Poisson distribution,” Biometrics, v. 12, 1956, p. 264-282. 


The authors investigate 3 problems associated with the Poisson distribution 
particularly if the parameter y» is small: (1) homogeneity of observations, (2) 
goodness of fit, and (3) deviation in the frequency of any class from its theoretical 
value, usually the zero class. For testing of homogeneity the authors compare the 
classical x? test with a newly derived likelihood ratio test. They conclude that the 
x? test is adequate for 7/f > 3 where T is the total of the f observations. Four 
tables are provided for exact tests and in all tables T = 3(1)10, f = 3(1)10(10) 100. 
In the first table the critical values of }>x,?, x; being the value of the i-th observa- 
tion, are given to 1 or 2 S such that the significance level is less than or equal to 
.05. The actual significance level is tabled. Similarly, the critical values of }>x; 
log.x;, the likelihood test for homogeneity, are given to 2 or 3 S. The third table 
lists the critical values of }- f-log.(f,r!), f- = observed frequency of the variate 
with value r. This is the likelihood ratio test for the observed distribution to be 
Poisson. In the fourth table, the authors give the critical value of the frequency 
of the zero class in the Poisson. The authors also investigate truncated Poisson 
distributions and illustrate the theory with interesting examples. 


L. A. AROIAN 


Systems Development Laboratories 
ughes Aircraft Company 
Culver City, California 


21(K].—HeEnry J. Horn, “Simplified LDs0 (or EDso) calculations,” Biometrics, 
v. 12, 1956, p. 311-322. 


Let F(x) bea cumulative distribution function and x; < - - + < x,,and let p; = r;/n; 
be a binomial estimate of F(x;) based on n; independent trials. Thompson [1] has 
given a moving average method for estimating the value of m such that F(m) = 3, 
ie., the LDg or EDso, in dose-response assays. In this method only the k values of 
P; most nearly bracketing 4 are used. Weil [2] gives tables to facilitate the calcu- 
lation of m and approximate confidence limits as a function of ras1, -** , Taye in 
the special case nm; = m and log x.41 = log x, + d. If k = 4, = 4 or 5, and the 
values of x are in geometric progression with x = 10°” or x = 10°” (s an integer), 
the present tables give directly the estimate of m and approximate confidence 
limits, both to 3.S, for the more probable of the n* possible outcomes, (fa41, --* , 
Tak) . 


PauL MEIER 


The Johns Hopkins University 
Baltimore, Maryland 


1. W. R. THomrson, ‘‘Use of moving averages and interpolation to estimate median effective 
dose,” Bact. Rev., v. 11, 1947, p. 115-145. 

2. Carrot S. WEIL, “lables for convenient calculation of median-effective dose (LDso or 
EDso) and instructions in their use,’’ Biometrics, v. 8, 1952, p. 249-263. 
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22[K ].—Gortrriep E. NoeTHER, “Two sequential tests against trend,” Amer. 
Stat. Assn., Jn, v. 51, 1956, p. 440-450. 


Table I aids in applying the Wald sequential probability ratio test to a binomial 
population when Hy: p = 4 and Hi: p = 3 + «. The table is for a = 8 = .05 and 
e = .01(.01).49. Values of the slopes, intercepts, and average sample numbers are 
given. 

Table II contains values of 


Ce dem 
(x 2x Jo . 


for x = .01(.01) 1.00(.02) 2.00(.04) 3.00 to 5D. 
I. R. SAVAGE 


Center for Advanced Study 
in the Behavioral Sciences 
Stanford, California 


23[K ].—Dovuc tas G. CHAPMAN, “Population estimation based on change of compo- 
sition caused by a selective removal,” Biometrika, v. 42, 1955, p. 279-290. 


A population consists of two classes with sizes X; and Y; (X; + Y; = N;,) at 
times ¢; (¢ = 0, 1). Samples of mo and m are taken from this population at t and 
i, respectively. Between these two periods, R, = X; — Xo are removed from the 
X-class and R, = ¥; — Yo from the Y-class. Formulas are given to estimate No 
and X» and the variances of the estimates. Table 1 presents values to 2D of the 
optimum sample ratio m/mo to estimate No when R, = 0 for selected values of 
(R./Xo): .01, .02, .05, and .1(.1).5. Table 2 presents simliar values to 3D to esti- 
mate Xo for selected values of (Xo/No): .1(.1).9; and (R./No): .01, .02, and 
.05(.05).25. 

The author also compares the cost of the above method relative to the familiar 
tagging or capture-recapture procedure. Table 3 presents minimum values to 2D 
of the relative cost of tagging to the cost of classification such that the variance of 
No is greater for the tagging procedure; this table is given for the selected values 
(Xo/No) and (R./No) used in Table 2. 

R. L. ANDERSON 


North Carolina State College 
Raleigh, North Carolina 


24[K ].—RotF BarRGMANN, “Signifikanzuntersuchungen der einfachen struktur in 
der Faktoren-analyse,” Mitteilungsblatt fiir Math. Stat., v. 7, 1955, p. 1-24. 
Let P, = I 0((k = 1) /2, 3), 

P(k,r,n) = >(" 1-9 i) PH, — P= Ip(r —k + 1,2 —1 +1), 


where I,(a, b) is the incomplete Beta function. Table I gives the values of P; to 
4D for k = 2(1)16. (Ps; = .1495 and not the .1459 entered in the table, Py = .2549 
instead of .2548). The principal table gives integral values of r approximately 
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satisfying P(k, r, n) = « for a = .001, .01, .05, .10, .25, 50; & = 2(1)12. The 
tabled r is usually the largest integer satisfying P(k, r, m) < a. In some cases the 
author uses the next lower value of r when the resulting P(k, r, 2) — a@ is small; 
in these cases the value of r is printed in italics. The rules followed by the author 
in deciding whether P(k, r, m) — a@ is small are not clear. 

INGRAM OLKIN 


Michigan State University 
East Lansing, Michigan 


25[K ].—S. B. Cuaupuurt, “Statistical tables and certain recurrence relations 
connected with p-statistics,” Calcutta Stat. Assn., Bull., v. 6, 1956, p. 181-188. 


In samples of m, and m, from normal p-variate universes, in order to test the 
hypothesis that the two universes have the same set of variances and covariances 
the fundamental test statistics are the p roots, K;*, K:?, --- , K,’, called p-statistics, 
of the determinantal equation | A’;; — K?A";;| = 0 in which (A’;;) and (A’’%;) 
are the variance-covariance matrices of the two samples. With A = (m — 1)/ 

\K? 
(m= 4) 8 eT 
the 5% and 1% points of # for = 2 and m = n = 2.5(.5)5(1)11 where m = 
(m — p—2)/2 and n= (m— p— 2)/2. This extends the table of Nanda [1]. In 
addition for p = 3, 5% and 1% points are here given for # to 2D for m, m = 0(.5)2 
and m = n = 2.5(.5) (1)8. 


(the maximum # so defined), the authors gives to 2D 


... Cake 


1. D. N. Nanna, ‘Probability distribution tables of the larger root of a determinantal equation 
= pi all Indian Soc. Agr. Stat., Jn., v. 3, 1951, p. 175-177. [RMT 1082, MTAC, v. 7, 
1953, p. 95. 


26[K ].—N. MARAKATHAVALLI, “Unbiased test for a specified value of the para- 
meter in the non-central F distribution,” Sankhyd, v. 15, 1955, p. 321-330. 


Given a sample of m observations {x,;} from normal populations with common 
variance o* but unequal means {a,;}. Then )~*_, x2/o? is distributed as a non- 
central x’? with n degrees of freedom and non-centrality parameter,’ = )~"_, a2/os. 
The author considers the non-central F’-distribution, which is the ratio of a non- 
central x:’* and a central x? with respective degrees of freedom, » and v2, under the 
null hypothesis, Hy: 4 = do. A table is presented to determine the rejection region 
for an unbiased test of Ho against the two-sided alternative \ + Xo. The region is: 
F’ < a or > a, where Prob (F’ < a;| Ho) = a and Prob (F’ > a| Ho) = a, 
a + a = a = .05. This table presents values of a; and a, to 3D, a; generally to 
2D and a, to 1D, for » = 1(1)8, » = 4(4)20, 40, 60, © and 


= Xo = 
7) Vz +1 ae 


The author also indicates how this table can be used for analysis of variance, 
correlation ratio, multiple correlation, and discriminant analysis problems. 
R. L. ANDERSON 


North Carolina State College 
Raleigh, North Carolina 
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27[K ].—_K. J. RAMACHANDRAN, “On the simultaneous analysis of variance test,” 
Ann. Math. Stat., v. 27, 1956, p. 521-528. 


The author first considers the general problem of determining the a; in the 
equation 


>| < ad¢;/m;i = 1,2, --- | =1-—a, 
where the S; and S are mutually independent; S;/o* is distributed as x’ with ¢; 
d.f. and S*/o* as x? with m d.f. A method of evaluating P is presented and then a 
special method is developed for the case ¢; = #, for all ¢. In the latter case P becomes 


Pl = = < |, 


where S,,.: is the maximum of the S; and b = at/m (a; = a, for all 4). Upper 5% 
points of u are given to 2D for k = 2;¢ = 1(1)4(2)12, 16, 20; and 


m = 5(1)8(2)12(4)24, o. 


R. L. ANDERSON 


North Carolina State College 
Raleigh, N. C. 


28[K ].—Martin Fox, “Charts of the power of the F-test,” Ann. Math. Stat., v. 
27, 1956, p. 484-497. 


Write f; and fe for the degrees of freedom of numerator and denominator, re- 
spectively, of F. Let a be the size and £8 the power of F. 
Write 
@ = (S.4/(fi + 1)0*}}, 


where S,* is the value of S,? when the observable random variables are’ replaced 
by their expectations under the alternative hypothesis, and S,? is the sum of squares 
in the numerator of F. Eight charts show constant ¢ for all combinations of 8 = 
.5,.7, .8, .9, and a = .01, .05. The charts are designed to show readily the combina- 
tion of f; and fz required to obtain power 8 against a given alternative. A detailed 
numerical example, two nomograms facilitating interpolation in 8 for both values 
of a, and references to earlier tables are included. 
H. A. FREEMAN 


Massachusetts Institute of Technology 
Cambridge, Mass. 


29[(K ].—H. A. Davin, “On the application to statistics of an elementary theorem 
in probability,” Biometrika, v. 43, 1956, p. 85-91. 


The author is concerned with statistical applications of the following theorem 
in probability. Let Pi;... denote the joint probability of r(< m) events A, 


Aj, --+ , Ax; and S, the sum of the (") P’s with r different subscripts; then the 
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probability ,,,, of the realization of at least m events out of n is given by 


= etiohaeeie 
Pun = 2, ( 1) m—1 )o= 
This result is applied to the distribution of the extreme deviate from the sample 
mean, and the maximum F-ratio s2,,,/s*. The author's table 1 gives the value of 
P(A), the probability to 3D of rejecting the largest of m observations by the use of 
Xmax — £ at the 5 per cent level of significance, when all the observations are normal 
with unit variance, m — 1 have mean yu and one has mean yw + X, for A = 1(1)4, 
m = 3(1)10, 12, 15, 20, 25. 

In table 2 are given the probabilities to 4D for various designs of size / X 1, 
(randomized blocks, Latin squares, and Graeco-Latin squares) with which the 
ordered F-ratio F(,) exceeds the upper percentage point F, of the corresponding 
random F-ratio where / = 3,¢ = 1, 2, 3;1 = 5, 7,9, ©;¢ = 1(1)4, a = .05, and 
0.01. 

In table 3 are given the upper 100a percentage points, to 3S usually, of the 
maximum F-ratio followed by the corresponding 100a/n percentage points of F 
for the m-factor designs of table 2 where / = 3, 5, 6(1)9, m = 2, 3, 4, a = 0.05 
and 0.01. 

L. A. AROIAN 
Systems Development Laboratories 


ughes Aircraft Company 
Culver City, California 


30[K ].—K. C. S. Pixar, “On the distribution of the largest or the smallest root 
of a matrix in multivariate analysis,” Biometrika, v.43, 1956, p. 122-127. 


The joint distribution of the roots of a certain matrix is well known, 


Ps, +++ 4 84) = C(s, m,n) TT 9n(1 — 0° TT (0s — 09, 


0O<&A<S°-' 58 < 1. 


The author gives approximations for the upper percentage points for the largest 
(smallest) root for s = 2, 3, 4, 5. For s = 2, Nanda [1] tabulated the inverse 
function p (.95), P-! (.99) for m = 0(4)2, m =4(3)10, where P(x) = P(# < x). 
The author tabulates P“ (.95), P- (.99) to 3S for m = 0(1)4, m = 5(5)30, 
40(20) 100, 130, 160, 200, 300, 500, 1000 using the approximation 


1 — P(x) = x"™(1 — x)"4[2B(m + 1, + 1)B(2m + 2, 2n + 2). 


INGRAM OLKIN 
Biahigns State University 
East Lansing, Michigan 
1. D. N. Nanna, “Probability distribution tables of the larger root of a determinantal equa- 


tion with two roots’, Indian Soc. Agr. Stat., Jn., v. 3, 1951, p. 175-177. [RMT 1082, MTAC, 
v. 7, 1953, p. 95.] 





REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 77 


31{K ].—D. E. Barton & F. N. Davm, “Tests for randomness of points on a line,” 
Biometrika, v. 43, 1956, p. 104-112. 


The present paper is one of a set of papers by the same authors dealing with 
ordered intervals [1]. On a line of unit length, (m — 1) points are dropped randomly 
forming m intervals. Let the distances from one end of the line be x; (4 = 1, 2, 

-,n — 1) witha < x2 < +++ < Xp. Letd = nm, 


Let {g;} be {d,;} with g, the smallest d;, g. the next larger, etc. Let 7; be the solu- 


tion of d,, = g;. Let r;’ be the solution of d; = g,,,. The paper deals with distri- 
butions of (7) 


R= irl, (ii) Y = D ids Gil) G = dre 


t n 
(iii’) Gan = Sorgs, and (iv) Gt=n Dida. 
i=—s t—1 


Rand Y are disposed of quickly by references to papers by other authors and by 
references to the use of the normal approximations to R and to m — Y for moderate 
or large . Moments of G,, are discussed. Four moments of G are displayed. Table 3 
gives y:’/n to 4D, w,/n to SD, nf; to 3D and n(6, — 3) to 3D for m = 5(5)25, «; 
p = s/n = 1/5((1/5)4/S. uw’ is the 1st moment, 6; = js*/ps*, Be = ps/us*, ue, 
ps, wa are central moments. For G*, u:', ue, 8; and 6 are displayed as functions of n 
and values of 6; and 6 are given to 2D in Table 5 for » = 2(1)10, 20, 40, 100, 
400, 1000, «. Percentage points of G* to 3D and of standardized G* to 2D are 
given in Table 6 for » = 40, 100, 400, 1000, and for percentiles 0.5, 1, 2.5, 5, 95, 
97.5, 99, 99.5. 


K. J. ARNOLD 


Michigan State University 
East Lansing, Michigan 


1. D. E. Barton & F. N. Davin, ‘“‘Some notes on ordered random intervals,” Roy. Stat. Soc., 
Jn., ser. B, v. 18, 1956, p. 79-94. (MTAC, this issue, Rev. 36, p. 80.) 


32[K].—W. D. Ray, “Sequential analysis applied to certain experimental designs 
in the analysis of variance,” Biometrika, v. 43, 1956, p. 388-403. 


Given the one-way classification model 
Xu=a + by + ac(t = i, “igpies , R34 = .. the »), 


where the z,; are uncorrelated with zero means and equal variances, o?, the ; are 
treatment parameters (>. 5, = 0). We consider the non-centrality parameter 
8’ = >> b,/ko* and the testing statistic 


ee Treatment Sum of Squares 
Error Sum of Squares 
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Two hypotheses are considered: Hy(é = 0) and H; (4 some specified value ~ 0). 
A sequential sampling procedure is set up (consider successive values of ) so that 
after m samples per treatment one of these 3 decisions is made: 


Result of n samples per treatment Decision 
Me? > (1 — B)/a accept Hi 
Mer? < B/(1 — a) accept Ho 
B/(1—a) < Me < (1-—8)/e take more samples 


M is a confluent hypergeometric function, 


u|* tena + a)|; 
\ = kné; and a and £ are the usual Type I and Type II error rates. 

Let M(G)e™? = (1 — 8)/a and M(G)e™? = 8/(1 — a). Values of G and G 
to 3D have been computed for a = 6 = .05; selected values of n; for 6 = 0.5, 
k = 2(1)7; for 6 = 1.0, & = 2(1)10; for 6 = 2.0, k = 2(1)6. Certain of the values 
of G and G were computed by approximation methods. Tables 7 and 8 present some 
comparisons of exact and approximate results. 

The results are extended to a two-way classification model with k blocks and n 
treatments (note the unfortunate interchange of m and k from the one-way model). 
In this case kn — 1 is changed to k(m — 1) and k — 1 tom — 1 in M. Tables of G 
and G generally to 3D are presented as before, except that for 6 = 0.5 and 1 0, 
k = 2(1)8. Conjectural expected sample sizes are presented to go with these tables, 
based on a generalization of Wald’s formula. 


R. L. ANDERSON 


North Carolina State College 
Raleigh, N. C. 


33[K].—S. R. BRoapBENT, “Lognormal approximation to products and quo- 
tients,” Biometrika, v. 43, 1956, p. 404-417. 


Table 1 of this paper gives upper and lower 1% and 5% points of a variable 
whose logarithm obeys the normal distribution law in terms of its coefficient of 
variation, v, to 4D for v = 0(.0001).15. The actual percentage point is obtained 
by multiplying the table entry by the mean value. First differences are given. 


AB ede 4 


34[K].—D. J. BARTHOLOMEW, ‘‘A sequential test of randomness for events occur- 
ring in time or space,’ Biometrika, v. 43, 1956, p. 64-78. 


The paper is concerned with a sequential test of a = 0 against a = a» + 0 for 
the class of processes for which Pr {event in (7, T + dT)} = w(uT)dT + 0(dT), 
(a > —1). At the cost of one observation, a transformation eliminates the nuisance 
parameter yu. The Wald formulas for the operating characteristic and expected 
sample numbers are then applied. Table 3 gives E.(| ao) for various values of 
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a and a, and for (a, 8) = 0.01, 0.05 where 


— L(a) log B + [1 — L(a)] log A 
. Pel | a0) = Tog (+ a) = ae +0)’ 


A = (1 — 8)/a, B = B/(1 — a), L(a) = (A*™ — 1)/(A*™ — B®), h(a) is 
the non-zero solution of 1 + hao/(1 + a) = (1 + ae)*. The value a’ of a for 
which (1) is indeterminate is near that which maximizes E,(n | a)-E.-("| ao) = 
— log A log B/[log (1 + ao) }. Table 4 gives E,(m | ao) and E,,(n | ao) for various 
values of ao and for (a, 8) = 0.01, 0.025, 0.05. Integral values are given in each 
case. Table 5 lists percentage savings in sample size in using the sequential test 
rather than a corresponding fixed sample size test. 





K. J. ARNOLD 


Michigan State University 
East Lansing, Michigan 


35[K ].—H. A. Davin, “A note on moving ranges,” Biometrika, v. 42, 1955, p. 
512-515. 


This publication deals with some of the statistical properties of moving ranges, 
moving maxima, and moving minima which are of course based on overlapping data. 
In particular, an expression is obtained for the correlation between two ranges in 
overlapping samples in terms of means, variances, and covariances of order sta- 
tistics. With the derived result and the hypothesis of a normal population, the 
efficiency of mean moving ranges is studied. 

Letting x», x@), --- , Xun represent the observations ordered in time, we may 
say that x9, x(i4), -** » X(s+n—_» is the ¢-th moving sample of size m and designate 
m;(n), m,/(n), and w;(m) as the corresponding maximum, minimum, and sample 
or group moving range respectively. Writing the order statistics of a sample of 2 
as x; > x, > x3 *** > Xn, values are given to 4D by David for the correlation 
between #() and #144(n) for m = 2(1)5 and d < (m — 1) and for the correlation 
between two ranges w, and w, + d in overlapping samples in terms of known mo- 
ments and cross-product expected values of the x;. Now the first mean moving 
range is defined as 

oil aes @; 


i=1 N-—-n+1 


and a mean moving range w, can be constructed on every r-th sample. David then 
studies the efficiency of w,, i.e., 


Eff = var S’/var w, 
where S’ is the usual (total sample) estimate of population standard deviation 


,_ T(N=1/7) 


which is tabulated to 3D for N = 10(10)60, 50, 100, 200, ©, and for m = 2(1)5 
for w; for nm = 4 for w; for n = 5 for a, the ordinary mean range estimation ; and 
w* the most efficient (non-overlapping) weighted mean range estimator [1 ]. 
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In the absence of trends, w:(5) is found to be considerably more efficient than 
w,(2), the ordinary mean successive difference. The statistic w:(2) is less influenced 
by trends in observed data, however, than the statistics w,(m) based on larger n, 
Compared with @ for groups of 5 observations each and w*, the moving range esti- 
mators, w, are seen to gain in efficiency as the total sample size N increases. The 
relatively low efficiency of the mean moving range for small total sample size NV 
and small group size m is due no doubt to the lack of weighting. 


F. E. Grusss 
Ballistic Research Laboratories 
Aberdeen, Maryland 


1. Frank E. Grusps & CHALMERS L. WEAVER, ‘‘The best unbiased estimate of population 
standard deviation based on group ranges,” Amer. Stat. Assn., Jn., v. 42, 1947, p. 224-241. 


36[K ].—_D. E. Barton & F. H. Davin, “Some notes on ordered random intervals,” 
Roy. Stat. Soc., Jn., ser. B, v. 18, 1956, p. 79-94. 


In this paper, which is a further investigation of the distribution of the ordered 
intervals obtained by dropping m — 1 points at random on a line segment of unit 
length, some short tables are included. First the values of the skewness (third 
standard moment) and the kurtosis (fourth standard moment minus three) are 
given to 4D for the distribution of the r-th greatest interval for m = 5(5)25 and 
r = 1(1)n. Next the coefficient of correlation is given between the two shortest 
intervals and between the first and third shortest intervals to 3D for nm = 2(1)10, 
15, 25, 50, 100, «©. Then the limiting values of +/n times the skewness, of n times the 
kurtosis, and of m? times the variance divided by the expected value are given to 
4 and 5S for the deciles. Turning to the distribution of ratios of two intervals, the 
1%, 5%, 10%, 90%, 95%, and 99% points are given for the ratio of the smallest 
to the largest interval to 3S for nm = 3(1)10, 15, 20. Finally upper 5%, 2.5%, 1%, 
and 0.5% points for the interquartile differences of interval length are given 4 or 
5D for n = 7, 11, 15, 19, 23. 

ORT ay oe 


37[K ].—A. E. SarHAN, “Estimation of the mean and standard deviation by order 
statistics, Part III,” Ann. Math. Stat., v. 26, 1955, p. 576-592. 


The data considered consists of singly and doubly censored samples from speci- 
fied types of probability density functions. The best linear estimates of the popula- 
tion mean and standard deviation using sample order statistics, also their variances 
and their efficiencies (compared to best estimates using uncensored samples), are 
derived for rectangular and exponential distributions and arbitrary sample size n. 
The same results are given for sample sizes up to 5 for other types of distributions. 
The symmetric distributions considered are U-shaped, rectangular, parabolic, tri- 
angular, normal, and double (tailed) exponential. Let the smallest r; and the 
largest rz sample elements be missing. Table I contains coefficients of the best 
linear estimates of the mean and standard deviation for these symmetric distribu- 
tions in the singly censored case. Table II contains the variances and efficiencies 
for the estimates of Table I. For Tables I and II, 1 = 0, 72 = 1(1)” — 2,” = 
3(1)5. Tables III and IV contain results of the same nature as Tables I and II, 
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respectively, for the case of doubly censored samples. In Tables III and IV, 1 = 
1, ro = 1(1)m — 1, and m = 3, 4. Results are also computed for a third degree 
polynomial skewed distribution. Table V contains coefficients of the best linear 
estimates of the mean and standard deviation of this skewed distribution for 
singly and doubly censored samples. Table VI contains the variances and effi- 
ciencies of the estimates of Table V. Table VII contains a comparison of the vari- 
ances of the estimates for the skewed distribution with the variances of the corre- 
sponding estimates for the exponential distribution. For Tables V, VI, and VII, 
n = 3(1)5 and various combinations of 1:1, 72m — 2 are considered. With the excep- 
tion of efficiencies, the results in all Tables are expressed to 7D. Efficiencies are 
expressed as percentages and to 2D. 
J. E. Watse 


Lockheed Aircraft Corporation 
Burbank, California 


38[K].—V. N. Murty, “The distribution of the quotient of maximum values in 
samples from a rectangular distribution,”” Amer. Stat. Assn., Jn., v. 50, 1955, 
p. 1136-1141. 


These tables are designed to test the hypothesis that two samples have been 
drawn from the same rectangular population with a known lower bound. Let Li 
be the maximum in a sample of size m, and L, the maximum in a sample of size n. 
The distribution of u = L/L, is obtained, and is independent of the parameter. 
The power of the test is calculated and 5 percent points of u for samples of size 
m,n = 2(1)10 to 3D is computed. 

J. R. VATNsDAL 


State College of Washington 
Pullman, Washington 


39[K ].—ALBERT Ross EcKLER, “Rotation sampling,” Ann. Math. Stat., v. 26, 
1955, p. 664-685. 


This paper shows how to find minimum-variance estimates of the mean a(é,) 
of a time-dependent population, assuming that one is restricted to the class of 
linear unbiased estimates. Each minimum-variance estimate is based on a specified 
sample pattern (a set of sample values drawn from the population at one or more 
distinct times). Let the random variable X;; denote the value of element j of the 
population at time ¢;. The correlation between X;; and Xq is assumed to be p! 4; 
the correlation between X;; and X % is assumed to be zero; the variance of X,; is 
assumed to be o? independent of time. Iterative methods are developed; the esti- 
mate of the population mean a(t;-1) is used in determining the population mean 
a(t;). 

The paper discusses two important methods of sampling: in one-level rotation 
sampling, the statistician can add to the sample pattern only sample values that 
have been drawn from the population at the current time; in two-level (and higher- 
level) rotation sampling, the statistician can add earlier sample values as well as 
current ones to the pattern. Schematic sample patterns associated with these two 
methods are illustrated. 
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Formula 3.2 gives an iterative scheme for finding the minimum variance esti- 
‘mate of a(t;). This formula involves an unknown coefficient A;. Table 3.1 gives 
the values of A; to 4D as a function of p for p = 0(.1)1, .75, .85, .95 for the case 
where half of the elements in the sample at time 4—; are retained in the sample 
drawn at time ¢, and each of the remaining elements is replaced by a new element. 
Formula 4.2 gives an iterative scheme for estimating a(¢;). This formula involves 
an unknown constant a;. Table 4.1 gives values of a; for 7 = 1(1)11 as a function 
of p for the same values of p as before. 


BENJAMIN EPSTEIN 
Wayne State University 
Detroit, Michigan 


40[K].—P. ArmsEn, “Tables for significance tests of 2 X 2 contingency tables,” 
Biometrika, v. 42, 1955, p. 494-511. 


In 2 X 2 tables it is assumed the frequencies result from a double dichotomy of 
N individuals according to the properties X and Y. The null hypothesis to be 
tested is: There is no association between X and Y. As is well known, the whole 
table—for fixed marginal totals—is determined by one of the entries and the 
probability of any set of frequencies is given by the exact hypergeometric formula. 
The object of this paper is to provide tables of 5 and 1 per cent significance levels 
for both one- and two-tailed tests containing up to N = 50 observations. 

Two possible procedures are given, called D,: Arrange the possible events in 
ascending order of the size of their probabilities under the null hypothesis, include 
in the 100 @ per cent two-tailed rejection region whose events for which the cumu- 
lative sum of these ordered probabilities is smaller than or equal to a. D;: Define 
F(E) or the “‘first-tail’”’ probability as the cumulative sum of the probabilities under 
the null hypothesis of all possible events which are more extreme, in the same 
direction, than a given event E, including E itself. Define S(Z) or the ‘“‘second-tail” 
probability as the cumulative sum of probabilities, starting with that of the event 
most extreme in the opposite direction as compared with Z, and cumulating up to 
but not exceeding the value of F(Z). If, and only if, F(Z) + S(Z) < a, include E 
in the rejection region for the two-tailed 100 a per cent level of significance. 

Tables are computed for both these criteria, there being no difference in most 
cases. For every possible set of frequencies the tables give significance levels for 
one- and two-tailed tests for a = .01 and .05. Explicit directions are given for the 
use of the tables and examples worked. 

J. R. VATNsDAL 


State College of Washington 
Pullman, Washington 


41[K ].—A. R. Kamat, “A two-sample distribution-free test,”” Biometrika, v. 43, 
1956, pp. 377-387. 


The data consists of two independent samples of sizes m and n with m > n. 
These sample values are unequal. The null hypothesis asserts that these samples are 
from the same population. The alternative hypothesis of interest is that the samples 
are from populations with different dispersions. The samples are pooled and their 
values ranked. Let R, be the range of the ranks for the sample of m and R,, this 
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range for the sample of m. The test statistic is D,.. = R, — Rn — m. Table 1 
contains the largest integer values D, such that P(D,,_ < Da) < aand the smallest 
integer values D, such that P(D,,m > Da) < a, for 100a = .5,1,2.5,5andm +n = 
7(1)20 with » < (m + n)/2. Table 2 contains 3D values of the mean (4D), 
standard deviation, p2,+/ Bi, and 62 of Da, for m + = 16and n = 2(1)8. Table 3 
contains 3D values of the standard deviation, ps2, B2 of Da,» form = m = 2(1)5, 7, 
10, 15, ©. Table 4 contains moments and standardized deviates for the limiting 
distribution of d = Da.» — m for m/(m + n) = p = .5(.1).9. The contents of 
Table 4 are useful in approximating percentage points of D,,., when m + n > 20. 
Table 5 contains a comparison of the approximate and true 2.5 percent points when 
m+n = 20, p = .5(.1).9, m = 10(2)18. Table 6 contains corrected and stand- 
ardized upper and lower 2.5 percent points for d to 2D when p = .5(.05).75. - 


J. E. Watsx 
Military tions Research Division 
Lockheed Aircraft Corporation 
Burbank, California 


42[K ].—J. B. Douctas, “Fitting the Neyman type A (two parameter) contagious 
distribution,’’ Biometrics, v. 11, 1955, p. 149-173. 


While moment estimates of the parameters for a two-parameter Neyman Type 
A distribution are easily computed, the maximum likelihood estimates are obtained 
by solving a pair of equations known to require a very laborious iteration pro- 
cedure. The author rewrites these likelihood equations in a form which, with the 
aid of tables of auxiliary functions, makes it possible to obtain the solution by 
iterations rather quickly. The case of a two-parameter Neyman Type A distribu- 
tion truncated by omitting the zero class is also treated. Numerical examples are 
presented for the complete as well as the truncated distribution. The auxiliary 
functions p, and g, are defined by writing 


@ 
pe’ = €* Yor Ar/r! 
r=0 


and setting pz = u,41/u,; Ge = Pe(Per1 — Pz). Tables of p, and gq, are given to at 
least 3 and at most 5S for x = 0(1)19, X = .000(.001) .03(.01) .3(.1)3. 


Z. W. BiRNBAUM 
University of Washington 
Seattle, Washington 


43[K].—J. H. Bennett, “Partitions in more than one dimension,” Roy. Stat. 
Soc., Jn., ser. B, v. 18, 1956, p. 104-112. 


Though there is some discussion of partitions in more than two dimensions, the 
greater part of this paper is concerned with bipartitions, extending and clarifying 
the original work by Fisher [1]. The first table gives values of certain auxiliary 
functions and the second gives the numbers of bipartitions of m = 1(1)8 with given 
marginal partitions. 

Guid. C. 


1. R. A. Frsuer, ‘‘A class of enumerations of importance in genetics,” Roy. Soc., Proc., London, 
Ser. B, v. 136, 1950, p. 509-520. 
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44[K ].—Rosert Hooke, “Symmetric functions of a two-way array,” Ann. Math. 
Stat., v. 27, 1956, p. 55-79. 


Generalized symmetric functions, (g.s.m.’s) of the elements in an r X c array 
are defined as averages of restricted monomial symmetric functions; the restrictions 
apply to location of the elements in rows and columns. Thus, the g.s.m. of degree 
(a + 6 + c) with two elements in the same row and one in a distinct row and a 
distinct column is 


a b - 1 , b c 
bos gy ~ Faw Leh Xen 


where >.’ indicates summation over distinct subscripts. 

Bipolykays are defined by a non-commutative multiplication of the polykays, 
a family of polynomials in symmetric functions defined by Tukey [1] in generaliza- 
tion of Fisher’s cumulants. Therefore, the bipolykays are linear combinations, with 
integer coefficients, of the g.s.m.’s. 

Conversion formulas for the representation of bipolykays in terms of g.s.m.’s 
and the reverse are given in the text for functions of degree 1 and 2; Table 1 gives 
the coefficients for degree 3 and Table 2 for degree 4. Multiplication formulas for 
bipolykays of degree 1 by bipolykays of degree up to 3 are given in the text and, 
in Table 3, for products of two bipolykays both of degree 2. Table 4 gives the 
variances and covariances of the bipolykays 


Ft} GZ) (2 2) 
for matrix bisamples selected by taking from an R XC population matrix, inde- 
pendently, samples of r from the R rows and c from the C columns and forming the 
indicated submatrices. 
Leo Katz 


Michigan State University 
East Lansing, Michigan 


P e. Joan W. TuxeEy, “Some sampling simplified,” Amer. Stat. Assn., Jn., v. 45, 1950, p. 501- 
19. 


45[L].—A. D. Smirnov, Tabliisky Funkisit Eiri i Speisial’nykh Vyrozhdennykh 
Gipergeometricheskikh Funkisit dha Asimptoticheskikh Reshenit Differentsial’ 
nykh Uravnenit Vtorogo Portadka (Tables of Airy Functions and of Special 
Confluent Hypergeometric Functions for Asymptotic Solutions of Differential 
Equations of the Second Order). Izdatel’stvo Akad. Nauk SSSR (Press of the 
Academy of Sciences of the USSR), Moscow, 1955. 264 p. + two loose cards, 
26.6 cm. Price 28 rubles. 


This volume in the series of Mathematical Tables of the Computational Centre 
of the Academy of Sciences of the USSR contains tables of three types of functions. 
Apart from their occurrence in various physical investigations, the functions 
play an important part in work of A. A. Dorodnitsyn on asymptotic solutions for 





REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


large \ of second-order linear differential equations 
[k(x)y’] + D'r(x) + a(x) y = 0 


when r(x) and g(x) have singularities (zeros of poles) in an interval. The tables 
were in fact prepared on the initiative of Dorodnitsyn and in constant consultation 
with him. Some indications concerning this work on asymptotic solutions are given 
in the Introduction; for further details, reference is made to Dorodnitsyn’s mem- 


oir [1 ]. 


The functions first considered are solutions of Airy’s equation U’’(s) + sU(s) = 0. 
The particular solutions U;(s, 1), U2(s, 1) which are tabulated, along with their 
first derivatives U;'(s, 1), Us'(s, 1), are such that at s = 0 we have U; = U,’ = 1, 
U; = U;' = 0. Consequently 


Us(s, 1) = $5 (— 1)*1-4-7 «++ (3m — 2)s%/(3m)! 


= 3-187 (2) 52 J_1y5(3s*), 


= by (—1)"2-5-8 --- (3m — 1)s***"*/(3n + 1)! 
n=0 


i (3) s"* Jus ($s?) - 


The two functions and two derivatives are tabulated, with smoothed second 
differences, for s = 0(.01)10 to 5D on pages 27-46, and for —s = 0(.01)6 to 5D 
or 5s on pages 47-58. Also tabulated on pages 27-46, equally for s = 0(.01)10 and 
to 5D, but without differences, are three functions Vu(s), Vie(s), Vee(s) which are 
integrals of expressions involving the U functions (explicit formulas are given at 
the bottom of page 7). 

The functions tabulated invite comparison with the Ai(x), Bi(x) solutions of 
y” = xy which are employed in the British Association tables [2] and are now in 
common use. Formulas on page B17 of [2] show that Ai(x) and Bi(x) are linear 
combinations of two solutions, y:(x) and s(x), analogous to the U solutions; in 
fact, Ui(s, 1) = y:(—s), Us(s, 1) = —ye(—s). The Ai, Bi functions behave more 
neatly as x —» —o, whereas the y (or U) solutions behave more simply near 
x = 0 (or s = 0). Since cases do arise in which both boundary conditions relate 
to the origin, it is useful to have the U functions tabulated as well as the Ai, Bi 
functions. 

The second batch of functions are solutions Ui(s, a), Us(s, a) of a generalized 
Airy equation U’(s) + s*U(s) = 0, in which, if complications are to be avoided, 
a + 2 (supposed positive) must not be the reciprocal of an integer. The particular 
solutions denoted by Ui(s, a), Us(s, a) are again such that at s = 0 we have (for 
a> —1) U,; = U,’ = 1, U, = U1’ = 0. Power series are given, as also are expres- 
sions in terms of Bessel functions of order +1/(a + 2). On pages 61-156 the func- 
tions Ui(s, a), Us(s, a) with first differences, and also the derivatives U;'(s, a), 
U,'(s, a) without differences, are tabulated to 4D for s = 0(.01)6 and the following 
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16 values of a(6 pages for each a) : 
a= a3 4,3. 2) do oo BE 2, 
1 1 2 3 
~— ae ee —}, "He 7" e 


The third set of functions are solutions V:(s, po), V2(s, po) of the special confluent 
hypergeometric equation 


V"'(s) + [hbo(1 — 460)s* + s*]V(s) = 0 


in which complications arise when p is integral. The particular solutions denoted 
by Vi(s, po), Ve(s, oo) are normally the power series beginning with s'*. and 
s'12P, respectively; the appropriate modification when f) = 1 is indicated, since 
the tabulation extends to this value of fp. Although the fact does not appear to be 
stated, the solutions are simply expressible in terms of Bessel functions. When 
is not integral, we have 


Vi(s, Po) - T'(po) s'* J p,~1(2s"”), 
Vo(s, Po) = T'(2 — po)s!* Si_»,(25"”). 


When fp» = 1, the first formula remains valid, giving Vi(s, 1) = s'”Jo(2s'), while 
the second is replaced by V2(s, 1) = 42s! Y,(2s'”), where Y, (sometimes called No) 
is the Bessel function of the second kind in the form now usual. On pages 159-258 
the two functions Vi(s, po), Ve(s, #0) with first differences, and the first derivatives 
Vi'(s, po), Ve'(s, po) without differences, are tabulated to 4D for s = 0(.01)6, 
po = .1(.1)1. (In the sub-title of the table, on page 157, 0.1 < #) < 10 is a misprint 
for 0.1 < po < 1.0.) On page 259 are tabulated to 4D or SS the first five zeros, 
$1 to Ss, of both Vi(s, po) and Ve(s, po), again for Pp = .1(.1)1. 

On page 260 (and also on one loose card) coefficients $u(u + 1) of the second 
difference are tabulated to 4D without differences for —u = 0(.001)1. On page 
20 (and also on a larger scale on the other loose card) is a nomogram for evaluating 
$u(u + 1) de. On page 264 are given corrigenda to three other tables; see MTAC, 
v. 11, 1957, p. 125-126. 


A. F. 


1. A. A. DoropniTsyn, “Asymptotic laws of distribution of eigenvalues for certain forms of 


differential equations of the second order,” Uspekhi Matem. Nauk (N.S.), v. 7, 1952, p. 3-96. 
For Russian title see RMT 1113, MTAC, v. 7, 1953, p. 179. 

2. J. C. P. Miter, The Airy Integral, Part-volume B of British Assoc. Adv. Sc., Math. Tables, 
Cambridge University Press, 1946. 


46[L ].—NuMERICAL CoMPUTATION BuREAU, Report No. 8, ‘Tables of Whittaker 
Functions,” The Tsuneta Yano Memorial Society (1-9 Yuraku-cho, Chiyoda- 
Ku, Tokyo, Japan), 1954, 26 cm., 39 p. Price $1.00, and 


[L].—NumERIcAL CoMPUTATION BuREAU, Report No. 9, “Tables of Whittaker 
Functions (Wave Functions in Coulomb Field),”” The Tsuneta Yano Memorial 
Society (1-9 Yuraku-cho, Chiyoda-Ku, Tokyo, Japan), 1956, 66 p., 26 cm. 
Price $3.00. 
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These two reports contain tabulations of the solutions to the differential equa- 
tion 
ew k,i4-w\, _ 


known as Whittaker functions [1]. One form of the solutions, due to Whittaker, 
is denoted by Wi,m(z) and W;,-.(z) ; these are sometimes written in the form 


Wrm(2) = 7 *2tu, (Z) 


and the first derivative 


dWi..m(2) 
dz 


= e Miegky, (2). 


Another form of the solutions is commonly denoted by Mi,«(z) and M;,.(z); 
Mi,m(2) is related to the confluent hypergeometric function ;F;(a, 8, z) by 


Mi.m(z) = 2?*+™e | F,(4 + m — k, 2m + 1;2). 


The tabulation is for m = 3} and k in the range 0 to 1. Recurrence relations for 
generating solutions for other k and m are given in Report No. 8. The first derivative 
of the Whittaker functions is tabulated for use in the recurrence relations. 

Report No. 8 contains the following tables: 


(a) Mean(s), & Mian(z) with least significant figure 10-*, 2 = 0.0(0.1)5.0 and 


k = 0.00(0.05)1.00. Second central differences with z and & as variables 
are also given. 

(b) wzaa(z), tern(z) with least significant figure 10-*, 2 = 5.0(0.2)25.0 and 
k = 0.05(0.05)0.95; kuzi(2), kt 12(z) with least significant figure 10~*, 
1/z = 0.1(0.1)0.4 and k = 0.05(0.05)0.95. Second central differences with 
z and k as variables are also given. 

(c) Ex(a), E,(8) with least significant figure 10, a = 1 — 8 = 0.00(0.01)0.50, 
where E,(a) and £,(8) are the interpolation coefficients appearing in 
Everret’s interpolation formula. 

(d) Log s with least significant figure 10~, z = 1.00(0.01)2.00(1.00) 10.00, 100. 

(e) e~* with least significant figure 10, x = 0.00(0.01)1.00 and with four sig- 
nificant figures x = 1.00(1)10(10) 100. 

Report No. 9 is a table of Coulomb wave functions. It gives 


exp (— 5 )exo| — i(Jx — 1) | Wasa) 
where £ and x are real, 1 = m — 3, and G;,(x) and F;,(x) are defined as the real 
and imaginary parts of this function. As x — 0, Gz1,(x) — © so in the neighbor- 
hood of x = 0, a function defined by the authors as 

Pean(2) = Gean(2) + 5-$— log #Fran(«) 

2C*e.y2 

which remains finite at zero is tabulated rather than G;,(x). This formula for 
T':1a(x) appears inconsistent with other equations presented and with the tables 
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themselves. The correct equation for T';,:,2(x) appears to be 


—2F 


, g 
Teaa(x) = Gera(x) + (==) F:in(x) log x. 


For x > 1 the asymptotic form of Wi:12(—ix) is: 
Wian(—tx) = u(x) exp4i | 5 + + § log x + oo(&) + s(x) |} 


ase(—*) . 4 = 5, [Gean(*) + iFean(x)] 


pe dFrip(x) _ j@Gran(*) 
dx dx 


= 40(x) exp {i + & log x + oo() + x(a) |} 


where oo(€) = arg I'(1 — 2€). The functions u, ¢, v and x are tabulated for x > 5. 
Report No. 9 contains the following tables: 





(a) Pean(2), © Frean(x), Teara(x), £ Teaale) with least significant 10-*, x = 


0.00(0.05)1.00, € = 0.00(0.02)0.50(0.05)1.00. Second central differences 
with £ and x as variables are also given. 


(b) Frapn(x), 5 Fran(x), Gean(x), £Gean(z) with least significant figure 10-*, 


* = 1.00(0.05)2.00(0.10)5.00 and = 0.00(0.02)0.50(0.05)1.00. Second 

central differences with ¢ and x as variables are also given. 

(c) Urrn(x), dt1n(x), Vean(x), Xe.12(x) with least significant figure 10-*, 1/x = 
0.000(0.005)0.200 and ¢ = 0.1(0.1)1.0. Second central differences with 
¢ and x as variables are also given. 

These tables are among the most complete tabulations of these functions 
currently available. They are well laid out and easily read. Unfortunately, how- 
ever, they suffer from an inadequate explanation of how to use them There is in 
fact an error in the definition of I';,:2(x) in Report No. 9 as cited abo:'e. The nota- 
tion is also difficult to follow. Furthermore it should be noted that the functions 
F:in(x) and Gz1,(x) tabulated do not have the same normalization as the corre- 
sponding functions cited in the reference [2] given by the authors of the tables. An 
extensive discussion of the Coulomb wave functions with many formulae may be 
found in the National Bureau of Standards Tables [3]. No discussion of the method 
by which the tables were computed or of the estimated error is presented. 


Lioyp D. Fospick 
gital ep yo Laboratory 
Uelnarelt Illinois 
eer llinois 


a TAYLOR WHITTAKER & G. N. Watson, Modern Analysis, Cambridge University 
Pres, 194. 
bie dela & M. eres, “Excitation of atoms by electron collision,” Phys. Math. 
an "Proc., Japan, v. 22, 1940, 
3. NBS Applied Mathematics Soins oo "n Tables of Coulomb Wave Functions, v. 1, U. S. 
Gov. Printing Office, Washington, D. C., 
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47[L].—J. W. Heap & W. P. Witson, Laguerre Functions: Tables and Properties, 
Monograph. No. 183R, The Inst. of Elec. Engineers, London, 1956, 13 p. + 5 
tables + 4 p. of appendices, 27 cm. Copies obtainable from The Institution of 
Electrical Engineers, Savoy Place, London, W. C. 2. Price 1s.3d. 


This paper tabulates the Laguerre functions e* L,(x) for » = 0(1)20, x = 
0(.1)1(.2)3(.5)6(1)14(2)40(S)100. Most of the entries are given to four decimals. 
For large values of the argument two significant figures are given. In addition 
there are zeros for the same values of to six decimals. Comparison of these zeros 
with existing tables [1] for » = 1(1)16 serve to substantiate the author’s claim 
that the values are uncertain in the fifth decimal place. However, the discrepancy 
is never more than one unit. Accompanying the tables is an attractive computation 
of the most important mathematical properties which should certainly serve the 
needs of engineers. An omission noted is Tricomi’s expansion [2] for L,(x) in 
terms of Bessel functions J,(2v2x). 


M. ABRAMOWITZ 
National Bureau of Standards 
Washington, D. 


1. Herpert E. Satzer & RutH Zucker, “Tables of the zeros and weight factors of the first 
fifteen Laguerre polynomials,”” NBS Applied Mathematics Series No. 37, Tables of Functions and 
Zerosof Functions, 1954, p. 191-199. [Amer. Math. Soc., Bull.,v. 55, 1949, p. 1004-1012.] [MTAC, 
Rev. 104, v. 10, 1956, p. 249-250. 

2. Harry BATEMAN, Higher Transcendental Functions, v. 2 (A. Erdélyi, Director), McGraw- 
Hill Book Co., Inc., New York, 1953, p. 199. [MTAC, Rev. 67, v. 11, 1957, p. 114-116.] 


48[M ].—E. J. Scort, Transform Calculus with an Introduction to Complex Variables, 


Harper and Brothers, New York, 1955, viii + 330 p., 23 cm. Price $7.50. 


The present book is intended for those who wish to learn the application of 
transform theory to the solution of problems in engineering, physics, and chemistry. 
It is devoted mostly to the theory and application of the Laplace transform, al- 
though in the final chapter the finite Fourier and Hankel transforms are introduced. 
The bulk of the book consists of the application of the Laplace transform to a great 
variety of problems in diffusion, vibrations of rods, deflection of beams, electrical 
transmission, etc. The early portion contains a concise but readable account of the 
necessary mathematical background. 

Chapter I is a summary of the elementary portion of the theory of analytic 
functions centering around Cauchy’s integral theorem: Taylor and Laurent series, 
singularities, residues, and contour integration. In Chapter II the Laplace trans- 
form is defined, its existence established under suitable conditions, and the complex 
inversion formula obtained. This latter takes for granted and as point of departure 
the Fourier integral theorem. Chapter III is devoted to formal properties of the 
transform-translations, convolutions, etc. The following six chapters are devoted to 
a wealth of applications: ordinary linear differential equations and systems of such, 
difference equations, partial differential equations, special functions, asymptotic 
series, and integral equations. The final Chapter X defines and applies the finite 
Fourier and Hankel transforms. Tables of the several transforms are appended. 

Some, but not many, misprints were noted. The integrand of the last integral 
of equation 8 -2 -6 should have ¢* in its denominator instead of 2. 

Particularly because of the abundance of well worked out examples, the engi- 
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neer or physical scientist should find this a very useful book, and it should provide 
him with a broad background of problems to which the transform method is ap- 
plicable. 

The tables appended to the book consist of 

(a) A table of operations with Laplace Transforms (9 entries) 

(b) A table of Laplace Transforms (111 entries) 

(c) Tables of finite sine transforms (14 entries) and cosine transforms (11 entries) 


Joun W. GREEN 
University of California 
Los Angeles, California 


49[P].—T. Y. Cuen, C. P. Stress, & N. M. Newmark, “Moments in simply sup- 
ported skew I-beam bridges,” University of Illinois Engineering Experiment 
Station Bulletin No. 439, Studies of Slab and Beam Highway Bridges, Part VI, 
1957, 72 p., 28 cm. Price $1.00. 


This investigation is concerned with the behavior of simply-supported skew 
I-beam bridges. This type of structure is used extensively in this country for spans 
up to about 100 ft. It consists of a reinforced concrete slab in the form of a parallelo- 
gram, supported on a series of five uniformly spaced, flexible girders which are 
parallel to the direction of traffic and are supported on abutments which, in general, 
are not perpendicular to the girders. 

Analyses have been made for a number of structures of this type for various 
proportions and relative stiffnesses of the slab and the beams, and for different 
angles of skew. From the results of these analyses general relations have been ob- 
tained for the rational design of this type of structure. 

The mathematical expressions governing this problem are: the differential 
equation for flexure of slabs, 

(1) NV*(Vw) = p 
where 
EI 
1 — 2 


, 


a constant related to the properties of the slab, 
id rd 
"ont oy 
the Laplace operator, p is the intensity of the external load and w is the unknown 
deflection of the slab at a point defined by the rectangular coordinates x and y; 
the differential equations for the moments M,, M,, M,, are 
ow Pw 
M, 7 Ee ax? + =) 
rw Pw 
ay? +h =) 


(2) 
M, = - 


ow 
M,, = —N(1 — p) 8.0” 
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the differential equation for the deflection w of a beam 


dw _ 


(3) Euly = 0 


and the equation for the moment of the beam 


d*w 
(4) M, = Bolo 


where E;,J, is the flexual rigidity of the beam and q is the load per unit length of 
the beam. 

The problem treated in this report is the determination of w which satisfies 
equations (1) and (3) and which in addition satisfies the boundary conditions of 
the problem. Once w, the deflection of the structure has been obtained, the bending 
moments (from which the stresses are evaluated) may be determined from equations 
(2) and (4). 

Because of the nature of the boundary conditions, an exact solution of the 
governing differential equations does not appear to be possible. In this investiga- 
tion the problem was analyzed by the method of finite differences using a network 
of points in skew coordinates, in which the elemental mesh is a parallelogram 
formed by lines parallel to the sides of the skew panel. A total of 63 node points 
were considered. By writing a difference equation for each point, there results, in 
general, a set of 63 simultaneous algebraic equations, the solution of which gives 
the values of the deflection at the nodal points. In this study, the load on the bridge 
was broken into symmetrical and antisymmetrical components. In this way the 
analysis of each particular problem was reduced to the solution of two sets of equa- 
tions, one set having 32 equations and the other 31. The complete solution was 
obtained by adding the solutions for the symmetrical and antisymmetrical com- 
ponents of loading. The choice of the number of nodal points considered was gov- 
erned by the capacity of the ILLIAC which could solve 39 simultaneous equations 
at the time this investigation was conducted (1954). (At present the ILLIAC can 
solve as many as 142 simultaneous equations.) 

A total of 18 structures having different geometrical and physical characteristics 
were analyzed in this investigation. These analyses necessitated the solution of 306 
sets of simultaneous equations (each having 31 or 32 unknowns). It took the ILLIAC 
3.5 minutes to solve each set, or about 17 hours to solve the 306 sets. 

The coefficients of these equations were functions of the angle of skew, ¢, the 
ratio of the sides, b/a, and the relative stiffness of the slab and the beams, H. As 
the evaluation of these coefficients on a desk calculator and the preparation of 
tapes on which numerical values of about 1000 coefficients were to be punched for 
each set of equations would be most laborious and time-consuming, a special pro- 
gram was prepared which permitted the evaluation of these coefficients within 
the ILLIAC for specified values of ¢, b/a, and H. 

Contour maps, giving the coefficients for the bending moment at midspan of 
the second beam due to a single concentrated load applied at any point on the bridge 
were obtained. For other loads the value of the bending moment is obtained by 
multiplying the coefficients given on the diagram by Pa, where P is the magnitude 
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of the concentrated load and a is the span of the bridge. The bending moment due 
to multiple wheels can be found by superimposing the effect of each load. 

From influence surfaces of this type, maximum live load moments were obtained 
for standard trucks and the results were used to develop general relationships which 
provide a more rational basis for the design of skew I-beam bridges. 


A. H. T. 


50[S ].—Pamir M. Morse & Huseyin Yiimaz, Tables for the Variational Deter- 
mination of Atomic Wave Functions, The Technology Press, Massachusetts 
Institute of Technology, Cambridge, Mass., 1956, 85 p., 30 cm. Price $2.50 
plus 15c for mailing. 


This set of tables should prove useful in the many instances when one desires 
an atomic wave function sufficiently accurate to distinguish between multiplets 
but simpie enough to be expressed in closed form for easy use in subsequent calcu- 
lation. Such functions should be welcomed when a modern high-speed computer is 
either not available or practical. The wave junctions in point are to be constructed 
as Slater determinants of hydrogen-type orbitals; the energy of the atom is com- 
puted and minimized with respect to the parameters contained in the orbitals in 
order to obtain the values which characterize the most acceptable function. This 
yields a wave function which neglects inter-election correlation. 

The tables embrace ten hydrogen functions: the one 1s, the one 2s, the three 
2p, and the five 3d orbitals. Each group is characterized by one parameter. Ex- 
pressions are given for the kinetic energy, potential energy, and exchange energy 
contributed by any combination of such electrons to the total energy. It is these 
expressions which are tabulated as function of the parameters so that the minimiza- 
tion can be conveniently carried out on a desk calculator. Even those energies 
which are trivia! functions of the parameters are tabulated. The functions tabulated 
are expressible in terms of combinations of hypergeometric functions and were 
computed on IBM CPC equipment. No discussion of numerical methods or pro- 
cedures is given. All tables are given to at least 5 decimal places. A recipe for carry- 
ing out the construction and minimization procedures along with a rather complete 
example is given. No discussion of the errors to be expected in this type of varia- 
tional wave function is included nor is a comparison between a worked out function 
of this type and other better or worse types of function given for some particular 
case. 

James N. SNYDER 
Digital Computer Laboratory 


University of Illinois 
Urbana, Illinois 


51[S].—A. J. Fercuson & A. R. RUTLEDGE, Coefficients for Triple Angular Corre- 
lation Analysis in Nuclear Bombardment Experiments, CRP-615, Atomic Energy 
of Canada Limited, 1957, xvi + 469 p. 27 cm. Price $6.00. 


Here is an ambitious table, in a somewhat difficult format, for the nuclear bom- 
bardment problem of studying nuclear states by measuring the angular correlations 
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among the directions of an entering particle and two gamma rays emitted in cas- 
cade. The tabulated expansion coefficients depend on some thirteen parameters, 
which explains most of the difficulty in the tabular format. The remaining difficulty 
is attributable to the suppression of zero coefficients. 

The introduction is clear enough for students of the problem. It is quite opaque 
to others, but adequate references are given. There is a well-marked sample page. 

The ranges covered are as follows, where s is the channel spin, the /’s are allowed 
incident orbital momenta, the J’s are the spins of the first and second states, J is 
the spin of the final state, and the L’s are the allowed multipolarities of the gamma 
rays. s: 0(1/2)9/2, with 8/2 omitted. /: 0(1)3. J: 0(1/2)5. I: 0(1/2)7/2. L: 1(1)2. 


R. Lazarus 


Los Alamos Scientific Laboratory 
Los Alamos, New Mexico 


52[S].—A. V. H. Masxket, R. L. Mackin, & H. W. Scumitt, Tables of Solid 
Angles and Activations: 1. Solid Angle Subtended by a Circular Disc; 11. Solid 
Angle Subtended by a Cylinder; 111. Activation of a Cylinder by a Point Source, 
ORNL-2170, United States Atomic Energy Commission, Technical Information 
Service Extension, Oak Ridge, Tennessee, Nov. 1956, 18 + 57 p., 27 cm. 


There are three tables, furnishing approximations to the activation integral J = 
Jy dv exp(— 1/d)/R? in the limiting cases of //A > 1, (J = Q, the solid angle) and 
i/\ K 1 (UI = Sy dv/R*), where V is a right circular cylinder. The application in 
view is the activation of a cylindrical detector by a point radiation source (or a 
point detector by a uniform cylindrical source). R is the total distance, and / the 
distance inside V, from the point source to dv; is the mean free path in V. A cor- 
rection and error estimation procedure is given for //A small but not negligible. It 
is suggested that some correction can be made for elastic scattering of the inci- 
dent radiation by using certain values in Table III to estimate an effective in- 
creased path length. The introduction is reasonably complete; additional mathe- 
matical details can be found in [1] [2]. 

Table I contains the solid angles subtended by a circular disc of unit radius at a 
point in a plane & units above the disc plane and p units from the disc axis. The 
entries are for kh = 0.1(0.1)5.0(0.2)10.0(1.0)20.0(2.0)100, and p = 0.0(0.1)0.9 
(0.05) 1.0(0.1)6.2(0.3)6.5(0.5)16(1.0)36. 

Table II contains the solid angles subtended by a right circular cylinder of unit 
radius and height ’ units at a point, P, in the base plane p units from the axis. 
Solid angles subtended at points in a plane cutting the cylinder can be found by 
addition of the entries for the upper and lower parts of the cut cylinder; those at 
points in a plane below the base must be found by the difference of two values, 
corrected by addition of the base term from Table I. The entries are for the same 
values, excluding p < 1. 

Table III contains the values of the integral J = [y dv/R* over the cylinder of 
Table II, where R is the distance from P to dv. The entries are for the same values 
as in Table I, except for the inclusion of p = 0.85 and the omission of p = 6.2. 

All entries are rounded to six significant figures and are of the form 0872425003, 








O4 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 
meaning 0.872425 x + 10+*, + or — for o printed as O or F, respectively. The 
character size is unfortunately rather small. 


RoGER LAZARUS 


Los Alamos Scientific Laboratory 
Los Alamos, New Mexico 


1. A. V. H. Masxet, R. L. Macguin, & H. W. Scumitt, ‘Cylinder radiation and activation,” 
Rev. Sci. Insir., v. 28, 1957, 3 189. 
i 


2. A. V. H. MAsKEr, “ angle contour integrals, series, and tables,” Rev. Sci. Instr., v. 28, 
1957, p. 191. 


53(W, Z].—C. E. SHaNNon & J. McCartuy, Editors, Automata Studies, Annals of 
Mathematics Studies No. 34, Princeton University Press, 1956, viii + 286 p., 
23 cm. Price $4.00. 


This collection of papers is classified under three headings: (1) Finite Automata, 
(2) Turing Machines, and (3) Synthesis of Automata. The papers by Kleene, 
Moore, and von Neumann under (1) and those by Shannon, Davis, and De Leeuw, 
Moore, Shannon, Shapiro under (2) obtain carefully worked out results of mathe- 
matical interest. The papers under (3) are concerned with simluating aspects of 
animal behavior. 

S. C. Kleene in “Representation of Events in Nerve Nets and Finite Auto- 
mata” obtains characterizations of the behavior of automata first postulating past 
time as infinite and then postulating the contrary. The principal result for auto- 
mata with infinite past is that for every automaton there is an automaton without 
cycles which behaves in the same way (this follows from theorems 1, 6 and lemma 
1) and this behavior can be described by a (finite) set of tables all of the same 
finite dimensions. The top row of a table specifies that input state at any time ¢, 
the next row the input state at ¢ — 1, etc. These tables specify exactly those se- 
quences of input states to which the automaton responds. Conversely, given a set 
of such tables there exists an automaton whose behavior it represents. This leads 
to the more significant and more complicated theory of automata with finite past. 
The notion of “regular event” is defined and it is proved that automata (with 
finite past) can detect all regular events and only these. The class of k-input regular 
events (those which k-input automata detect) could be defined as the smallest class 
of sets of sequences of k-tuples of zeros and ones which includes all the unit sets 
and contains a U 8, a-8 and a* whenever the class contains a and 8. Here a-8 
is the set of sequences obtained by prolonging any sequence of a with a sequence 
from 8 and a* = a U a-a U a-a-a U---. Concrete examples of automata are 
obtained by interconnecting (models of) neurons. The resulting configuration is 
called a nerve net. 

The notion of indistinguishability is central for E. F. Moore “Gedanken- 
Experiments on Sequential Machines.” With each initial internal state q of an 
automaton A a transformation T,4 (g) may be defined which when applied to a 
finite sequence of input states has as value the determined sequence of output 
states. The author’s concept of indistinguishability of two states g and q’ of auto- 
mata A and A’, respectively, may be formulated as the requirement that 7, (¢) = 
T4/(q’). The author’s concept of indistinguishability of two automata seems odd. 
It may be formulated as follows: there exist functions f, g such that T,(q)[¢] = 
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Ta (f (i,q) C4] and Ta) 4] = Ta(eG, 7)) [4], where i varies over finite se- 
quences of input states, g over internal states of A and q’ over internal states of 
A’. If A is strongly connected (i.e., for any ordered pair of internal states there 
exists a sequence of input states which takes the first internal state into the second) 
then f may be chosen so as not to depend on i. Thus, when applied to strongly con- 
nected machines, the notion of indistinguishability becomes natural. The author 
shows that among all automata indistinguishable from a given strongly connected 
one there exists an (essentially) unique automaton with a minimal number of in- 
ternal states. Besides qualitative theorems of the kind described above, the author 
obtains quantitative theorems typified by the following result. If an automaton 
A has the property that each two of its internal states are distinguishable then for 
any pair g, ¢’ of internal states there exists a sequence s of input states of length 
n — 1 such that T, (g) ¥ Ta(q’), where n is the number of internal states of A. 

In ‘‘Probabilistic Logics and the Synthesis of Reliable Organisms from Unreliable 
Components,” J. von Neumann introduces a technique for constructing reliable 
finite automata out of elements whose probability of improper operation is positive. 
The technique is called multiplexing. Instead of using, say, an element with two 
input lines and one output line, one uses a network of elements (a generalized ele- 
ment) with two bundles of N (N is usually a large positive integer) input lines and 
an output bundle of WN lines. It is assumed a certain fiduciary level A is fixed (the 
same for all bundles) with 0 < A < 4. The bundle is stimulated if at least (1 — A) 
N lines of the bundle are stimulated. The bundle is unstimulated if at most AN 
lines of the bundle are unstimulated. Otherwise the bundle is in an uncertain state 
and the network malfunctions. The generalized element is designed to realize the 
same function as the basic element. In some cases the input bundles of a generalized 
element may each be stimulated or unstimulated while the output state is un- 
certain. To minimize this effect the author introduces, besides an “executive organ” 
intended to realize a given function, a “restoring organ” intended to transform 
(with high probability) an input bundle with stimulation level near zero (one) to 
an output bundle with stimulation level nearer zero (one). To perform this re- 
storative function it must be assumed that the stimulated (unstimulated) lines are 
distributed at random in a bundle. It is not considered, in detail, how this random- 
ness can be effected (and maintained). The author obtains detailed quantitative 
results for networks constructed out of a generalized Sheffer stroke element. If 
A = .07 and e < .0107 (e is the probability of error of a single line Sheffer stroke 
element independent of the state of the network) and if the generalized Sheffer 
stroke element is constructed out of the single line Sheffer elements, then the prob- 
ability of malfunction of the network may be made arbitrarily small by taking N 
sufficiently large. In certain examples, rough calculations show that with « = .005, 
N must be about 20,000 to 25,000 to obtain “‘good” reliability. 

In “Some Uneconomical Robots,” J. T. Culbertson shows how to realize arbi- 
trary truth functions by interconnecting neuron-like elements. He, apparently, 
regards his robots as “‘existing”’ for a finite amount of time only. Since he seeks (in 
this paper) to design only “uneconomical” robots he does not consider neuron nets 
with cycles. 

M. L. Minsky in “Some Universal Elements for Finite Automata” is con- 
cerned with constructing automata, out of certain types of elements, which will 
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realize certain, rather restricted, kinds of behavior. A ‘‘stimulus’’ is a finite set of 
(input) states at certain discrete instants of time. A “response” is analogous. A 
“response function” is defined on a finite set of stimuli with values which are re- 
sponses. An automaton “‘realizes’’ a response function, f, if whenever it receives a 
stimulus, s, which is in the domain of f, its response is f(s). The author’s main 
result is that a wide class of response functions (which he describes) can be realized 
(with some delay) by automata constructed from disjunctions, conjunctions, and 
any “non-monotonic” element. Inasmuch as the domain of the,response functions 
is finite, his problem is similar to that of finding complete sets of primitives for truth 
functions. The problem is complicated by the fact that his primitives all have a 
“built in” unit delay. 

M. D. Davis in “A Note on Universal Turing Machines’’ proposes a definition 
of “universal Turing machine.” Let Z be a Turing machine. An instantaneous de- 
scription of Z is a pair <q, E> consisting of a state g of the reading head of Z, 
and an expression E with a marked place, composed of symbols of Z. Let 5z be the 
set of Gédel numbers of all instantaneous descriptions <q, E> such that if E is 
printed on the tape of Z, and the reading head is put in state Q scanning the marked 
place of E, then Z will eventually halt. Now the machine Z is defined to be uni- 
versal if the set 5z is recursively enumerable and has the further property that to 
every recursively enumerable set R there is a recursive function o(x) such that 
R = {x | o(x)e 5z}. Using methods of the theory of recursive functions (particularly 
Kleene’s T-function) the author shows: 1. How an arbitrary machine Z can be so 
encoded on a universal machine U, as to make U carry out the computation per- 
formed by Z (this, of course, is relative to particular conventions concerning the 
use of machines to compute functions.) 2. That a Turing machine which would 
(relative to the adopted convention concerning computations) ordinarily be called 
universal (on more or less intuitive grounds) is universal in the precise sense of the 
definition. 3. That for a universal machine U there is a machine C which performs 
the encoding mentioned in 1, and which is not itself a universal machine. These 
results are intended to make it plausible that the given definition of universality is 
reasonable. In particular 3 shows that the encoding necessary to make a universal 
machine perform the computation of another machine is simple in a precise sense. 

C. E. Shannon in “A Universal Turing Machine with Two Internal States” 
investigates the number of tape symbols and internal states required for universal 
Turing machines. If U is a Turing machine, let st(U) denote the number of states 
which the reading head of U can take on, and let sy(U) denote the number of sym- 
bols which the reading head can print on the tape. Given a machine A to construct 
a machine B which acts essentially like A and such that st(B) < st(A) (or sy(B) < 
sy(A)). Part of this problem, of course, is to adopt reasonable interpretations of 
“essentially acts alike.’’ The specific results obtained are: 1. To every machine A 
one can construct a machine B such that st(B) = 2, sy(B) = sy(A) + 4sy(A) -st(A), 
and such that B essentially acts like A in the following sense: Let A:, ---, An be 
the symbols of A, let A:, --- , An, Bi, --- , B, be the symbols of B, let T> be any 
tape which carries symbols A; only, and let ¢ be a marked cell of To. By properly 
changing the symbols in ¢ and one of the neighboring cells to symbols B; one can 
obtain a tape T,’ such that if the cells c of T>) and Ty’ at time 1 are placed under the 
reading heads of A and B respectively, there will be a sequence 1 = h < fp < fs --- 
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of time instances such that for any i the tape of A at time i carries the same sym- 
bols in corresponding cells as the tape of B at time ¢;, except for the cells under the 
reading heads and one of the neighboring cells. 2. There is a universal Turing ma- 
chine U such that st(U) = 2. This is clearly a corollary to 1. It is shown also that 
a machine W with st(W) = 1 cannot be universal. 3. To every machine A one can 
construct a machine C such that sy(C) = 2 and st(C) < 8-st(A)-sy(A), and such 
that C acts essentially the same as A in roughly the sense that if the symbols of A 
are coded by use of the two symbols of C, then the machine C acts on the coded 
tape T’ like the machine A acts on T. These results suggest that it is possible, 
within limits, to exchange symbols for states and vice versa, with about a factor 8 
of change in the product st(M)-sy(M). It is conjectured that the loss factor 8 
could be decreased if “essentially acts like” is interpreted as “having the same 
calculating abilities in the large” rather than as “behaving alike in a local sense.” 

K. de Leeuw, E. F. Moore, C. E. Shannon, and N. Shapiro in an article entitled 
“Computability by Probabilistic Machines” seek to answer the question: “Is there 
anything that can be done by a machine with a random element but not by a de- 
terministic machine?” A machine is defined in such a way that if its input is a 
function f on the natural numbers to zero, one that is an infinite sequence of zeros 
and ones then, if f is recursive, its output may be interpreted as an enumeration of 
a recursively enumerable set of natural numbers. For f fixed, given a recursively 
enumerable set there exists a machine which will enumerate it. If the machine is 
now modified so as to receive as input a function f such that f(m) = 1 with proba- 
bility p (independent of n), 0 < p < 1, the machine is called a p-machine. A set 
is strongly enumerated by a p-machine if the probability that the set will occur as 
the output of the machine is non-zero. A set is weakly enumerated by a p-machine 
if the probability that x occurs in the output is greater than 4 if and only if x is 
an element of the set. The main result of the paper may be stated as the equivalence 
of the following three statements concerning a set S: 


1) S is enumerated by some machine whose input is the binary expansion of p. 
2) S is strongly enumerated by some p-machine. 
3) S is weakly enumerated by some p-machine. 


From this it follows readily that if p is computable, i.e., the sequence representing 
p is recursive, then p-machines can only strongly (weakly) enumerate recursively 
enumerable sets while if p is not computable some p-machines will strongly (weakly) 
enumerate non-recursively enumerable sets. The authors define a computable 
stochastic machine (c-s machine) as an object having a countable (finite or de- 
numerable) set of states, a distinguished initial state, a countable set of output 
symbols and an effective function which when applied to the sequence x1, %2, ---, 
%n41, S, m (n > 0) where the x’s are states, s is an output symbol and m is a positive 
integer, yields the first m bits of the binary expansion of the probability that the 
machine will “next” be in state x4: and “print” s if it has, up to the “present,” 
been in the states x, x2, ---, %,- The authors show that the c-s machines are 
equivalent (in a precisely defined sense) to p-machines with p-computable. From 
this result it follows that the outputs of c-s machines represent exactly the re- 
cursively enumerable sets. 

J. McCarthy in “The Inversion of Functions Defined by Turing Machines” 
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considers the task of designing machines to solve well-defined intellectual problems. 
This task is identified with inverting functions defined by Turing machines. (The 
identification could be made plausible by referring to formal theories.) Suppose 
fm is the function calculated by the m-th Turing machine. Let g(m, m) be called 
an inversion-function if (1) f,(g(m, m)) = m in case there exists an x such that 
fm (x) = m, and (2) g(m, n) is undefined otherwise. The problem is to construct 
machines which “efficiently” calculate values of an inversion-function. Suggestions 
are made as to what efficiency should be taken to mean. No theorems are proved. 

The papers of W. Ross Ashby, “Design for an Intelligence Amplifier,” and 
D. M. MacKay, “The Epistemological Problem for Automata,” are parallel in 
many respects: i) each is concerned with synthesizing automata which can per- 
form functions until now reserved to the mental processes of the higher animals, 
e.g., solving difficult problems or forming concepts, ii) each presents descriptions 
of specific mechanisms and then suggests ways of constructing a general device 
using principals embodied in the specific device, iii) each makes random fluctua- 
tions a basic part of the mechanism (although in different ways), iv) each makes 
comparison a certral part of the process (again in different ways). Ashby’s process 
requires that the universe of potential solutions to a given problem be specified 
either by means of a model or by letting the mechanism operate directly on the 
appropriate part of the real world. The mechanism then makes random seiections 
from the potential solutions and compares them with the conditions to be satisfied. 
The basic unit of the process proposed by MacKay consists of a comparator, a 
“control-signed generated and mixer’ M, and an effector. The comparator com- 
pares input (from outside the unit) with effector output and produces a signal 
which gives a measure of different between the two corresponding to the error 
signal in a servomechanism. This error signal is used to activate, via M, a new 
effector output. MacKay’s process then consists of a hierarchy of such units with 
the output of M at one level acting as input to the next (as well as activating the 
effector at its own level). The input of the lowest level unit comes from the environ- 
ment of the whole device. MacKay concludes: “...an automaton designed on 
statistical principles, which can evolve an internal organizing routine to respond 
adaptively to regularities of its sensory input, is capable in principle of developing 
its own symbols for concepts of any order of abstraction, including metalinguistic 
concepts, without prior instruction ...”” MacKay’s mechanism seems to be more 
sophisticated than the one proposed by Ashby in that it allows for learning. That is, 
Ashby’s upon being given the same problem a second time can be expected to take 
as long as before; not so MacKay’s. MacKay’s mechanism when presented with a 
problem of searching for regularities in the input at any level will improve its per- 
formance after several runs in the sense of becoming more efficient at reducing the 
error signal. 

In “Conditional Probability Machines,”’ Albert M. Uttley discusses a machine 
whose response to stimuli are similar in some ways to those of animals and gives a 
set of sufficient design requirements for such a machine. The design incorporates 
the mathematical principles of set inclusion and conditional probability, i.e., 
Bayes’ theorem. Each unit computes the probability of occurrence of the stimulus 
J which corresponds to it when another stimulus K occurs. This probability de- 
pends on the past history of joint occurrences of stimuli. “Temporal and Spatial 
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Patterns in a Conditional Probability Machine” discusses the problem of dis- 
tinguishing between temporal patterns by a method of delayed inputs. The dis- 
cussion on spatial patterns is rather vague. 


The University of Michigan 
Ann Arbor, Michigan 


54[X].—Her Mayjesty’s NauticaL ALMANAC OFFICE, Interpolation and Allied 
Tables, Her Majesty’s Stationery Office, London, England, 1956, 80 p. 24.5 
cm. Price 5 s. net. 


This is a new edition of a booklet with the same title, first issued in 1936 as a 
reprint of certain pages of the British Nautical Almanac, and which was frequently 
reprinted, sometimes with additions. It derives largely from the experience, over 
several decades, of the late L. J. Comrie, Superintendent of H. M. Nautical Al- 
manac Office, his successor, D. H. Sadler, and their disciples. In its present form 
it provides an excellent practical introduction to a large part of classical numerical 
analysis and can be thoroughly recommended as an inexpensive compact text, 
complete with tables. The typography is excellent, and the reviewer’s only complaint 
concerns the redundant decimal points in Table 6. 

I. Part I contains a section on elementary computing principles, which is fol- 
lowed by an account of the principles and practices of interpolation and some 
remarks on numerical differentiation and integration. There are worked examples 
showing the use of the tables in the second part. 

II. Part II contains various tables of Bessel and Everett coefficients, some 
critical and some direct. To see the scope of these tables it will be worth while sketch- 
ing formally current (British) practice in the interpolation. The latest developments, 
of applying the ideas of Chebyshev economization to the problems of interpolation 
are discussed briefly in Part III. They will be discussed briefly in a forthcoming 
review of the booklet of L. Fox [5] in this journal. 

One method is to truncate the Bessel formula 


fo = fo + Pdin + Ba( do? + 5°) + Bsdy® + Ba(dot + d:*) + --- 
at an appropriate place. Here 
B:(p) = B,(1 — p) = 36(p — 1)/21, 
Bs(p) = —B(1 — p) = o(b — 1)(b — 4)/3}, 
Bp) = Bi — p) = 3(p + 1)P(b — 1) (b — 2)/41. 
The modified form of this, with fourth and sixth differences thrown back is 
¥ fo = fot pore + Be(5no + bm) + Bady* 
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or, if 5, 5 are varying slowly, 


fo = fo + porn + Ba(Smo + 8%:) + Bs(Si1 — Fro) 
where 


bn? = & — 0.1845 + 0.0386°. 
Otherwise we use the Everett formula 

fo = (1 — P)fo + Pf: + (Eade? + Fei?) + (Eads + Fads’) + --- 
truncated appropriately. Here 

E,(p) = Fa(1 — p) = —p(p — 1) (p — 2)/3! 

E(p) = Fi(l — p) = —(b + 1)P(b — 1) (@ — 2) (b — 3)/S!. 
Various modifications of this are possible. Corresponding to * we have 

fo = (1 — b)fo + fi + (Laing + Fd) 


with 52, defined as above. If we throw the sixth and eighth differences back into 
the fourth according to 


bm? = 54 — 0.2075 + 0.0458, 
keeping the second differences as they are, we have 
fo = (1 — P)fo + bf: + (Hah? + Fabi®) + (Eads, + Fadl). 

We can also use a simultaneous throw-back of 6, 5°, 5 into 6*, 5* according to 

5,7 = & — 0.013125 + 0.00435* — 0.0015" 

Smt = 5¢ — 0.278278 + 0.06858 — 0.0165" 
and use 
m fe = (1 — b)fo + dfi + (Eating + Fain) + (Easing + Fudan): 
Finally, we can use a generalized throw-back 

6,2? = &# — 0.18454 + 0.0380825° — 0.008308* + --- 
1074 = 0.015 — 0.0027835* + 0.000685* + --- 


(here +‘ is a modified fourth difference (divided by 1000) ). With this notation we 
can use 


fo = (1 — d)fo + Phi + (Eadno + Fabia) + (Mert + Neri‘) 
or, if +“ is varying slowly 


fo = (1 — Dd) fo + dfi + (Eabng + Fed.) + Talo + vi‘) 
where 


M, = 1000(E, + 0.1842), N, = 1000( F, + 0.184F;), 
T, = 3(M, + M). 
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In recent British tables, e.g., those of the Airy integral [1] and the Weber 
functions [2],-the y’s are tabulated where necessary. Modified differences have 
been used in many tables since about 1930. 

The twelve tables given, which we shall now describe in detail, are sufficient 
for most applications of these methods of interpolation. 

1. First-difference correction piip. This gives to the nearest integer, the value of 
piin for p = .1(.1)9 and 5:2 = 1(1)128. The table is arranged in blocks corre- 
sponding to eight consecutive values of 51/2 and on each block is given the appropriate 
contribution of the second decimal of 9, i.e., 51/2 for p = .01(.01).09 for an average 
512 (2 pages). 

2. Second-difference correction B2(i? + 46;*). This gives to the nearest integer 
the values of —Bz(do? + 6,*) for p = 0(.01).5 and for 5? + 8,2 = 10(5)200. 

3. Bessel coefficients Bz, B;, By. There is a critical table of Bz, proceeding in units 
of the third decimal, with p given to 4D, together with 3D values of B; and B, 
at interval .1. 

4. Third-difference correction B;6,*. This gives to the nearest integer, the values 
of Byd;* for p = 0(.01).5 for 5;* = 100(100) 1000. 

5. Fourth-difference correction B,(o* + 5;*). This gives to the nearest integer, 
the values of By( dot + 4:*) for p = 0(.01).5 and dot + 5,4 = 100(100) 1900. 

6. Second-difference coefficients B:, E,, F2. This is a critical table, proceeding 
in units of the fourth decimal; is given to 5D. Such a table is feasible since 2B, = 
E, + Fz. (7 pages). 

7. Bessel coefficients Bz, B;, By. This gives B; to 6D, B; to 5D and B, to 4D for 
p = 0(.001).5. First differences of B, are given (5 pages). 

8. Throw-backs 0.184 5‘ and 0.0385*. This is a critical table, proceeding in units, 
with integral values of 5‘ and 4 up to 5 = 1002, 5° = 250. 

9. Throw-backs 0.2075* and 0.0458*. This is a critical table of 0.207 5° proceeding 
by 5 units up to 220, with integral values of 5*, and also giving 0.2074 to the nearest 
unit for 6* = 1000(1000) 10000. There is a critical table of 0.0455* proceeding by 
5 units up to 40, giving integral values of 6°. 

10. Residual fourth-difference correction Ti(70' + i). This gives the values 
of the correction, to the nearest integer, for p = 0(.02).5 and (yo' + v:‘) = 
001 (dot + 8:4) = 2(2)24. 

11. Short table of Bessel and Everett coefficients. This gives B, for vy = 2 (1) 11 
and E2,, F2, for »y = 1(1)5, to 6D for p = .1(.1).9. 

12. Coefficients for Everett's formula. This gives Ex, and F; to 7D, M, and N, 
to 3D, T, to 2D, E,and F, to 6D, E,and F;, to 5D, all for p = 0(.001).5 (10 pages). 

III. Part III is a collection of relevant formulae, subdivided into eight sections, 
some of which, e.g., that on the solution of ordinary differential equations, contain 
much valuable commentary. We deal briefly with each section. 


A. Operators. The finite difference operators are defined and their interrelations 
are recorded. 


B. Differences. The notation is discussed and there are short tables useful in 
connection with checking by differencing. 

C. Interpolation. In addition to the material summarized in our account of 
Part II (supplemented, of course, by working limits for applicability of the 
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various formulae) there are sections on the magnitude of the Bessel co- 
efficient, on interpolation to halves, on the Lagrangian interpolation method 
and Aitken’s scheme for carrying it out and on adjusted divided differences. 
There is a discussion of the central difference formulae of Stirling and 
Steffensen for |p| < 1/2. A polynomial form of ** is given and there is a 
discussion, with tables, of the use of economized polynomials for interpola- 
tion [3, 4, 5]. The problem of sub-tabulation is not discussed: it will be the 
subject of a companion booklet [6] (5 pages). 

D. Differentiation. This includes expressions for the r-th durivative; r=1(1)10, 
evaluated at a tabular point, or at a halfway point, with coefficients given 
both as rational numbers and as decimals. Expressions in terms of ordinates 
and in terms of forward and backward differences (for use near the ends 
of a table) are also given, in addition to expressions for differences in terms 
of derivatives. 

A table of coefficients of the differences in the expressions for the first 
three derivatives at a non-tabular point p is given, for p=0(.01)5. These 
expressions are obtained by differentiating the Bessel formula. 

The formulae appropriate for tables with modified differences are also 
given (5 pages). 

E. Integration. Various formulae for evaluating definite integrals in terms of 
various differences, derivatives, and ordinates are given. In particular, the 
coefficients (of the ordinates) for use in Gregory’s formula, the correspond- 
ing central difference formula, and the Gaussian quadrature are given to 
8, 8, and 10 decimals. 

The case of simple and double integration to non-tabular points is dis- 
cussed and suitable tables are provided. There is a selection of starter, 
predictor, and corrector formulae for use in the numerical integration of 
differential equations (5 pages). 

F. Differential equations. This account is restricted to initial value problems 
for ordinary differential equations. Five methods are mentioned: (1) local 
Taylor series, (2) Runge-Kutta, (3) Adams-Bashforth, (4) Central differ- 
ence extrapolation methods, (5) Deferred-correction method (5 pages). 

G. The precision of calculation. This includes a table giving the absolute errors 
and the variances of errors resulting from certain arithmetical operations 
on quantities subject to error, e.g., 2-th difference, or the result of a direct 
or inverse interpolation. There is a short table of ¢(x) = e@”%2/\/2x and 


P(t) = 200 i ¢(x)dx, the percentage area of the tails of the normal curve. 
t 
H. Constants. A two-page table of constants to a reasonable number of places 
is given. 
The booklet closes with a one-page bibliography. 
Joun Topp 


California Institute of Technology 
ees California 


1. J.C. P. Mrtier, The Airy Integral, Part-volume B of British Assoc. Adv. Sc., Math. Tables, 
Cambridge University Press, 1946. 
2. NATIONAL PHysICAL LABORATORY, Tables of Weber Parabolic Cylinder Functions, Her Ma- 
jesty’s Stationery Office, London, 1955. 
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3. C. W. CLensHaw & F. W. J. Over, ‘‘The use of economized polynomials in mathematical 
tables,” Camb. Phil. Soc., Proc., v. 51, 1955, p. 614-628. 

4. ZpEnéK Korat, Numerical Analysis, John Wiley and Sons, Inc., New York, 1955. 

5. L. Fox, The Use and Construction of Mathematical Tables, National Physical Lab., Math. 
Tables, v. 1, Her Majesty's Stationery yen 1956. 

6. HER Majesty's NavuticaL ALMANAC OFFICE, Sublabulation, A Companion Booklet to Interpo- 
lation and Allied Tables, Her Majesty s Stationery Office, London, 1958. 


55[(X ].—J. L. Syncz, The Hypercircle in Mathematical Physics: A Method for the 
Approximate Solution of Boundary Value Problems, Cambridge Univ. Press, 
New York, 1957, xii + 424 p., 24 cm. Price $13.50. 


There are, roughly speaking, two lines of approach to the problem of the nu- 
merical solution of partial differential equations: The methods of finite differences 
and the methods based on ideas of Functional Analysis or the Calculus of Varia- 
tions. There exist many connections between the theories of the two procedures, 
but most of the work done in this field falls distinctly into one of those two cate- 
gories. 

The present book gives a complete and detailed account of a variational method 
developed, during the past ten years, by the author and several other mathe- 
maticians. The underlying ideas are connected with the names of Thomson, Ritz, 
Rayleigh, and Trefftz, but the terminology, the mode of the presentation and the 
numerical techniques are original contributions of the author and a group of more 
recent mathematicians. 

Here is, in a nutshell, a description of the hypercircle method: The problem is 
the approximate numerical solution of a boundary value problem of elliptic type 
in a region V. A suitable Hilbert space F is introduced such that to every function 
in V there corresponds a point in F. If F is suitably defined the functions that 
satisfy the boundary conditions and the functions that satisfy the differential 
equation define two linear subspaces L’ and L’’ of F, respectively, that are or- 
thogonal. They intersect at the point S corresponding to the solution of the bound- 
ary value problem. The position of S is unknown, but if one chooses r points of 
L’ and s points of L’’ one can determine, by linear combination, the hyperplane 
L', of dimension r in L’ through the first r points and, similarly, a hyperplane 
L”, in L”’ of dimension s. It is a simple matter to calculate the end points V’, V”’ 
of the shortest straight segment joining L’, and L”,. If these r + s points, i.e., the 
corresponding functions in the domain V, are skillfully chosen the segment V’V” 
in Hilbert space will be short, and its midpoint C will be a good approximation to 
the true solution S. The error can be appraised by observing that S lies on the 
hypersphere with diameter V’V’’, thanks to the orthogonality of L’ and L”. Also, 
the r + s known points of F limit the position of S to a known hyperplane defined 
by r + s linear equations. The intersection of the hypersphere with this hyper- 
plane is the hypercircle from which the book takes its name. Simple geometric 
arguments based on the triangle inequality lead then to convenient inclusion 
theorems by means of which the error in the Hilbert norm can be appraised. Point- 
wise bounds on the error are then obtainable from an inequality derived with the 
help of a fundamental solution of the differential equation. 

The larger part of the book is devoted to applications based on the technique 
of so-called pyramid functions. These are functions that are piecewise linear in V 
with certain discontinuities admitted along straight lines (or planes, in three di- 
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mensions). The author shows how approximating subspaces L’,, L’”’, can be con- 
structed by linear combination of such pyramid functions. Approximation by piece- 
wise linear functions is a natural bridge between the variational and the finite 
difference methods. This aspect, which in the reviewer’s opinion deserves further 
study, is only briefly discussed in the book. 

This book is intended for engineers and physicists. The author is aware of the 
fact that many readers in this group shy away from abstract arguments in Func- 
tional Analysis. He therefore develops his topic very slowly and carefully assuming 
no more mathematical knowledge than a good course in Advanced Calculus will 
impart. The systematic appeal to geometric intuition assisted by a large number of 
diagrams, as well as many numerical examples that are described in great detail, is 
intended to make the subject palatable to the applied scientist. On the other hand, 
the proofs for all statements are mathematically complete. 

The author is a skillful expositor and writes in a very attractive, lively style. 
Whether he has steered the safest course between the Scylla of excessive mathe- 
matical conciseness and the Charybdis of tiresome repetitiousness in his explana- 
tions is hard for this reviewer to decide. A book of 424 pages seems a rather lengthy 
account of one particular numerical method. It is true that the reader will learn 
much that is of more general importance, such as the most elementary properties 
of Hilbert space, the simpler aspects of the theory of elliptic differential equations, 
and some elements of the theory of elasticity with particular emphasis on the 
torsion problem. 

The examples are, with one exception, worked out with the help of desk calcu- 
lators only. The question as to the suitability of this method for automatic com- 
puting machines is, unfortunately, not discussed. 


WOLFGANG Wasow 


Mathematics Research Center 
U. S. Army 

University of Wisconsin 
Madison, Wisconsin 


56[X].—J. KuntzMann, “Evaluations d’Erreur dans les Representations Ap- 
prochees de Derivees,”’ Societé d’ Electronique et d’Automatisme, 138 Boulevard 
de Verdun, Courbevoie (Seine), France, 1955, 49 p. (mimeographed), 30.5 cm. 


These tables serve to supplement an earlier report [1] by the same author. 
The error R,,¢,, in the numerical differentiation formulas 


f” (x¢) - 3 A ip.a.nVi + Ryans 
i=0 


is given by 
Roan =f Kyaldfo (oa 
where 
Kye(t) = ae Ain id Se (xi — #)° 1 if ¢ is between x; and xq, 


n! i 


~ 0 otherwise. 
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The upper bounds for the remainder terms 


P 5,408" Mays 
where 
| for (t) | << Mans 


and the quantities 
+= . 
Asan = / | Ky,¢.n(t) | dé 


are tabulated for m = 1(1)11, p, q = 0(1)n — 1. 


M. ABRAMOWITZ 
National Bureau of Standards 
Washington, D. C. 


1. J. Kuntzmann, “Formules de Dérivation Approchée au Moyen de Points Equidistants,” 
Report No. 1.373/1, Societé d’Electronique et d’Automatisme, 138 Boulevard de Verdun, Cour- 
bevoie (Seine), France, 1954. (MTAC, Rev. 51, v. 10, 1956, p. 171-172.) 
57[Z].—Preston C. Hammer, Editor, The Computing Laboratory in the Uni- 

versity, The University of Wisconsin Press, Wisconsin, 1957, xv + 236 p. 23 

cm. Price $6.50. 


A set of 31 papers dealing with various aspects of the computing laboratory in 
the University. These papers, each presumably presenting the individual views of 
its author, were presented at a conference held by the University of Wisconsin in 
August 1955. The papers vary in style and in content, but in total they convey an 
idea of many of the problems, many of the pleasures and some of the expenses of 
running a computing laboratory in a University. 

The University of Wisconsin should be thanked for sponsoring this conference 
and for making the material presented available in printed form. 

The table of contents follows: 

Computing and the University. The computing laboratory in the University by 
C. A. Elvehjem; The new significance of computation in higher education by J. H. 
Curtiss; Equipmental aids to computing by Jay W. Forrester. 

Applications of Computing in Science. Weather prediction by Philip Duncan 
Thompson; Computing in astronomy by W. J. Eckert; Applications of computing 
to fiuid dynamics problems by Harwood G. Kolsky; Applications of high-speed 
computing to chemical problems by Joseph O. Hirschfelder; The impact of fast 
computers on physics by Marshall Rosenbluth; The use of desk calculators by 
Paul S. Dwyer. 

Applications of Computing in Industry. The computer laboratory in industry 
by H. R. J. Grosch; Applications of computing in the Aircraft Industry by H. S. 
Wolanski; Computers improve power system performance by L. K. Kirchmayer; 
Assignment, programing, and scheduling by David F. Votaw, Jr. 

Future Demands for Personnel in Computing. Future demands for trained per- 
sonnel by E. K. Ritter; Supply and demand in computational mathematics by 
Forman S. Acton; The future demand for mathematicians in the computing field 
by R. E. Gaskell; Future demands for engineers and scientists in the field of comput- 
tation by Eldred C. Nelson. 
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Curriculums for Computing. The contribution of the computing laboratory to 
the University curriculum by Charles W. Adams; The educational program in 
Numerical Analysis of the Department of Mathematics, U.C.L.A. by George E. 
Forsythe; Curriculum needs in the computing field by Vincent C. Rideout; The 
Numerical Analysis Program at the University of Maryland by David M. Young, Jr. 

Equipping a University Computing Laboratory. Equipping the University com- 
putation laboratory by John W. Carr III; Equipping a University computing 
laboratory by C. C. Gotlieb; Equipping a University laboratory to satisfy the 
computational demands by H. O. Hartley; Equipping a University computing 
laboratory by Ralph E. Meagher; Equipping the university computing laboratory 
by Alan J. Perlis. 

Organizing and Financing a University Computing Laboratory. Organizing and 
financing a University computing laboratory by J. P. Nash; On organizing and 
financing a laboratory by Carl F. Kossack; The University computation laboratory 
as a cooperative industry-education project by Arvid W. Jacobson; The Cornell 
Computing Center, a minimum University computing laboratory by R. J. Walker. 

Computation in Perspective. Dangerous gulfs: some reflections on the social 
implications of computing machines by J. H. Van Vleck. 


ee 


58[Z].—Proceedings of the Third Annual Computer Applications Symposium, held 
October 9-10, 1956, sponsored by the Armour Research Foundation of the 
Illinois Institute of Technology, Chicago, Illinois, 1957, 148 p., 23 cm. Price 
$3.00. 


This conference had as its theme the application of medium scale digital com- 
puters. It is stated that ‘‘a primary purpose of the symposium was to illustrate 
diversity among computing machines and among applications, to provide a true 
picture of medium-scale computer applications as existing today.”” Following a de- 
scription of the medium scale computers to be considered and introductory articles, 
“The election and the UNIVAC” by C. Collingwood, and ‘“‘Model making problems 
in electric forecasting” by M. A. Woodbury, there are two sections of papers, one 
devoted to Business and Management Applications and the other to Engineering and 
Research Applications. The table of contents of these two sections follows: 
Business and Management Applications 

Development of a Products Pipe Line Simulator on an NCR 102A by J. H. 

Mallas 

Application of the IBM 650 to Stock Brokerage Operations by V. Lazzaro 

The Elecom 125 in Personnel Classification Research by J. M. Leiman 

Programming Ordinary Life Insurance Operations for the Datatron by J. S. Hill 

Manufacturing Data Processing on the IBM 650 by H. H. Marlow, Jr. 

The Univac File-Computer Applied to General Accounting Functions by 

T. R. Lyon 
Panel Discussion—Business and Management Applications 
Engineering and Research Applications 
Solution of Rotating Electric Machinery Problems with ALWAY by C. G. 
Veinott 
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The IBM 650 Applied to Problems of the Electrical Industry by R. Haberman, 
Jr. and F. J. Maginniss 
The NCR 102A as an Aid in Training and Research by E. J. Stewart 
Optical Calculations using the Burroughs E101 by A. Cox 
Use of the Datatron in the Petroleum Industry by J. S. Aronofsky 
Panel Discussion—Engineering and Research Applications 


Cc. B. T. 


59[Z].—R. Hunt Brown, Office Automation, Automation Consultants, Inc., New 
York, 1955, looseleaf, xviii + 283 p. Price $12.50. “Service to keep manual 
up-to-date” offered at $25.00 per year. 


This book starts with a readable section on the “what,” “why,” and “how,” 
of office automation, goes on to extensive discussions of hardware, especially 
electronic computers, briefly considers integrated-data-processing accounting sys- 
tems, and closes with short sections on sociological implications, scientific frame- 
works, and probable future developments related to office automation. The manner 
of presentation is aimed at the nontechnician, but computer experts not widely 
familiar with business applications will find much of interest in this book. 

The core of Office Automation seems to this reviewer to be the “Hardware Sec- 
tion” (p. 25-204), which indeed comprises over 60 per cent of the book’s total 
wordage, under chapter headings which can be condensed as follows: The common 
language concept; Integrated data processing; Electrical communication facilities; 
Common language machines; Native language machines; Tabulating and com- 
puting equipment; Punched card machine accounting; Fundamentals of electronic 
computers (four chapters); Principal electronic business computers; Memory 
systems; Output printers; Programming; Computer capabilities; Choosing a com- 
puter; Rent or buy; Miscellaneous topics (three chapters). The greatest value of 
this “Hardware Section” lies in its compact but expressive descriptions of specific 
equipment. If the promised up-dating service is well maintained, this book will 
continue to be a handy reference for everyone who is concerned with commercial 
computing equipment. 

One criticism, of a type which can only too often be levelled at semipopular 
treatments, is that too much of the material reads like a “sales-pitch.” For example, 
on p. 8, integrated data processing and electronic equipment is described as “the 
answer”’ (for the present purpose it matters not to what question). However, the 
“show-me”’ attitude of the typical business reader probably demands such a treat- 
ment, and almost surely guarantees that little or no harm will result from it. 

On the positive side, having discussed concepts of automation with several 
business administration students before and after they read Office Automation, the 
reviewer feels that it does get the main points across, and does so effectively and 
without undue effort on the reader’s part. 


James R. JAcKson 


University of California, 
Angeles 





TABLE ERRATA 


TABLE ERRATA 
The following erratum is mentioned in this issue: 


J. Peters, Ten-Place Logarithms of the Numbers from 1 to 100,000 Review 1 

page —. 

259.—MATHEMATICAL TABLES Project, New York, A. N. Lowan, technical dis 
rector, Tables of Lagrangian Interpolation Coefficients, Columbia University 
Press, New York, 1944 and 1948 editions. [MTAC, Rev. 162, v. 1, 1943-1945 
p. 314-315.] 


Page xxx in Introduction—IJn formulas (12), (13), (12a), and (13a), the fa 
1/n! should be deleted. 
GERTRUDE BLANCH 


214 W. Beechwood Avenue 
Dayton 5, Ohio 


260.—I. M. VinoGrapDov, Elements of Number Theory, Translated from the fift 
edition by Saut Kravetz, Dover Publications, Inc., New York, 1954. 


Four of the five errata in the English translation [1] occur in this edition. Th 
prime 2753 is listed correctly. In [1] the first erratum should read: 


p = 71. For ind 57 = 43, read ind 57 = 42. 
ANDRZEJ MAKOwsKI 


Warszawa 10 


ul. Bagatela 15 m. 49 
POLAND 


1. MTAC, Rev. 59, v. 9, 1955, p. 125. 


NOTES 
Research Potential and Training in the Mathematical Sciences 


The University of Chicago has published results of a survey of research poten 
tial and training in the mathematical sciences undertaken under grants to 
University by the National Science Foundation. The Chairman of the committe 
in charge was A. Adrian Albert, and the Investigator was John W. Green. Leon 
Cohen represented the National Science Foundation. The investigation include 
visits (and a few letters) to principal mathematics departments in American unit 
versities and colleges and questionnaires submitted to a large majority of th 
holders of the Ph.D. degree in mathematics. The report contains tabulations ¢ 
some of the more interesting information gathered covering the opinions of many 
mathematicians concerning the ease with which mathematical research can b 
carried out under present conditions and their opinions of changes which might b 
made to improve these conditions. 

Some undistributed copies of the report remain and may be available to li 
braries of properly interested organizations from the Program Director for Mathe 
matics, The National Science Foundation, Washington 25, D. C. 

C Bat 








